Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,545,608 C→A 100% R177I (AGA→ATA)  yjiC ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,545,6080CA99.2% 806.5 / ‑2.0 246R177I (AGA→ATA) yjiCuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (108/136);  total (110/136)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 246 total aligned reads displayed.

TCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTA  >  CP009273/4545464‑4545752
                                                                                                                                                |                                                                                                                                                
tCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATAttt                                                                                                                                               <  1:1162550/148‑1 (MQ=255)
   ataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTg                                                                                                                                            <  1:534848/148‑1 (MQ=255)
   ataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTg                                                                                                                                            <  2:1743507/148‑1 (MQ=255)
   ataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTg                                                                                                                                            <  1:2000530/148‑1 (MQ=255)
    taCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGc                                                                                                                                           >  1:2543881/1‑148 (MQ=255)
     aCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGcc                                                                                                                                          <  2:3464135/148‑1 (MQ=255)
      ccAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCt                                                                                                                                         >  2:1613041/1‑148 (MQ=255)
         attattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTgag                                                                                                                                      <  2:2864736/148‑1 (MQ=255)
          ttattGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTgaga                                                                                                                                     <  2:3486342/148‑1 (MQ=255)
            attGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAg                                                                                                                                   <  1:1247098/148‑1 (MQ=255)
             ttGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGc                                                                                                                                  <  1:297842/148‑1 (MQ=255)
             ttGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGc                                                                                                                                  <  2:1767808/148‑1 (MQ=255)
             ttGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGc                                                                                                                                  <  1:2154400/148‑1 (MQ=255)
              tGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCt                                                                                                                                 >  2:190108/1‑148 (MQ=255)
               gggTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTg                                                                                                                                <  2:1117413/148‑1 (MQ=255)
                 gTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGa                                                                                                                              >  1:3360270/1‑148 (MQ=255)
                  tCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGaa                                                                                                                             >  2:1454359/1‑148 (MQ=255)
                  tCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGaa                                                                                                                             >  2:1454258/1‑148 (MQ=255)
                    aTTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAgg                                                                                                                           <  1:2528677/148‑1 (MQ=255)
                     ttACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAggg                                                                                                                          >  2:2930707/1‑148 (MQ=255)
                       aCCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGtt                                                                                                                        <  1:3404065/148‑1 (MQ=255)
                            aCTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAat                                                                                                                   >  2:3095509/1‑148 (MQ=255)
                             cTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAata                                                                                                                  >  1:1678521/1‑148 (MQ=255)
                              tCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAatat                                                                                                                 <  2:127982/148‑1 (MQ=255)
                                 atatTAGCATTAAGTAATGCGGCAAGGGCGACTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              >  1:1289939/1‑148 (MQ=255)
                                 atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              <  1:2787956/148‑1 (MQ=255)
                                 atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              >  2:1067557/1‑148 (MQ=255)
                                 atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                                              <  1:3056214/148‑1 (MQ=255)
                                  tatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATcc                                                                                                             <  1:512995/148‑1 (MQ=255)
                                   atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            >  2:1837572/1‑148 (MQ=255)
                                   atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGCCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            <  2:3255566/148‑1 (MQ=255)
                                   atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATAAATATAGATATCAATCCTTGGCTGCTCATTGGTCAATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            <  2:948524/148‑1 (MQ=255)
                                   atTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGACCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                                                                                            <  1:2407674/148‑1 (MQ=255)
                                    ttAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTa                                                                                                           <  2:2633824/148‑1 (MQ=255)
                                     tAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAt                                                                                                          >  1:2777001/1‑148 (MQ=255)
                                     tAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAt                                                                                                          <  2:1043890/148‑1 (MQ=255)
                                     tAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAt                                                                                                          <  2:2746056/148‑1 (MQ=255)
                                     tAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAt                                                                                                          <  2:2345041/148‑1 (MQ=255)
                                     tAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAt                                                                                                          >  1:2092050/1‑148 (MQ=255)
                                          ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTgg                                                                                                     <  1:1608329/148‑1 (MQ=255)
                                             aGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTc                                                                                                  <  2:169939/148‑1 (MQ=255)
                                                 aTGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGc                                                                                              >  1:2845812/1‑148 (MQ=255)
                                                  tGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                                             >  2:1338882/1‑148 (MQ=255)
                                                  tGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                                             >  2:2447677/1‑148 (MQ=255)
                                                  tGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                                             <  1:1855172/148‑1 (MQ=255)
                                                    cGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATa                                                                                           <  1:459731/148‑1 (MQ=255)
                                                          gggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                                                                     >  2:2545779/1‑148 (MQ=255)
                                                           ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa                                                                                    >  1:3145769/1‑148 (MQ=255)
                                                             cGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAc                                                                                  >  2:1222145/1‑148 (MQ=255)
                                                               aaTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAt                                                                                >  1:238505/1‑148 (MQ=255)
                                                                aTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAtt                                                                               >  1:1021204/1‑148 (MQ=255)
                                                                aTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAtt                                                                               >  2:487178/1‑148 (MQ=255)
                                                                 tCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                                              <  1:294203/148‑1 (MQ=255)
                                                                     gaagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAc                                                                          >  2:3146087/1‑148 (MQ=255)
                                                                        gaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAg                                                                       >  1:667157/1‑148 (MQ=255)
                                                                           tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGc                                                                    <  2:235453/148‑1 (MQ=255)
                                                                            ttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGcc                                                                   <  2:648970/148‑1 (MQ=255)
                                                                               cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTaa                                                                <  1:732562/148‑1 (MQ=255)
                                                                                  ccccTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGt                                                             >  2:2536402/1‑148 (MQ=255)
                                                                                     cTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                                          <  1:3062324/148‑1 (MQ=255)
                                                                                     cTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                                          <  2:1987301/148‑1 (MQ=255)
                                                                                      tGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGACATAACCAGGAACATTGCAACCAGGGCCTAATGTTAtt                                                         <  2:1645440/148‑1 (MQ=255)
                                                                                       gTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAttt                                                        <  1:3194791/148‑1 (MQ=255)
                                                                                         tACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCa                                                      >  2:517092/1‑148 (MQ=255)
                                                                                         tACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCa                                                      <  1:3459570/148‑1 (MQ=255)
                                                                                            cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGa                                                   >  1:2559049/1‑148 (MQ=255)
                                                                                              tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAAc                                                 <  1:518656/148‑1 (MQ=255)
                                                                                               cATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACt                                                >  2:3309238/1‑148 (MQ=255)
                                                                                                aTTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTg                                               <  2:3360277/148‑1 (MQ=255)
                                                                                                      tGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGt                                         <  1:1338882/148‑1 (MQ=255)
                                                                                                       gATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTc                                        <  1:3336868/148‑1 (MQ=255)
                                                                                                           aaTATAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCa                                    <  2:504575/148‑1 (MQ=255)
                                                                                                            atatAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa                                   >  1:889524/1‑148 (MQ=255)
                                                                                                            atatAGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa                                   <  1:1242121/148‑1 (MQ=255)
                                                                                                             tataGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAAc                                  >  2:2243565/1‑148 (MQ=255)
                                                                                                              ataGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACa                                 <  1:1313862/148‑1 (MQ=255)
                                                                                                                aGATATCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAgg                               >  1:2920335/1‑148 (MQ=255)
                                                                                                                  atatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGAt                             >  2:2665841/1‑148 (MQ=255)
                                                                                                                  atatCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGAt                             >  1:2658423/1‑148 (MQ=255)
                                                                                                                     tCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGa                          >  2:2793912/1‑148 (MQ=255)
                                                                                                                     tCAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGa                          >  2:2422705/1‑148 (MQ=255)
                                                                                                                      cAATCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGaa                         >  2:585538/1‑148 (MQ=255)
                                                                                                                        aTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAg                       <  2:1369190/148‑1 (MQ=255)
                                                                                                                        aTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAg                       >  2:26471/1‑148 (MQ=255)
                                                                                                                        aTCCTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAg                       >  1:440140/1‑148 (MQ=255)
                                                                                                                          ccTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAggg                     <  2:2845801/148‑1 (MQ=255)
                                                                                                                          ccTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAggg                     <  2:1333396/148‑1 (MQ=255)
                                                                                                                          ccTTGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAggg                     <  2:2928815/148‑1 (MQ=255)
                                                                                                                            ttGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGct                   >  2:3431863/1‑148 (MQ=255)
                                                                                                                            ttGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGct                   <  2:520964/148‑1 (MQ=255)
                                                                                                                             tGGCTGCTCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGctc                  <  1:3438681/148‑1 (MQ=255)
                                                                                                                               gctgctCATTGGTCCATATTTGTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCtt                <  2:409925/148‑1 (MQ=255)
                                                                                                                                ctgctCATTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTc               <  1:4262/148‑1 (MQ=255)
                                                                                                                                ctgctCATTGGTCCATATATTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTc               >  2:883209/1‑148 (MQ=255)
                                                                                                                                      aTTGGTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCataaat         <  1:2930650/148‑1 (MQ=255)
                                                                                                                                         ggTCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCAtaaataaa      <  1:1392632/148‑1 (MQ=255)
                                                                                                                                           tCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAtt    >  1:107294/1‑148 (MQ=255)
                                                                                                                                           tCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAtt    >  1:1946046/1‑148 (MQ=255)
                                                                                                                                           tCCATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAAtt    <  1:2836565/148‑1 (MQ=255)
                                                                                                                                             cATATTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTa  <  1:2714547/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATCGAAAGGGCTCTTCATAAATAAATTTA  >  CP009273/4545464‑4545752

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: