Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,470,017 G→A 100% M157I (ATG→ATA yjgN → DUF898 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,470,0170GA99.4% 1054.6 / ‑2.6 319M157I (ATG→ATAyjgNDUF898 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (137/180);  total (139/180)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-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 322 total aligned reads displayed.

TTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTT  >  CP009273/4469872‑4470160
                                                                                                                                                 |                                                                                                                                               
 tttGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATaggg                                                                                                                                              >  1:1679126/1‑148 (MQ=255)
  ttGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGa                                                                                                                                             <  2:1216543/148‑1 (MQ=255)
   tGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGAc                                                                                                                                            >  1:1572062/1‑148 (MQ=255)
      ttATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGt                                                                                                                                         <  2:2104607/148‑1 (MQ=255)
       tataTTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGtt                                                                                                                                        <  2:1995552/148‑1 (MQ=255)
       tataTTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGtt                                                                                                                                        <  2:2909564/148‑1 (MQ=255)
        atatTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGttt                                                                                                                                       <  2:2608233/148‑1 (MQ=255)
          atTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTc                                                                                                                                     <  2:2641888/148‑1 (MQ=255)
          atTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTc                                                                                                                                     >  1:2691422/1‑148 (MQ=255)
            ttATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCtt                                                                                                                                   >  2:903004/1‑148 (MQ=255)
                ggCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTAtc                                                                                                                               >  1:1013941/1‑148 (MQ=255)
                ggCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTAtc                                                                                                                               <  1:1297264/148‑1 (MQ=255)
                ggCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTAtc                                                                                                                               <  1:1298412/148‑1 (MQ=255)
                   agcaAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCg                                                                                                                            <  1:3658059/148‑1 (MQ=255)
                     caAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGAc                                                                                                                          >  2:2538076/1‑148 (MQ=255)
                     caAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGAc                                                                                                                          >  2:2279084/1‑148 (MQ=255)
                       aaaGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACaa                                                                                                                        <  1:542036/148‑1 (MQ=255)
                         aGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAg                                                                                                                      <  1:3377569/148‑1 (MQ=255)
                         aGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAg                                                                                                                      <  1:1055714/148‑1 (MQ=255)
                          ggACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGa                                                                                                                     >  1:1327719/1‑148 (MQ=255)
                           gACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGAt                                                                                                                    >  2:1739185/1‑148 (MQ=255)
                           gACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGAt                                                                                                                    >  2:438317/1‑148 (MQ=255)
                            aCTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATg                                                                                                                   <  1:2036009/148‑1 (MQ=255)
                             cTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGc                                                                                                                  >  2:1867428/1‑148 (MQ=255)
                              tGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCt                                                                                                                 <  2:1378837/148‑1 (MQ=255)
                                   catcaGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTg                                                                                                            >  2:3207279/1‑148 (MQ=255)
                                      caGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCa                                                                                                         >  1:358520/1‑148 (MQ=255)
                                       aGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCaa                                                                                                        <  1:502889/148‑1 (MQ=255)
                                        ggCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                       >  2:2947731/1‑148 (MQ=255)
                                        ggCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                       <  1:601649/148‑1 (MQ=255)
                                        ggCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAAt                                                                                                       <  1:600571/148‑1 (MQ=255)
                                         gCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATa                                                                                                      >  1:1187663/1‑148 (MQ=255)
                                         gCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATa                                                                                                      >  2:539312/1‑148 (MQ=255)
                                                 gACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagt                                                                                              >  1:203229/1‑148 (MQ=255)
                                                  aCCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCGGGGGGGTGTTGCTACCTGCCAATAGTTAAagta                                                                                             >  2:3635999/1‑148 (MQ=255)
                                                   ccAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagtag                                                                                            <  1:26512/148‑1 (MQ=255)
                                                   ccAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAagtag                                                                                            <  2:23985/148‑1 (MQ=255)
                                                          cAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttatt                                                                                     <  2:962344/148‑1 (MQ=255)
                                                           aaCGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGttatta                                                                                    >  1:3464870/1‑148 (MQ=255)
                                                            aCGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTtattat                                                                                   >  2:1445382/1‑148 (MQ=255)
                                                             cGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATtata                                                                                  <  1:1315162/148‑1 (MQ=255)
                                                                cGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTTATGTCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATcc                                                                               <  2:208668/148‑1 (MQ=255)
                                                                 gTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCa                                                                              <  2:2739375/148‑1 (MQ=255)
                                                                  tAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCAt                                                                             <  1:3418495/148‑1 (MQ=255)
                                                                  tAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCAt                                                                             <  2:1133750/148‑1 (MQ=255)
                                                                   aaGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGTCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATg                                                                            >  2:2367117/1‑148 (MQ=255)
                                                                   aaGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATg                                                                            <  1:1764039/148‑1 (MQ=255)
                                                                    aGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATgg                                                                           <  1:3111825/148‑1 (MQ=255)
                                                                    aGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATgg                                                                           >  2:578946/1‑148 (MQ=255)
                                                                     gATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGt                                                                          >  2:3010373/1‑148 (MQ=255)
                                                                      aTTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTATTGTTATTATATCCATGGtt                                                                         <  2:1524398/148‑1 (MQ=255)
                                                                      aTTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGtt                                                                         >  1:969411/1‑148 (MQ=255)
                                                                       tttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTc                                                                        >  2:1085096/1‑148 (MQ=255)
                                                                       tttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTc                                                                        <  1:3071633/148‑1 (MQ=255)
                                                                        ttAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCt                                                                       >  1:1976179/1‑148 (MQ=255)
                                                                          aGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGa                                                                     <  2:313027/148‑1 (MQ=255)
                                                                           gTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGAt                                                                    <  1:3101261/148‑1 (MQ=255)
                                                                           gTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGAt                                                                    <  2:1609292/148‑1 (MQ=255)
                                                                             tttAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATgg                                                                  >  1:3545257/1‑148 (MQ=255)
                                                                              ttAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGc                                                                 >  1:320254/1‑148 (MQ=255)
                                                                                 aTTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAAt                                                              <  1:591845/148‑1 (MQ=255)
                                                                                  ttGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATa                                                             <  2:3343432/148‑1 (MQ=255)
                                                                                   tGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAg                                                            >  2:3090822/1‑148 (MQ=255)
                                                                                         aTGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTa                                                      <  2:521951/148‑1 (MQ=255)
                                                                                          tGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTAt                                                     >  1:2755636/1‑148 (MQ=255)
                                                                                           gAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAATATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGATGGTAtt                                                    >  1:1487742/1‑148 (MQ=255)
                                                                                            aaaGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTg                                                   <  1:1078692/148‑1 (MQ=255)
                                                                                             aaGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGt                                                  >  1:845682/1‑148 (MQ=255)
                                                                                               gggTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGttt                                                <  1:2991070/148‑1 (MQ=255)
                                                                                               gggTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGttt                                                <  2:3361366/148‑1 (MQ=255)
                                                                                                ggTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTc                                               <  1:1281268/148‑1 (MQ=255)
                                                                                                   tCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCAt                                            >  2:823098/1‑148 (MQ=255)
                                                                                                   tCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGACGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCAt                                            <  1:3635999/148‑1 (MQ=255)
                                                                                                       gtgggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGt                                        >  1:2806993/1‑148 (MQ=255)
                                                                                                         gggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAt                                      <  2:3136452/148‑1 (MQ=255)
                                                                                                         gggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAt                                      <  1:2589721/148‑1 (MQ=255)
                                                                                                          ggtgACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTAtt                                     >  2:3451578/1‑148 (MQ=255)
                                                                                                              aCCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTa                                 <  1:1906292/148‑1 (MQ=255)
                                                                                                              aCCTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTa                                 <  1:2521127/148‑1 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                <  2:1680437/148‑1 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                >  1:2059740/1‑148 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                <  1:822611/148‑1 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                <  2:1283073/148‑1 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTCGGTATTGTTTCCATTGGTATTTTTaa                                >  1:1427506/1‑148 (MQ=255)
                                                                                                               ccTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGGTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTaa                                <  1:2382372/148‑1 (MQ=255)
                                                                                                                cTTTTTCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAAt                               >  1:1270055/1‑148 (MQ=255)
                                                                                                                  ttttCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATgg                             <  1:2878069/148‑1 (MQ=255)
                                                                                                                   tttCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGt                            >  1:1119683/1‑148 (MQ=255)
                                                                                                                    ttCTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTa                           <  2:1327719/148‑1 (MQ=255)
                                                                                                                      cTTGCCGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACt                         <  1:748603/148‑1 (MQ=255)
                                                                                                                          ccGATTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTtat                     <  2:1808605/148‑1 (MQ=255)
                                                                                                                             aTTTTAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTtatata                  >  2:224280/1‑148 (MQ=255)
                                                                                                                              ttttAATGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAg                 >  1:2348203/1‑148 (MQ=255)
                                                                                                                                   aTGGCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTgg            <  2:3218347/148‑1 (MQ=255)
                                                                                                                                     ggCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTa          <  1:1947746/148‑1 (MQ=255)
                                                                                                                                      gCCATTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTaa         <  2:3551139/148‑1 (MQ=255)
                                                                                                                                         aTTGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGa      <  1:2228017/148‑1 (MQ=255)
                                                                                                                                           tGGGATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGt    >  2:927514/1‑148 (MQ=255)
                                                                                                                                            gggATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGtt   <  2:667757/148‑1 (MQ=255)
                                                                                                                                             ggATAGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGttt  >  2:3063454/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTTTGCTTATATTTATGGCAGCAAAAGGACTGCGTCATCAGGCCTTAATGACCAGTCTCAACGGCGTAAGATTTAGTTTTAATTGCTCTATGAAAGGGTTCTGGTGGGTGACCTTTTTCTTGCCGATTTTAATGGCCATTGGGATGGGGACTGTTTTCTTTATCTCGACAAAGATGCTACCTGCCAATAGTTCAAGTAGTGTTATTATATCCATGGTTCTGATGGCAATAGTTGGTATTGTTTCCATTGGTATTTTTAATGGTACTTTATATAGTCTGGTAATGAGTTT  >  CP009273/4469872‑4470160

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: