Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,493,437 A→T 23.7% intergenic (‑189/‑438) insG ← / → yjhB IS4 transposase/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,493,4370AT23.7% 581.0 / 192.6 286intergenic (‑189/‑438)insG/yjhBIS4 transposase/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (101/117);  new base T (27/41);  total (128/158)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 291 total aligned reads displayed.

CAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATG  >  CP009273/4493290‑4493584
                                                                                                                                                   |                                                                                                                                                   
caGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACTAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGa                                                                                                                                                     >  2:963198/1‑148 (MQ=255)
 aGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGAt                                                                                                                                                    <  1:755289/148‑1 (MQ=255)
       tGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGtc                                                                                                                                              >  1:3333510/1‑147 (MQ=255)
        ggATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGtat                                                                                                                                             <  1:2192885/148‑1 (MQ=255)
        ggATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGtat                                                                                                                                             <  1:2788826/148‑1 (MQ=255)
         gATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGtata                                                                                                                                            <  2:1454917/148‑1 (MQ=255)
         gATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGtata                                                                                                                                            <  2:2302496/148‑1 (MQ=255)
         gATCGGTCAACTGACCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGtata                                                                                                                                            <  2:3457242/148‑1 (MQ=255)
          aTCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATaa                                                                                                                                           >  1:329595/1‑148 (MQ=255)
           tCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAAc                                                                                                                                          >  1:1995632/1‑148 (MQ=255)
              gTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAgg                                                                                                                                       <  2:2186027/148‑1 (MQ=255)
                cAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTCTACCAGGta                                                                                                                                     >  2:1912528/1‑148 (MQ=255)
                cAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGta                                                                                                                                     >  2:504012/1‑148 (MQ=255)
                   gatctCCTTAACTTATCGGCATTAATCTTGGTTCTGGTTTTTGTAAAAAAATATCAGCTACAAAAATATTCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATaa                                                                                                                                  <  2:1939381/144‑1 (MQ=255)
                    tGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATaaa                                                                                                                                 <  2:1330361/148‑1 (MQ=255)
                       tCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGt                                                                                                                              <  1:1788693/148‑1 (MQ=255)
                              cTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGa                                                                                                                       <  1:1290051/148‑1 (MQ=255)
                              cTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGa                                                                                                                       <  1:1529168/148‑1 (MQ=255)
                               tGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGaa                                                                                                                      >  1:1018318/1‑148 (MQ=255)
                                   cGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGc                                                                                                                  <  1:3002164/148‑1 (MQ=255)
                                   cGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGc                                                                                                                  <  1:2878716/148‑1 (MQ=255)
                                      cATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCaaa                                                                                                               >  2:3028631/1‑148 (MQ=255)
                                       aTTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAg                                                                                                              >  2:3545015/1‑148 (MQ=255)
                                       aTTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAg                                                                                                              <  2:173084/148‑1 (MQ=255)
                                       aTTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAg                                                                                                              <  2:869404/148‑1 (MQ=255)
                                        ttAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAAAAGGTATAAAGGTATACAGAAAAGCAAAGa                                                                                                             >  2:1215364/1‑148 (MQ=255)
                                         tAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGaa                                                                                                            <  2:2278261/148‑1 (MQ=255)
                                           aTCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAAt                                                                                                          <  2:2004934/148‑1 (MQ=255)
                                            tCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATa                                                                                                         >  1:3094327/1‑148 (MQ=255)
                                            tCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGAAAAGAAATa                                                                                                         >  2:1805563/1‑148 (MQ=255)
                                             cTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAc                                                                                                        <  1:572901/148‑1 (MQ=255)
                                             cTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATAc                                                                                                        <  2:3020122/148‑1 (MQ=255)
                                              ttGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACt                                                                                                       <  2:3562642/148‑1 (MQ=255)
                                                  ttCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCaa                                                                                                   >  2:5498/1‑148 (MQ=255)
                                                           ttGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACa                                                                                          >  2:1419992/1‑148 (MQ=255)
                                                            tGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAg                                                                                         <  2:888197/148‑1 (MQ=255)
                                                              tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt                                                                                       <  2:1077429/148‑1 (MQ=255)
                                                              tAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCt                                                                                       <  2:2857961/148‑1 (MQ=255)
                                                               aaCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCta                                                                                      <  1:707166/148‑1 (MQ=255)
                                                                 cAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCtata                                                                                    <  1:907335/148‑1 (MQ=255)
                                                                 cAAAATATCAGCTACAAAAATATGATAAATTTGTGACTTCAGTAACAAAACGCGTTTTTTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCtata                                                                                    <  1:2693785/148‑1 (MQ=39)
                                                                   aaCTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAAc                                                                                  >  1:2598764/1‑148 (MQ=255)
                                                                    aCTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACg                                                                                 >  1:351732/1‑148 (MQ=255)
                                                                     cTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGt                                                                                >  1:3316250/1‑148 (MQ=255)
                                                                     cTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGt                                                                                <  1:1719690/148‑1 (MQ=255)
                                                                       aTCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTaa                                                                              <  2:2608802/148‑1 (MQ=255)
                                                                         cAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGc                                                                            >  2:941517/1‑148 (MQ=255)
                                                                           gCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTa                                                                          <  2:3214105/148‑1 (MQ=255)
                                                                           gCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTa                                                                          >  1:2768791/1‑148 (MQ=255)
                                                                             tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTaaa                                                                        >  1:1874638/1‑148 (MQ=255)
                                                                             tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTaaa                                                                        <  1:782236/148‑1 (MQ=255)
                                                                             tACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTaaa                                                                        <  2:754612/148‑1 (MQ=255)
                                                                              aCAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAg                                                                       >  1:997899/1‑148 (MQ=255)
                                                                                aaaaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGta                                                                     >  2:2112246/1‑148 (MQ=255)
                                                                                  aaaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAATCTAAAGtaat                                                                   >  2:2989055/1‑148 (MQ=255)
                                                                                   aaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaata                                                                  <  1:2501801/148‑1 (MQ=255)
                                                                                   aaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaata                                                                  <  1:2249399/148‑1 (MQ=255)
                                                                                   aaTATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaata                                                                  <  2:2676800/148‑1 (MQ=255)
                                                                                    atatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaataa                                                                 >  1:3114835/1‑148 (MQ=255)
                                                                                     tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAAc                                                                >  1:1426292/1‑148 (MQ=255)
                                                                                        gCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctc                                                             <  1:1005266/148‑1 (MQ=255)
                                                                                        gCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctc                                                             <  1:2982496/148‑1 (MQ=255)
                                                                                         cTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctct                                                            >  1:2384669/1‑148 (MQ=255)
                                                                                         cTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctct                                                            <  2:3367244/148‑1 (MQ=255)
                                                                                           cAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCa                                                          >  2:801463/1‑148 (MQ=255)
                                                                                              tttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTc                                                       >  1:1169203/1‑148 (MQ=255)
                                                                                              tttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTc                                                       <  2:1300004/148‑1 (MQ=255)
                                                                                               ttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCt                                                      <  2:2383743/148‑1 (MQ=255)
                                                                                                 gtgACATCAGTAACAAAACTCGTTTTGTTATGTGTATTGCTTGTTTTTCTATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCttt                                                    <  2:270319/148‑1 (MQ=255)
                                                                                                   gACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTcc                                                  >  2:3117927/1‑148 (MQ=255)
                                                                                                      aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCa                                               <  1:504012/148‑1 (MQ=255)
                                                                                                       tCAGTAACATAACGCTTTTTGTTATGTGGATTGCTTGTTTTTGGATCTGGCTTAACAGGTATATATGTATTCTGATAAGCAAAGAAATACTGAAAAGGAAAACAGCTTTTTAGTAAGTTATAGTAATCACCTCTATGTTTTTCCtatt                                              <  2:187344/148‑4 (MQ=255)
                                                                                                        cAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCAtt                                             <  2:1618083/148‑1 (MQ=255)
                                                                                                          gTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTg                                           <  2:2462117/148‑1 (MQ=255)
                                                                                                              cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGa                                       >  1:2066624/1‑148 (MQ=255)
                                                                                                               aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGaa                                      <  2:1116755/148‑1 (MQ=255)
                                                                                                               aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGaa                                      <  1:1934517/148‑1 (MQ=255)
                                                                                                               aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTCGACGaa                                      <  1:2436651/148‑1 (MQ=255)
                                                                                                               aaaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGaa                                      >  1:531342/1‑148 (MQ=255)
                                                                                                                aaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAg                                     <  1:3186478/148‑1 (MQ=255)
                                                                                                                 aaCGCGTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAgg                                    >  2:158706/1‑148 (MQ=255)
                                                                                                                   cgcgTTTTGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGa                                  >  1:1075800/1‑148 (MQ=255)
                                                                                                                         ttGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTa                            <  2:1690478/148‑1 (MQ=255)
                                                                                                                          tGTTATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTAt                           <  1:1893314/148‑1 (MQ=255)
                                                                                                                            ttATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTc                         >  1:1316808/1‑148 (MQ=255)
                                                                                                                             tATGTGGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCa                        >  2:689407/1‑148 (MQ=255)
                                                                                                                                gtgGATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAAcc                     <  1:8445/148‑1 (MQ=255)
                                                                                                                                 tgGATTGCTTGTTTTTCGTTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCt                    >  2:333233/1‑148 (MQ=255)
                                                                                                                                   gATTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGa                  <  2:2751760/148‑1 (MQ=255)
                                                                                                                                    aTTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCACAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGaa                 <  2:2088990/148‑1 (MQ=255)
                                                                                                                                    aTTGCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGaa                 <  1:3028631/148‑1 (MQ=255)
                                                                                                                                       gCTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACgg              <  1:1470100/148‑1 (MQ=255)
                                                                                                                                        cTTGTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGa             <  1:3402218/148‑1 (MQ=255)
                                                                                                                                           ttttttcGATCTGGTATACCAGGTATAAAGGTATAAATAAAAGCAAAGAAATACTGCAAAGGAAAACAGGTATAACTTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTa          <  1:1215364/147‑1 (MQ=255)
                                                                                                                                           gTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTa          <  1:1755321/148‑1 (MQ=255)
                                                                                                                                           gTTTTTCGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTa          <  2:2735590/148‑1 (MQ=255)
                                                                                                                                            tttttcGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTAc         >  1:1358162/1‑148 (MQ=255)
                                                                                                                                               ttcgATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGaa      <  2:2405072/145‑1 (MQ=255)
                                                                                                                                                tcgATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGaaa     <  2:402914/146‑1 (MQ=255)
                                                                                                                                                   aTCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATg  <  1:2870183/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCCACTGGATCGGTCAACTGATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATG  >  CP009273/4493290‑4493584

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