Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,893,446 T→A 93.1% K63* (AAA→TAA)  ygcW ← putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,893,4460TA93.1% 869.8 / 41.2 282K63* (AAA→TAA) ygcWputative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (9/11);  new base A (136/126);  total (145/137)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 282 total aligned reads displayed.

CACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGA  >  CP009273/2893303‑2893590
                                                                                                                                               |                                                                                                                                                
 aCCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCaa                                                                                                                                             >  1:1167703/1‑148 (MQ=255)
    ttATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCa                                                                                                                                          <  1:1488481/148‑1 (MQ=255)
       ttCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAttt                                                                                                                                       >  1:2190307/1‑148 (MQ=255)
          aGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCt                                                                                                                                    <  1:2901401/148‑1 (MQ=255)
           gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAAtcatttcatt                                                                                                                                   >  2:2138854/1‑148 (MQ=255)
           gCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCtt                                                                                                                                   >  1:2808162/1‑148 (MQ=255)
            cTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCttt                                                                                                                                  >  2:1465240/1‑148 (MQ=255)
             ttACAAAGACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTg                                                                                                                                 <  2:1339390/148‑1 (MQ=255)
              tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGt                                                                                                                                >  1:712619/1‑148 (MQ=255)
              tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGt                                                                                                                                >  1:235341/1‑148 (MQ=255)
               aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGtt                                                                                                                               <  1:1001852/148‑1 (MQ=255)
                cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGttt                                                                                                                              >  1:1330262/1‑148 (MQ=255)
                cAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGttt                                                                                                                              <  1:2024044/148‑1 (MQ=255)
                 aaaTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                             >  1:2580905/1‑148 (MQ=255)
                  aaTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                            >  2:2717283/1‑148 (MQ=255)
                  aaTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCg                                                                                                                            >  1:3321375/1‑148 (MQ=255)
                   aTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGc                                                                                                                           <  1:1439633/148‑1 (MQ=255)
                   aTACCGGCATTGTTAACCAGAATATCAACTGTAACTAAACGCTCACAGCAAGCAGCGATAAACTTCTGCGGCGCGCCTTCTGCGGTGATAACCACATGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGc                                                                                                                           >  2:2100631/1‑148 (MQ=255)
                    tACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGcc                                                                                                                          >  1:1633379/1‑148 (MQ=255)
                    tACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGcc                                                                                                                          <  1:2709869/148‑1 (MQ=255)
                     aCCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCg                                                                                                                         <  2:1944739/148‑1 (MQ=255)
                       cGGCATGGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGtt                                                                                                                       <  2:614548/148‑1 (MQ=255)
                          cATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATc                                                                                                                    <  2:1651040/148‑1 (MQ=255)
                          cATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATc                                                                                                                    >  1:2271035/1‑148 (MQ=255)
                           aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCt                                                                                                                   <  1:331878/148‑1 (MQ=255)
                             tGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                 >  2:1030987/1‑148 (MQ=255)
                             tGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                 >  2:2615262/1‑148 (MQ=255)
                              gTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCTCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTg                                                                                                                <  1:1467525/148‑1 (MQ=255)
                               ttAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGa                                                                                                               >  2:200150/1‑148 (MQ=255)
                                tAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGAc                                                                                                              >  2:1128089/1‑148 (MQ=255)
                                  aCCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGa                                                                                                            >  1:3304736/1‑148 (MQ=255)
                                  aCCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGa                                                                                                            >  2:2263562/1‑148 (MQ=255)
                                    cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaca                                                                                                          <  2:554578/148‑3 (MQ=255)
                                    cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaaa                                                                                                          <  1:6250/148‑1 (MQ=255)
                                    cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaaa                                                                                                          <  2:781322/148‑1 (MQ=255)
                                    cAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGaaa                                                                                                          <  1:2944177/148‑1 (MQ=255)
                                     aGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAAc                                                                                                         >  2:275428/1‑148 (MQ=255)
                                      gAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACt                                                                                                        >  1:1597015/1‑148 (MQ=255)
                                      gAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACt                                                                                                        >  2:2326485/1‑148 (MQ=255)
                                       aaTATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCATTATCTTTGACGAAACTa                                                                                                       >  2:3242749/1‑148 (MQ=255)
                                        atatCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAg                                                                                                      <  2:1248361/148‑1 (MQ=255)
                                         tatCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAgg                                                                                                     <  2:161812/148‑1 (MQ=255)
                                          atCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGa                                                                                                    >  1:1817849/1‑148 (MQ=255)
                                          atCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGa                                                                                                    <  2:1424361/148‑1 (MQ=255)
                                          atCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGa                                                                                                    >  1:396472/1‑148 (MQ=255)
                                              aCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaa                                                                                                >  1:1291464/1‑148 (MQ=255)
                                                tGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAg                                                                                              <  2:1759133/148‑1 (MQ=255)
                                                 gTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGa                                                                                             <  1:3038939/148‑1 (MQ=255)
                                                  tACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGat                                                                                            <  2:1516022/148‑1 (MQ=255)
                                                  tACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGat                                                                                            <  2:2650405/148‑1 (MQ=255)
                                                     cGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATAtt                                                                                         >  1:856699/1‑148 (MQ=255)
                                                       aaaCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTg                                                                                       <  2:3145559/148‑1 (MQ=255)
                                                        aaCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgc                                                                                      >  2:1168382/1‑148 (MQ=255)
                                                         aCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgcg                                                                                     >  1:2512620/1‑148 (MQ=255)
                                                          cGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgcgc                                                                                    >  2:3219968/1‑148 (MQ=255)
                                                           gCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGcc                                                                                   <  1:3372609/148‑1 (MQ=255)
                                                            cTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCa                                                                                  >  1:2255415/1‑148 (MQ=255)
                                                             tCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAg                                                                                 >  2:2498548/1‑148 (MQ=255)
                                                              cacaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGc                                                                                >  2:733589/1‑148 (MQ=255)
                                                              cacaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGc                                                                                <  2:535768/148‑1 (MQ=255)
                                                                caGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCtt                                                                              <  2:3058013/148‑1 (MQ=255)
                                                                 agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                             <  1:1308171/148‑1 (MQ=255)
                                                                 agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                             <  2:486101/148‑1 (MQ=255)
                                                                     agcagcGATAATCTTCTGCGGCGCGCCTTCAGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGcc                                                                         >  2:3468217/1‑148 (MQ=255)
                                                                      gcagcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCa                                                                        >  1:474400/1‑148 (MQ=255)
                                                                         gcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGACGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACg                                                                     >  1:599573/1‑148 (MQ=255)
                                                                              aaTCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATg                                                                >  1:349923/1‑148 (MQ=255)
                                                                              aaTCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATg                                                                >  1:1732734/1‑148 (MQ=255)
                                                                                 cttctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCa                                                             >  1:3086723/1‑148 (MQ=255)
                                                                                      gcggcgCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGc                                                        >  2:1838328/1‑148 (MQ=255)
                                                                                      gcggcgCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGc                                                        <  1:1696045/148‑1 (MQ=255)
                                                                                       cggcgCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGcc                                                       <  1:1238358/148‑1 (MQ=255)
                                                                                           gcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGc                                                   <  2:390949/148‑1 (MQ=255)
                                                                                           gcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGc                                                   <  1:2538575/148‑1 (MQ=255)
                                                                                                ttCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTaaa                                              <  1:1856425/148‑1 (MQ=255)
                                                                                                   tGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCg                                           <  2:3337942/148‑1 (MQ=255)
                                                                                                        tGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAtt                                      <  1:2001948/148‑1 (MQ=255)
                                                                                                        tGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTAtt                                      <  2:2799224/148‑1 (MQ=255)
                                                                                                         gATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTCAACCGCTATTc                                     >  2:3019026/1‑148 (MQ=255)
                                                                                                           tACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTccc                                   >  2:2531551/1‑148 (MQ=255)
                                                                                                           tACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTccc                                   <  1:890101/148‑1 (MQ=255)
                                                                                                              ccACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                <  2:2454744/148‑1 (MQ=255)
                                                                                                              ccACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                >  1:429673/1‑148 (MQ=255)
                                                                                                              ccACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCAcc                                >  1:2993346/1‑148 (MQ=255)
                                                                                                                aCCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCgg                              <  2:2067116/148‑1 (MQ=255)
                                                                                                                    gCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAAAAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGAAAATGCCTGGCCTAAACCGCTATTCCCACCGGTaaa                          >  1:2795484/1‑147 (MQ=255)
                                                                                                                     cATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACa                         <  1:1168382/148‑1 (MQ=255)
                                                                                                                     cATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACa                         >  2:1295247/1‑148 (MQ=255)
                                                                                                                       tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAAt                       <  2:301287/148‑1 (MQ=255)
                                                                                                                       tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAAt                       <  2:2580346/148‑1 (MQ=255)
                                                                                                                        gAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAAtt                      <  1:1651451/148‑1 (MQ=255)
                                                                                                                             ccACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGt                 >  2:2413945/1‑148 (MQ=255)
                                                                                                                               aCCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt               <  2:3203719/148‑1 (MQ=255)
                                                                                                                                ccTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGtttt              <  1:1922538/148‑1 (MQ=255)
                                                                                                                                 cTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTa             >  1:216794/1‑148 (MQ=255)
                                                                                                                                    aaCACCCTGTTTTTCAATCATTTCCTTGGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCt          <  2:903898/148‑1 (MQ=255)
                                                                                                                                      cacCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACAttt        <  1:2138465/148‑1 (MQ=255)
                                                                                                                                         ccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAg     >  1:1375260/1‑148 (MQ=255)
                                                                                                                                           tGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg   <  2:1602101/148‑1 (MQ=255)
                                                                                                                                           tGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg   <  1:412852/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGA  >  CP009273/2893303‑2893590

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: