Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,568,632 A→T 11.4% L62Q (CTA→CAA)  eutP ← putative P‑loop NTPase ethanolamine utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,568,6320AT11.4% 549.5 / 51.8 189L62Q (CTA→CAA) eutPputative P‑loop NTPase ethanolamine utilization protein
Reads supporting (aligned to +/- strand):  ref base A (80/87);  new base T (7/14);  total (88/101)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 190 total aligned reads displayed.

TTGCGGCGACGTCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTC  >  CP009273/2568486‑2568779
                                                                                                                                                  |                                                                                                                                                   
ttGCGGCGACGTCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAg                                                                                                                                                    >  2:1005041/1‑148 (MQ=255)
      cgACGTCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGt                                                                                                                                              >  2:1721658/1‑148 (MQ=255)
       gACGTCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCACCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTa                                                                                                                                             <  2:1341125/148‑1 (MQ=255)
         cGTCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATCTCTACATCCTGTAGCGTGGTAAt                                                                                                                                           >  2:1878886/1‑148 (MQ=255)
           tCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTa                                                                                                                                         >  1:2577501/1‑148 (MQ=255)
            cAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTaa                                                                                                                                        <  1:313755/148‑1 (MQ=255)
            cAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTaa                                                                                                                                        <  1:1355583/148‑1 (MQ=255)
            cAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTaa                                                                                                                                        <  1:1722712/148‑1 (MQ=255)
            cAGCATCTGGCATGTCAGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCCTGGTAATTaa                                                                                                                                        >  2:1739846/1‑148 (MQ=255)
              gCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAgg                                                                                                                                      >  1:1088925/1‑148 (MQ=255)
               cATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGc                                                                                                                                     >  2:141729/1‑148 (MQ=255)
                aTCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCg                                                                                                                                    <  2:2363170/148‑1 (MQ=255)
                aTCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCg                                                                                                                                    <  2:2363148/148‑1 (MQ=255)
                  cTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTg                                                                                                                                  <  1:1807129/148‑1 (MQ=255)
                    ggCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGAt                                                                                                                                >  2:2254823/1‑148 (MQ=255)
                      cATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATAc                                                                                                                              <  2:2603774/148‑1 (MQ=255)
                         gTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAg                                                                                                                           >  2:2698784/1‑148 (MQ=255)
                          tCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGc                                                                                                                          <  2:2630662/148‑1 (MQ=255)
                           ccGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGGCATTTGCGCCATGGACATAGATCAGCCTGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCg                                                                                                                         >  1:2782349/1‑148 (MQ=255)
                            cGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCgg                                                                                                                        <  2:883822/148‑1 (MQ=255)
                              ttttACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCgggg                                                                                                                      >  2:1084949/1‑148 (MQ=255)
                               tttACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGa                                                                                                                     >  2:1246497/1‑148 (MQ=255)
                                 tACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATg                                                                                                                   >  2:250697/1‑148 (MQ=255)
                                  aCTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTCCATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATgg                                                                                                                  >  2:2666118/1‑148 (MQ=255)
                                  aCTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATgg                                                                                                                  <  1:1707771/148‑1 (MQ=255)
                                   cTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGt                                                                                                                 <  2:1727857/148‑1 (MQ=255)
                                   cTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGt                                                                                                                 <  2:1271163/148‑1 (MQ=255)
                                      atgaCGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTaa                                                                                                              <  2:2931856/148‑1 (MQ=255)
                                       tgaCGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTaaa                                                                                                             >  2:479570/1‑148 (MQ=255)
                                        gaCGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTaaaa                                                                                                            >  2:1587887/1‑148 (MQ=255)
                                          cGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAata                                                                                                          <  1:2769986/148‑1 (MQ=255)
                                          cGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTGTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAata                                                                                                          <  2:2024402/148‑1 (MQ=255)
                                           ggCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAatat                                                                                                         >  2:1335110/1‑148 (MQ=255)
                                             cGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTc                                                                                                       <  1:2601109/148‑1 (MQ=255)
                                               aTTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCAc                                                                                                     >  1:1326665/1‑148 (MQ=255)
                                               aTTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCAc                                                                                                     <  2:885131/148‑1 (MQ=255)
                                                tttGCCGTTTTCTGACGCCAATATCCAGCAACCCGGCAGGCAACCCACTTTCCGGGTCATTTGCGCCATGGACATAGATCACCAATTCTACAACCTGTAGCGTGGTAAATAAGGCGTGATAACAGCGGGGATGGTTAAAATATTCAcc                                                                                                    >  1:2910588/1‑148 (MQ=255)
                                                tttGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCAcc                                                                                                    >  2:2746259/1‑148 (MQ=255)
                                                 ttGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCAccc                                                                                                   >  1:397231/1‑148 (MQ=255)
                                                        ttGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTAt                                                                                            <  1:1569466/148‑1 (MQ=255)
                                                            tGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAat                                                                                        <  1:2159327/148‑1 (MQ=255)
                                                            tGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAat                                                                                        <  2:666263/148‑1 (MQ=255)
                                                              aCGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAatat                                                                                      <  2:2813715/148‑1 (MQ=255)
                                                                gCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCg                                                                                    <  1:122500/148‑1 (MQ=255)
                                                                 ccAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGc                                                                                   >  2:2922666/1‑148 (MQ=255)
                                                                  cAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGcc                                                                                  <  1:2721046/148‑1 (MQ=255)
                                                                  cAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGcc                                                                                  <  1:1068316/148‑1 (MQ=255)
                                                                   aaTATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGccc                                                                                 <  1:1235654/148‑1 (MQ=255)
                                                                   aaTATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGccc                                                                                 <  2:2831810/148‑1 (MQ=255)
                                                                   aaTATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGccc                                                                                 <  1:2842618/148‑1 (MQ=255)
                                                                      atCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTa                                                                              <  2:2286991/148‑1 (MQ=255)
                                                                      atCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTa                                                                              >  1:2207489/1‑148 (MQ=255)
                                                                        ccAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATc                                                                            >  1:817860/1‑148 (MQ=255)
                                                                         cAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCa                                                                           <  2:3306526/148‑1 (MQ=255)
                                                                             aaCCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTaa                                                                       <  1:657130/148‑1 (MQ=255)
                                                                                ccGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAAtt                                                                    <  1:1059853/148‑1 (MQ=255)
                                                                                 cGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTc                                                                   <  2:2950087/148‑1 (MQ=255)
                                                                                 cGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTc                                                                   <  1:1816205/148‑1 (MQ=255)
                                                                                 cGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTc                                                                   <  1:1816044/148‑1 (MQ=255)
                                                                                 cGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTc                                                                   <  1:2733754/148‑1 (MQ=255)
                                                                                  tgcaggcaACAGACTTTACGGGTCATTTGAGCCATGGACATAGATCAGCATGTATACATCCTGTAGCGTGTTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCAACCGGCTTATCAATATCTCCCTTATCATTAAATTcc                                                                  <  1:250624/147‑1 (MQ=255)
                                                                                  ggcaggcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTcc                                                                  <  1:250767/148‑1 (MQ=255)
                                                                                   gcaggcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATGAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCa                                                                 >  1:2630398/1‑148 (MQ=255)
                                                                                   gcaggcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCa                                                                 <  1:1949974/148‑1 (MQ=255)
                                                                                      ggcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACgg                                                              >  2:1067678/1‑148 (MQ=255)
                                                                                       gcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACggt                                                             <  2:2577549/148‑2 (MQ=255)
                                                                                       gcaACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCAACggt                                                             >  2:1691656/1‑147 (MQ=255)
                                                                                          aCCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtg                                                          >  2:2462383/1‑148 (MQ=255)
                                                                                           ccGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCtgt                                                         >  2:566891/1‑148 (MQ=255)
                                                                                             gACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTgtgt                                                       <  2:2188442/148‑1 (MQ=255)
                                                                                              aCTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGtt                                                      >  1:2730535/1‑148 (MQ=255)
                                                                                                 ttCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGttttt                                                   >  1:2426011/1‑148 (MQ=255)
                                                                                                   ccGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCt                                                 <  2:2413821/148‑1 (MQ=255)
                                                                                                   ccGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCt                                                 <  1:2741314/148‑1 (MQ=255)
                                                                                                      ggTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGc                                              >  1:16015/1‑148 (MQ=255)
                                                                                                      ggTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGc                                              >  1:1692422/1‑148 (MQ=255)
                                                                                                           tttGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAggg                                         <  2:922948/148‑1 (MQ=255)
                                                                                                           tttGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAggg                                         <  2:3061103/148‑1 (MQ=255)
                                                                                                            ttGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGt                                        <  1:1037550/148‑1 (MQ=255)
                                                                                                                   aTGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATCATTc                                 >  2:1738697/1‑148 (MQ=255)
                                                                                                                     ggACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTccc                               <  1:3242635/148‑1 (MQ=255)
                                                                                                                     ggACATAGATCAGCATGTCTACATCCTGTTGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTccc                               <  1:566891/148‑1 (MQ=255)
                                                                                                                     ggACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTccc                               <  2:1635959/148‑1 (MQ=255)
                                                                                                                           aGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCaa                         >  1:786016/1‑148 (MQ=255)
                                                                                                                            gATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAAc                        >  2:1663971/1‑148 (MQ=255)
                                                                                                                              tCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGc                      <  2:2910588/148‑1 (MQ=255)
                                                                                                                               cAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCa                     >  2:3222088/1‑148 (MQ=255)
                                                                                                                                aGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCAt                    >  2:107938/1‑148 (MQ=255)
                                                                                                                                  cATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTa                  >  2:2787850/1‑148 (MQ=255)
                                                                                                                                     gTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAAt               <  1:2496114/148‑1 (MQ=255)
                                                                                                                                     gTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAAt               >  1:2152232/1‑148 (MQ=255)
                                                                                                                                     gTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAAt               <  2:2252988/148‑1 (MQ=255)
                                                                                                                                       cTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATaa             <  2:1629574/148‑1 (MQ=255)
                                                                                                                                          cATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgt          >  2:1645998/1‑148 (MQ=255)
                                                                                                                                            tCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgttg        <  1:498552/148‑1 (MQ=255)
                                                                                                                                            tCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgttg        >  1:2356669/1‑148 (MQ=255)
                                                                                                                                             ccTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACgttgt       <  2:3117993/148‑1 (MQ=255)
                                                                                                                                               tGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGttt     <  2:1332529/148‑1 (MQ=255)
                                                                                                                                                gTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGtttt    >  2:2547496/1‑148 (MQ=255)
                                                                                                                                                gTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGTCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGtttt    >  2:891841/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGCGGCGACGTCAGCATCTGGCATGTCCGTTTTACTGATGACGGCGATTTGCCGTTTGCTGACGCCAATATCCAGCAACCCGGCAGGCAACCGACTTTCCGGGTCATTTGCGCCATGGACATAGATCAGCATGTCTACATCCTGTAGCGTGGTAATTAAGGCGTGATACCAGCGGGGATGGTTAAAATATTCACCCGGCGTATCAATATCGCCCTTATCATTAAATTCCACGGCCTGTGTTTTTCTGGCGAGGGTATAATTCCCCTGCAACGCATTAAATAACGTTGTTTTTC  >  CP009273/2568486‑2568779

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: