Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,630,704 C→G 15.2% intergenic (‑80/‑309) nohQ ← / → ynfO pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,630,7040CG15.2% 310.6 / 14.0 144intergenic (‑80/‑309)nohQ/ynfOpseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (23/79);  new base G (19/2);  total (63/81)
Fisher's exact test for biased strand distribution p-value = 6.02e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

TGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGT  >  CP009273/1630557‑1630850
                                                                                                                                                    |                                                                                                                                                     
tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttct                                                                                                                                                      <  2:4178474/149‑1 (MQ=255)
 gCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGACACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttctt                                                                                                                                                     <  2:1330551/149‑1 (MQ=255)
    aGTTCTTATTGGTACGGATACTCGCACCTAAAAGGTCAGCCAGCTGCTTTTTGTTGACTGCCATTGTTAATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTc                                                                                                                                                  <  1:4507644/149‑1 (MQ=14)
     gTTCTGAATGGTACTGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGCCCAAAA‑ACACAGCACGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCa                                                                                                                                                 <  2:1265562/149‑1 (MQ=14)
     gTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGGTTTAttttga                                                                                                                                                 >  1:1729353/1‑147 (MQ=255)
     gTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGGTTTAttttga                                                                                                                                                 >  1:1729609/1‑147 (MQ=255)
     gTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTCttttta                                                                                                                                                 >  2:3990250/1‑147 (MQ=255)
     gTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCa                                                                                                                                                 <  2:568916/149‑1 (MQ=255)
     gTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCa                                                                                                                                                 <  2:3781596/149‑1 (MQ=255)
      ttCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTCCATTCATTTCAg                                                                                                                                                >  2:1935186/1‑149 (MQ=255)
         tGAATGGTACGGATACTCGCACCTAAATTGTCAGCCAGCTGATTTTTGTTGACTTCAATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAgggg                                                                                                                                             <  1:3755702/149‑1 (MQ=255)
          gAATGGTACGGATACTCGCACCGAAAATGTCCGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAggggg                                                                                                                                            <  2:3197736/149‑1 (MQ=255)
           aaTGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCTGTATTTTAAgggggg                                                                                                                                           >  1:2943516/1‑148 (MQ=14)
            aTGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTa                                                                                                                                          <  2:4138715/149‑1 (MQ=255)
            aTGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTa                                                                                                                                          <  1:2655139/149‑1 (MQ=255)
             tGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTAt                                                                                                                                         <  2:2676180/149‑1 (MQ=255)
             tGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTAt                                                                                                                                         <  2:2676143/149‑1 (MQ=255)
              gttactcaTACTCGCACAGAAATTGTCAGCCTGATTTTTTTTGTTGACTTCCATGTTACATCCCCCGCCCAAAA‑ACAGATAAAGCAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTAtt                                                                                                                                        <  1:4266849/142‑1 (MQ=11)
                 aCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTTTACATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTta                                                                                                                                     <  2:2041252/149‑1 (MQ=255)
                    gATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCGTTGTTTTTAGGGGGGATTTttaat                                                                                                                                  >  1:2582622/1‑149 (MQ=14)
                      tACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTtaaataa                                                                                                                                <  2:213554/149‑1 (MQ=255)
                         tCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTGTTTTAAGGGGGTATTTTAAATaaaaa                                                                                                                             >  2:936752/1‑149 (MQ=255)
                         tCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGGTTTTTTTTAATGGGGGGTTTTTAATaaaaa                                                                                                                             >  1:3867820/1‑149 (MQ=14)
                         tCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGGTTTGTTTTGTTGGGGGTTTTTTTATaaaaa                                                                                                                             >  1:412243/1‑149 (MQ=14)
                         tCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGGTTTATTTTTAATGGGGATTTTATATaaaaa                                                                                                                             >  1:3868138/1‑149 (MQ=14)
                          cGCACCTAAATTGTCAGCCAGCTGCTTTTTGTTGACTTCATTTGTTCATTCCACGGCCAAAA‑ACAGTGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAc                                                                                                                            <  2:16866/149‑1 (MQ=14)
                          cGCACCTAAATTGCCAGCCAGCTGCTTTTTGTTGACTTCCTTTGTTCATCCCACGGCCAAAA‑ACAGATAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAc                                                                                                                            <  2:16868/149‑1 (MQ=14)
                            cACCTAAAATGTCAGCCTGTTGCTTTTTGTTCACTGCCATTGTTCATACCACGCCAAAAA‑ACAGATAAAGCACACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACAt                                                                                                                          <  1:4339873/149‑1 (MQ=14)
                            cACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTGTTTTGAGGGGGTATTTTTAATaaaaaaaa                                                                                                                          >  2:1578709/1‑146 (MQ=14)
                            cACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACAt                                                                                                                          <  2:741496/149‑1 (MQ=255)
                            cACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTAAATTCATTTAAGGGGGGATTTTAAATaaaaaaaa                                                                                                                          >  2:1898082/1‑146 (MQ=14)
                             aCCTAAATTGTCAGCCTGCTGCTTTTTGTTGAATTCAATTTTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACAtt                                                                                                                         <  1:3712586/149‑1 (MQ=14)
                              cctaaaaTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTa                                                                                                                        <  2:1827827/146‑1 (MQ=255)
                              ccGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAATGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTa                                                                                                                        <  2:1981504/149‑1 (MQ=255)
                              ccGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTa                                                                                                                        <  2:1828068/149‑1 (MQ=255)
                               cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCTTTGTTCATTCCCCGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                                                       <  2:3227362/149‑1 (MQ=255)
                               cGAAAATGTCAGCCAGCTGATTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                                                       <  1:2463440/149‑1 (MQ=255)
                                 aaaaTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTATTTTTATGGGGGTTTTTATAAAAAAAAATTaatt                                                                                                                     >  1:4076300/1‑147 (MQ=14)
                                   aaTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTATTTTCAGGGGGTATTTTAAATAAAAACAATAAGttt                                                                                                                   >  2:3693144/1‑148 (MQ=255)
                                    ttGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCAGTCCACGGCCAAAA‑ACAGAGAATGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTac                                                                                                                  <  2:1396654/148‑1 (MQ=255)
                                     tGTCAGCCATCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacg                                                                                                                 <  2:1455349/149‑1 (MQ=255)
                                     tGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAGACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCGTTGTTTTTGGGGGGTTTTTTTTATAAAAAAAATAAGTTaag                                                                                                                 >  1:571346/1‑147 (MQ=14)
                                     tGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacg                                                                                                                 <  2:330756/149‑1 (MQ=255)
                                      gTCAGCCTGCTGCTTTTTGTTCACTTCAATTGTACATTCCACGCCAAAAA‑ACAGATAACGCACACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacg‑a                                                                                                               <  1:99303/149‑1 (MQ=14)
                                       tccgccagcTGCTTTTTTTTGACTTCATTTTTTCATTCCACGGAAAAAA‑ACATACAAAGCAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacg‑ac                                                                                                              <  1:881447/146‑1 (MQ=14)
                                       tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacg‑ac                                                                                                              <  2:1094519/149‑1 (MQ=255)
                                          gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGGTTTATTTTTAGGGGGTATTTTTAATAAAAAAAATAATTTaataagaa                                                                                                            >  1:2622586/1‑142 (MQ=14)
                                          gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCGTTGTTTTAAGGGGGTTTTTTTAATAAAAAAATATATTTAAG‑ACgaa                                                                                                           >  1:331813/1‑149 (MQ=14)
                                          gccagcTGATTTTTGTTGAATTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACgaa                                                                                                           <  2:3928944/149‑1 (MQ=255)
                                              gctgctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACgaagaag                                                                                                       <  2:272061/149‑1 (MQ=255)
                                                tgctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTATTTTAACGGGGTATTTTAAATAAAAAAAATAAATTAAG‑ACcaataat‑aa                                                                                                    >  2:571442/1‑149 (MQ=14)
                                                tgctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCGTTCATTTAAGGGGGGATTTTAAAAAAAAACAGTAAATTACC‑ACgaagaag‑aa                                                                                                    >  1:2076889/1‑149 (MQ=14)
                                                 gtttttttttGACTTCCATTGTTCATTCCACGCCAAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑c                                                                                                  <  1:2900425/147‑1 (MQ=255)
                                                 gctTTTTTTTGACTTCCATTGTTCATACCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑c                                                                                                  <  1:2901495/149‑1 (MQ=255)
                                                  ctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTATTTTCAGGGGGTATTTTAAATAAAAAAAATAAATTAAG‑AACCAGAAG‑AA‑cc                                                                                                 >  2:1591758/1‑148 (MQ=14)
                                                        gttgACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCGTTCTTTTTAGGGGGGATTTTAAATAAAAAAAATAAATTAAGAACAAAAACA‑AA‑CGGAAAt                                                                                            >  2:1272515/1‑149 (MQ=14)
                                                        gttgACTTCCATTGTTCATTCCACGCCAAAAA‑ACAGATAAAGCAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATg                                                                                           <  1:2484429/149‑1 (MQ=255)
                                                        gttcACTTCCATTGTACATTCCACGCCAAAAA‑ACAGATAAAGTAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATg                                                                                           <  1:680182/149‑1 (MQ=14)
                                                             cTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTGAGGGGGTATTTTAAATAAAAACATTAATTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTa                                                                                      >  2:1125118/1‑149 (MQ=255)
                                                              ttCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCTTTCTTTTAAGGGGGGATTTTAAATAAAAACATTAATTTACG‑AAGAAGAAG‑AA‑AGGAAATGCCttta                                                                                     >  2:219184/1‑147 (MQ=17)
                                                               tcaaTTTTTCATACCCCGCCCAAAA‑ACAGATAATGTACACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTaaa                                                                                    <  1:959058/146‑1 (MQ=14)
                                                                  aTTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCTTTCGGTTTAGGGGGGATTTTAAATAAAAAAAAAAAGTTACG‑ACCAAAAAG‑AACCGAAAAGGCCTAAAAcc                                                                                  >  1:3380942/1‑149 (MQ=14)
                                                                  aTTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCTTTATTTTTAGGGGGGATTTTTAAAAAAAAAAATAAATTAAG‑ACGAAGAAGAAA‑CGAAAATACATAAAAcc                                                                                  >  1:1071091/1‑149 (MQ=14)
                                                                   ttGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTCCTTTCTTTTGTGGGGGTATTTTTAATAAAAAAAATAAGTTACG‑ACGGAGAAG‑AA‑CGGGAATGCGTTAAAcccg                                                                                >  2:2349727/1‑147 (MQ=17)
                                                                      ttCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTATTTTGATGGGGTATTTTAAAAAAAAAAATTAATTTAAG‑AAGAAGAAG‑TA‑ACGAAATGCCTGAAAACCGaaa                                                                             >  2:4351156/1‑149 (MQ=14)
                                                                       tcagaccACGCCCAAAA‑ACAGATAACGCAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGaaaa                                                                            <  1:2449513/144‑1 (MQ=14)
                                                                       tCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTCCGTTCTTTTCAGGGGGGATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAAAGCCGTAAACCCGcaaa                                                                            >  2:382279/1‑149 (MQ=14)
                                                                       tCATCCCCCGGCAAAAA‑ACAGAGAACTGACACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCGGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGaaaa                                                                            <  1:158807/149‑1 (MQ=14)
                                                                         aTTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCGTTGTTTTTAGGGGGTGTTTTAAATAAAAAAAATAATTTAAA‑AAGAAGAAG‑AA‑ACGAAAAGCCATAAAACCGaaaaat                                                                          >  1:3313986/1‑147 (MQ=14)
                                                                         aTTCCACGCCAAAAA‑ACAGACAAAGCACACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAAtt                                                                          <  1:1596956/149‑1 (MQ=14)
                                                                            cccccgCCAAAACACAGAGACCGCAAACGACAGAGGCCCAACAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAAtttt                                                                        <  1:3484434/144‑1 (MQ=14)
                                                                               cGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATa                                                                    <  2:3900891/149‑1 (MQ=255)
                                                                                   cAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGGTTTGGTTTAAGGGGGTTTTTTTAAAAAAAAAAATAATATAAG‑AACAATAAA‑AA‑AAGAAAAACCATAAAACAGAAAAATTTAATaaaaa                                                                >  1:4418980/1‑147 (MQ=2)
                                                                                                aaGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGTGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCcgc                                                    <  1:993565/149‑1 (MQ=255)
                                                                                                    aaaCTAAAGAGGCCAAAAAGCTGGTTTTCAGAACCTGTCGTTTCCTTTGTTTTCAGGGGGTATTTTAAATAAAAAAATTATTTTACG‑ACGAAGATG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAgg                                                <  1:1236524/149‑1 (MQ=11)
                                                                                                            gagGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTATTTTGAGGGGGTATTTTAAATAAAAAAAATAAGTTAAG‑ACCAAGAAG‑AA‑CGGAAAATCCTGAAACCCGAAAAATTTAAAAAATATCGAAAAACCCCGAGGTCgccgcc                                        >  2:1754132/1‑149 (MQ=11)
                                                                                                            gagGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTACTTTATTTTAATGGGGTATTTTAAATAAAAAAAATTAGTTACG‑AAGAAGAAA‑AA‑ACGAAATGGCTTAAAACGGAAAAATTTCATAAATAGAGAAAACCCGCGGGGTCgccccc                                        >  2:2405459/1‑149 (MQ=255)
                                                                                                             agGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGccc                                       <  2:858640/149‑1 (MQ=255)
                                                                                                              ggCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTCCTTTATTTTAAGGGGGTATTTTAAATAAAAAAAATAAATTAAG‑ACCAAGAAG‑AA‑AGGAAAAGCCTTAAACCCGCAAATATTCATAAAAAACGAAAAACCCCCCGGTCGCcccccc                                      >  2:4396928/1‑149 (MQ=9)
                                                                                                                cccAAAAGCTCGTTTTCAGCACCTGTCGTGTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATACCGAAAACCCGCGAGGTCGCCGCCCCGt                                    <  2:419357/149‑1 (MQ=255)
                                                                                                                  cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa                                  <  2:2434122/149‑1 (MQ=40)
                                                                                                                  cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa                                  <  2:2434060/149‑1 (MQ=40)
                                                                                                                     aaGCTCGTTTTCAGCACCTGTCGTTTGCATTCTTTTGAGGGGGGATTTTAAATAAAAAAAATAAATTAAC‑ACGGAGAAG‑AACCGGAAAT‑CCGTCAACCCGAAAAATTTCAAAAAAAGCCAAAACCCCCCCGGTCGCCCCCCCCGAAcca                               >  2:2091333/1‑148 (MQ=255)
                                                                                                                         tCGTTTTCAGCACCTGTCGTGTCCGTTGTTTTCAGGGGGTATTTTAAATAAAAAAAATAAGTTACG‑ACGAACAAG‑AA‑CGGAAAAGCCATAAACCCGCCAAATTTCATAAATAGCGCAAACCCCCCAGGGCGCCCCCCCCCAAccagacc                           >  2:3564665/1‑144 (MQ=2)
                                                                                                                            ttttCTGAAACTTTCGTTTACTTTGTTTTCTGGGTTTATTTTAAAAATTTAGTTTTAGTTAGT‑AGGATGAAG‑AA‑CGGAAATGACATAAACCGTAAAATTTTCATAAATCGCGAAAACCCGCGAGGTCGCCGCCCCGTTACCTGTCGGAt                        <  1:40576/149‑1 (MQ=2)
                                                                                                                            ttttCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt                        <  2:3421806/149‑1 (MQ=38)
                                                                                                                            ttttCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGAt                        <  2:3666624/149‑1 (MQ=38)
                                                                                                                               tCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCAc                     <  1:3345470/149‑1 (MQ=38)
                                                                                                                                      cTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg              <  1:198660/149‑1 (MQ=39)
                                                                                                                                      cTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg              <  2:4364230/149‑1 (MQ=39)
                                                                                                                                            tttCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCg        >  1:4513353/1‑149 (MQ=40)
                                                                                                                                            tttCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCg        >  1:4513249/1‑149 (MQ=40)
                                                                                                                                            tttCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCg        <  2:3345099/149‑1 (MQ=40)
                                                                                                                                               ccTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaa     >  1:2749511/1‑149 (MQ=255)
                                                                                                                                                ctttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaaa    >  1:835281/1‑149 (MQ=255)
                                                                                                                                                 tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAg   >  2:1720973/1‑149 (MQ=255)
                                                                                                                                                 tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAg   >  2:1721010/1‑149 (MQ=255)
                                                                                                                                                 tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAg   >  1:4029959/1‑149 (MQ=255)
                                                                                                                                                  ttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt  <  2:412243/149‑1 (MQ=255)
                                                                                                                                                  ttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt  >  1:4551449/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAA‑ACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACG‑ACGAAGAAG‑AA‑CGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGT  >  CP009273/1630557‑1630850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑