Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,775,093 G→A 5.3% intergenic (+72/‑452) yfjW → / → ypjI CP4‑57 prophage; putative inner membrane protein/pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,775,0930GA 5.3% 596.1 / 13.9 182intergenic (+72/‑452)yfjW/ypjICP4‑57 prophage; putative inner membrane protein/pseudogene, CP4‑57 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (90/82);  new base A (5/5);  total (95/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 182 total aligned reads displayed.

AGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAAGCGCCAG  >  NC_000913/2774946‑2775240
                                                                                                                                                    |                                                                                                                                                    
aGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGc                                                                                                                                                     >  2:1022652/1‑149 (MQ=255)
      gCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCtt                                                                                                                                               <  3:330994/149‑1 (MQ=255)
        tAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGa                                                                                                                                             >  1:532810/1‑149 (MQ=255)
        tAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGa                                                                                                                                             <  1:498233/149‑1 (MQ=255)
        tAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGa                                                                                                                                             <  3:888246/149‑1 (MQ=255)
          gCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGc                                                                                                                                           >  2:709328/1‑149 (MQ=255)
             aTTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGCTTTAAAAACAGGCttaacttag                                                                                                                                        <  1:876093/149‑8 (MQ=255)
                 cTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGc                                                                                                                                    <  2:337246/149‑1 (MQ=255)
                      aaaTTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGaca                                                                                                                               >  1:32365/1‑149 (MQ=255)
                       aaTTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGacac                                                                                                                              >  1:816400/1‑149 (MQ=255)
                        attcattcattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACt                                                                                                                             >  3:1157648/1‑149 (MQ=255)
                           acatcattcTTAAAGTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCAt                                                                                                                          >  1:180478/4‑149 (MQ=255)
                             ttcattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATAc                                                                                                                        >  3:746401/1‑149 (MQ=255)
                             ttcattcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATAc                                                                                                                        >  4:596648/1‑149 (MQ=255)
                                 ttcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGc                                                                                                                    >  1:1233913/1‑149 (MQ=255)
                                 ttcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGc                                                                                                                    <  3:1490312/149‑1 (MQ=255)
                                 ttcTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGCTTTAAAAACAgg                                                                                                                                                  <  4:1007640/119‑1 (MQ=255)
                                       cTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTaaa                                                                                                              <  2:1382580/149‑1 (MQ=255)
                                          tAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                                                                           <  3:1418250/149‑1 (MQ=255)
                                          tAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTg                                                                                                           <  1:1086647/149‑1 (MQ=255)
                                             gAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCg                                                                                                        >  1:679260/1‑149 (MQ=255)
                                             gAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCg                                                                                                        >  1:816453/1‑149 (MQ=255)
                                               atgatgTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGtt                                                                                                      <  3:424685/149‑1 (MQ=255)
                                                  atgTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGttttt                                                                                                   >  2:687284/1‑149 (MQ=255)
                                                    gtgAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGt                                                                                                 <  2:254060/149‑1 (MQ=255)
                                                     tgaatgaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGaca                                                                                                                               >  3:588746/1‑118 (MQ=255)
                                                     tgaatgaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGaca                                                                                                                               <  4:588746/118‑1 (MQ=255)
                                                     tgaatgaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTc                                                                                                >  4:140026/1‑149 (MQ=255)
                                                         tgaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTa                                                                                            <  2:958338/149‑1 (MQ=255)
                                                         tgaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTa                                                                                            <  2:1435688/149‑1 (MQ=255)
                                                          gaaAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTAc                                                                                           <  2:1158774/149‑1 (MQ=255)
                                                           aaaaTTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCTTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTAcc                                                                                          >  2:757044/1‑149 (MQ=255)
                                                           aaaaTTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTAcc                                                                                          <  1:320863/149‑1 (MQ=255)
                                                            aaaTTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTAccc                                                                                         <  2:654564/149‑1 (MQ=255)
                                                              aTTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCAt                                                                                       >  4:1532152/1‑149 (MQ=255)
                                                               tttATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGa                                                                                                                                             >  4:952927/1‑94 (MQ=255)
                                                               tttATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGa                                                                                                                                             <  3:952927/94‑1 (MQ=255)
                                                               tttATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGa                                                                                                                                 >  4:298732/1‑106 (MQ=255)
                                                               tttATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATAc                                                                                                                        >  1:690141/1‑115 (MQ=255)
                                                                      ttGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTatat                                                                              >  3:28504/1‑149 (MQ=255)
                                                                          aaCGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACTCATA‑CCTATATTTTc                                                                          <  4:1157648/149‑1 (MQ=255)
                                                                          aaCGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTc                                                                          <  4:28504/149‑1 (MQ=255)
                                                                          aaCGCATTTGGCTCCAAGTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTc                                                                          <  4:1325020/149‑1 (MQ=255)
                                                                           aCGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCt                                                                         >  3:774876/1‑149 (MQ=255)
                                                                             taaaaaggCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCtgt                                                                       >  1:1320683/7‑149 (MQ=255)
                                                                              cATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCtgtt                                                                      <  2:1343943/149‑1 (MQ=255)
                                                                               aTTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCtgttg                                                                     >  1:1488390/1‑149 (MQ=255)
                                                                                tttGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCtgttgt                                                                    <  2:180478/149‑1 (MQ=255)
                                                                                 aaggCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTgttgtt                                                                   >  3:320819/3‑149 (MQ=255)
                                                                                   ggCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTaaa                                                                                                              >  1:930410/1‑105 (MQ=255)
                                                                                   ggCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTgttgttca                                                                 >  2:461360/1‑147 (MQ=255)
                                                                                    gCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTaa                                                                <  2:679260/149‑1 (MQ=255)
                                                                                      tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGa                                                                                                                                 >  3:1485164/1‑83 (MQ=255)
                                                                                      tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCt                                                                         >  1:1014578/1‑138 (MQ=255)
                                                                                      tCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCt                                                                         <  2:1014578/138‑1 (MQ=255)
                                                                                          aTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCa                                                          <  4:428142/149‑1 (MQ=255)
                                                                                          aTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCa                                                          >  2:1482723/1‑149 (MQ=255)
                                                                                          aTTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTCGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCa                                                          >  4:14995/1‑149 (MQ=255)
                                                                                               aGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACt                                                     >  1:309646/1‑149 (MQ=255)
                                                                                                tctttttttTGACAATTT‑TTTAAATATTCTTAGTTAAAGCCTGTTTTTAAGGCGGGCTGGGCGGTGAAGCGACCACTAAT‑CCCGCCTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGa                                                   <  4:1477377/146‑1 (MQ=255)
                                                                                                gCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTg                                                    >  3:754524/1‑149 (MQ=255)
                                                                                                  cTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAGGT‑GCGGTTTTGGCCTTAACCCTACCCAATATTTTTTTTGGTTAACCTCCCGGAACtaga                                                  >  3:1477377/1‑146 (MQ=255)
                                                                                                   tttttATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGc                                                 >  1:1167243/1‑149 (MQ=255)
                                                                                                   tttttATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGc                                                 >  2:537638/1‑149 (MQ=255)
                                                                                                      ttATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGt                                              <  4:1482694/149‑1 (MQ=255)
                                                                                                       tATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                             <  3:1008235/149‑1 (MQ=255)
                                                                                                       tATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTGCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                             <  1:757044/149‑1 (MQ=255)
                                                                                                          gACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaa                                          >  4:693112/1‑149 (MQ=255)
                                                                                                           aCTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaa                                         <  2:108669/149‑1 (MQ=255)
                                                                                                           aCTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaa                                         >  1:714689/1‑149 (MQ=255)
                                                                                                            cTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaaa                                        <  2:803958/149‑1 (MQ=255)
                                                                                                            cTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGaaaa                                        <  1:461360/149‑1 (MQ=255)
                                                                                                             tatttattAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGaca                                                                                                                               <  1:664624/62‑1 (MQ=255)
                                                                                                             tatttattAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGGTGAAAAg                                       >  1:1562818/1‑149 (MQ=255)
                                                                                                                 tattaAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAg                                   >  3:1208558/1‑149 (MQ=255)
                                                                                                                  attaAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGa                                  <  2:309646/149‑1 (MQ=255)
                                                                                                                   ttaAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCtgtt                                                                      >  2:13010/1‑112 (MQ=255)
                                                                                                                   ttaAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCtgtt                                                                      <  1:13010/112‑1 (MQ=255)
                                                                                                                      aaaTATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAata                              >  4:645349/1‑149 (MQ=255)
                                                                                                                       aaTATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAatat                             <  2:908999/149‑1 (MQ=255)
                                                                                                                         tatCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGc                           <  3:261475/149‑1 (MQ=255)
                                                                                                                             cTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAg                       >  1:982717/1‑149 (MQ=255)
                                                                                                                              tAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGa                                                   >  4:602429/1‑120 (MQ=255)
                                                                                                                              tAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGc                                      <  3:1535996/133‑1 (MQ=255)
                                                                                                                              tAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGc                      >  3:993698/1‑149 (MQ=255)
                                                                                                                              tAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGc                      >  4:1557673/1‑149 (MQ=255)
                                                                                                                                   tcaaGCCTGCTTTAAAAACAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCa                                    <  1:1105941/129‑1 (MQ=255)
                                                                                                                                   tcaaGCCTGCTTTAAAAACAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCa                                    >  4:760602/3‑131 (MQ=255)
                                                                                                                                       gCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGa                        >  3:563585/1‑139 (MQ=255)
                                                                                                                                       gCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGa                        <  4:563585/139‑1 (MQ=255)
                                                                                                                                        ccTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCAt             <  4:1305914/149‑1 (MQ=255)
                                                                                                                                        ccTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCAt             <  2:714689/149‑1 (MQ=255)
                                                                                                                                        ccTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCAt             <  2:816453/149‑1 (MQ=255)
                                                                                                                                         cTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGa                        <  4:450328/137‑1 (MQ=255)
                                                                                                                                         cTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGa                        >  3:450328/1‑137 (MQ=255)
                                                                                                                                            tttttAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATACCGAGCCAGCGGCATGGaa         <  2:360615/149‑1 (MQ=255)
                                                                                                                                              tttAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAAgc       <  2:1233913/149‑1 (MQ=255)
                                                                                                                                                tAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                             >  2:607586/1‑109 (MQ=255)
                                                                                                                                                tAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGtt                                             <  1:607586/109‑1 (MQ=255)
                                                                                                                                                tAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAAgcgc     <  3:140026/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAATGAAAATTTATTCTTGTAACGCATTTGGCTCCAATTGGAGCCTTTTTATGACTATTTATTAAAATATCCTAAGTT‑AAGCCTGTTTTTAAAGCAGGCTTGAGCGTGAAGCGGACACTCATACCTGCTTTAAAGTGGCGTTTTTGTCGTTACCCATA‑CCTATATTTTCTGTTGTTTAACTGCCAGAACTGAGCGGTTGAAAAGCCAGAAATATGCCGAGCCAGCGGCATGGAAGCGCCAG  >  NC_000913/2774946‑2775240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑