Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 332,897 G→T 92.8% D137Y (GAT→TAT)  yahF → putative NAD(P)‑binding succinyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273332,8970GT92.8% 618.1 / 36.1 223D137Y (GAT→TAT) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base T (101/106);  total (110/113)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 223 total aligned reads displayed.

AACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTGGTAAC  >  CP009273/332753‑333043
                                                                                                                                                |                                                                                                                                                  
aaCACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTAtc                                                                                                                                                 >  2:376691/1‑148 (MQ=255)
aaCACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTAtc                                                                                                                                                 <  1:288563/148‑1 (MQ=255)
aaCACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTAtc                                                                                                                                                 >  2:2467335/1‑148 (MQ=255)
aaCACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTAtc                                                                                                                                                 <  2:1655649/148‑1 (MQ=255)
 acacCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGTCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTAtct                                                                                                                                                >  1:2424589/1‑148 (MQ=255)
 acacCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTAtct                                                                                                                                                >  1:3069380/1‑148 (MQ=255)
    ccAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCaa                                                                                                                                             <  2:2554646/148‑1 (MQ=255)
      aaagcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCCGGCGCTGCAAAACTATCTCAACg                                                                                                                                           >  2:1970901/1‑148 (MQ=255)
      aaagcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACg                                                                                                                                           <  1:2507418/148‑1 (MQ=255)
        agcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTg                                                                                                                                         >  2:906140/1‑148 (MQ=255)
        agcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTg                                                                                                                                         <  1:3344147/148‑1 (MQ=255)
         gcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGa                                                                                                                                        >  2:1378680/1‑148 (MQ=255)
         gcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGa                                                                                                                                        <  2:3159504/148‑1 (MQ=255)
          cgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGAt                                                                                                                                       >  2:3190610/1‑148 (MQ=255)
           gcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATg                                                                                                                                      <  1:2504190/148‑1 (MQ=255)
             aaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCt                                                                                                                                    >  2:430325/1‑148 (MQ=255)
             aaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCt                                                                                                                                    <  1:2060544/148‑1 (MQ=255)
              aagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTg                                                                                                                                   <  1:1115945/148‑1 (MQ=255)
                gcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGtt                                                                                                                                 <  1:156586/148‑1 (MQ=255)
                       gCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGAt                                                                                                                          <  1:305641/148‑1 (MQ=255)
                        cACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGTCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATa                                                                                                                         <  2:2424589/148‑1 (MQ=255)
                         aCGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATaa                                                                                                                        <  2:635970/148‑1 (MQ=255)
                              gcgcGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                                                                                   <  2:723118/148‑1 (MQ=255)
                              gcgcGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                                                                                   <  1:151547/148‑1 (MQ=255)
                              gcgcGTTACGCCACTAGTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGtct                                                                                                                   <  2:723191/148‑1 (MQ=255)
                                gcgTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCa                                                                                                                 <  1:591960/148‑1 (MQ=255)
                                 cgTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAg                                                                                                                >  1:272002/1‑148 (MQ=255)
                                 cgTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAg                                                                                                                >  1:272143/1‑148 (MQ=255)
                                 cgTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAg                                                                                                                >  1:953009/1‑148 (MQ=255)
                                   ttACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGtt                                                                                                              >  1:925563/1‑148 (MQ=255)
                                   ttACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGtt                                                                                                              >  2:1522761/1‑148 (MQ=255)
                                    tACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCACAACGATCTCAACGTGATGCTGTTTTACGATAACGTCTCAGTTg                                                                                                             >  2:142620/1‑148 (MQ=255)
                                     aCGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGa                                                                                                            <  2:1450058/148‑1 (MQ=255)
                                        ccACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGa                                                                                                         >  1:2068999/1‑148 (MQ=255)
                                         cACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGAt                                                                                                        >  1:2392419/1‑148 (MQ=255)
                                         cACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGAt                                                                                                        >  1:3392438/1‑148 (MQ=255)
                                            tATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGaa                                                                                                     <  1:3367091/148‑1 (MQ=255)
                                             aTTGGCAGCGCCAAAAAGCCTATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTGTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAAc                                                                                                    <  1:1524015/148‑1 (MQ=255)
                                             aTTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAAc                                                                                                    <  1:1524323/148‑1 (MQ=255)
                                              ttGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACt                                                                                                   <  2:2455738/148‑1 (MQ=255)
                                                 gCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGgc                                                                                                >  1:1779654/1‑148 (MQ=255)
                                                  cAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGgcg                                                                                               <  1:1390763/148‑1 (MQ=255)
                                                    gcgcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCt                                                                                             <  2:3431381/148‑1 (MQ=255)
                                                    gcgcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCt                                                                                             >  2:1859538/1‑148 (MQ=255)
                                                    gcgcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCt                                                                                             >  2:1859754/1‑148 (MQ=255)
                                                      gcCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCa                                                                                           <  2:2652970/148‑1 (MQ=255)
                                                       ccAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCaa                                                                                          >  2:609186/1‑148 (MQ=255)
                                                        cAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAg                                                                                         >  1:326275/1‑148 (MQ=255)
                                                         aaaaaGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGc                                                                                        >  1:1335783/1‑148 (MQ=255)
                                                         aaaaaGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGc                                                                                        <  1:3505608/148‑1 (MQ=255)
                                                         aaaaaGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGc                                                                                        <  2:2683862/148‑1 (MQ=255)
                                                          aaaaGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCa                                                                                       >  1:1854746/1‑148 (MQ=255)
                                                           aaaGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCaa                                                                                      <  2:2393116/148‑1 (MQ=255)
                                                               cATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTgg                                                                                  >  1:2905324/1‑148 (MQ=255)
                                                                 tatCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGcc                                                                                >  1:446734/1‑148 (MQ=255)
                                                                 tatCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGcc                                                                                >  2:1477069/1‑148 (MQ=255)
                                                                        gAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGaaa                                                                         <  1:3015406/148‑1 (MQ=255)
                                                                         aaaGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGaaaa                                                                        <  2:894887/148‑1 (MQ=255)
                                                                           aGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAg                                                                      <  2:3264792/148‑1 (MQ=255)
                                                                           aGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAg                                                                      <  2:206168/148‑1 (MQ=255)
                                                                              aaCCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGgc                                                                   <  2:44909/148‑1 (MQ=255)
                                                                                ccTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGgctg                                                                 <  2:790853/148‑1 (MQ=255)
                                                                                ccTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGgctg                                                                 <  1:1378680/148‑1 (MQ=255)
                                                                                   ggCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGctgctg                                                              <  1:3418060/148‑1 (MQ=255)
                                                                                     cGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCtgat                                                            >  2:2294785/1‑148 (MQ=255)
                                                                                      ggTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCtgatg                                                           <  1:2982176/148‑1 (MQ=255)
                                                                                        tGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTgatgat                                                         >  1:2203187/1‑148 (MQ=255)
                                                                                         gATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGatgatg                                                        >  1:469779/1‑148 (MQ=255)
                                                                                          aTTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATgg                                                       >  1:2760210/1‑148 (MQ=255)
                                                                                          aTTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATgg                                                       <  2:3079896/148‑1 (MQ=255)
                                                                                            ttCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATgggg                                                     <  2:2336000/148‑1 (MQ=255)
                                                                                            ttCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATgggg                                                     >  2:2783395/1‑148 (MQ=255)
                                                                                                 tCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGa                                                >  1:368621/1‑148 (MQ=255)
                                                                                                  cAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGAc                                               >  2:3034405/1‑148 (MQ=255)
                                                                                                  cAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGAc                                               >  1:2915771/1‑148 (MQ=255)
                                                                                                    aCGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACtg                                             >  2:1294098/1‑148 (MQ=255)
                                                                                                     cGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACtgt                                            <  1:1455396/148‑1 (MQ=255)
                                                                                                     cGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACtgt                                            >  1:416950/1‑148 (MQ=255)
                                                                                                         cTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCa                                        >  1:656306/1‑148 (MQ=255)
                                                                                                            tttGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACgg                                     <  2:363809/148‑1 (MQ=255)
                                                                                                              tGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCg                                   <  1:2153921/148‑1 (MQ=255)
                                                                                                               gCCGCTCGCGACGCGCGTCAGGCGCTGCAAAACTATCTCACCGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGa                                  <  2:2158091/148‑1 (MQ=255)
                                                                                                                ccGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGAt                                 <  2:765108/148‑1 (MQ=255)
                                                                                                                 cGCTCGCGAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGAtt                                <  2:1746385/148‑1 (MQ=255)
                                                                                                                     cgcgAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCa                            <  1:1818022/148‑1 (MQ=255)
                                                                                                                       cgAAGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAAc                          >  1:1320354/1‑148 (MQ=255)
                                                                                                                         aaGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAAcgg                        >  2:2590634/1‑148 (MQ=255)
                                                                                                                         aaGCGCGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAAcgg                        >  2:796372/1‑148 (MQ=255)
                                                                                                                           gcgcGTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAAcggcg                      <  2:1564207/148‑1 (MQ=255)
                                                                                                                            cgcgTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAAcggcgc                     >  2:3363542/1‑148 (MQ=255)
                                                                                                                              cgTCAGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAAcggcgcgg                   <  1:3190610/148‑1 (MQ=255)
                                                                                                                                  aGGCGCTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGct               >  1:1637066/1‑148 (MQ=255)
                                                                                                                                      gcTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGt           <  2:3275191/148‑1 (MQ=255)
                                                                                                                                      gcTGCAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGt           <  1:1693708/148‑1 (MQ=255)
                                                                                                                                          cAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTg       >  2:1619022/1‑148 (MQ=255)
                                                                                                                                          cAAAACTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTg       >  2:1080594/1‑148 (MQ=255)
                                                                                                                                           aaaaCTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTgg      >  1:1279738/1‑148 (MQ=255)
                                                                                                                                             aaCTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTGGTa    <  1:1974221/148‑1 (MQ=255)
                                                                                                                                             aaCTATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTGGTa    <  1:2304028/148‑1 (MQ=255)
                                                                                                                                               cGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTGGTAAc  >  1:1033840/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTCTCAGTTGAAGATGAACTGGCGCTCAAGCAACTGGCCCACGAAAAAGGGCTGCTGATGATGGGGCCAGACTGTGGCACGGCGATTATCAACGGCGCGGCGCTCTGTTTTGGTAAC  >  CP009273/332753‑333043

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: