Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,800 A→G 7.4% S202P (TCG→CCG) ‡ ynjF ← CDP‑alcohol phosphatidyltransferase family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,840,8000AG 7.4% 194.5 / 5.6 68S202P (TCG→CCG) ‡ynjFCDP‑alcohol phosphatidyltransferase family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base G (2/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCC  >  NC_000913/1840652‑1840948
                                                                                                                                                    |                                                                                                                                                    
aTGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTAGGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGa                                                                                                                                                      <  2:781270/149‑1 (MQ=255)
       cacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                                                                                                                                               <  1:1150823/149‑1 (MQ=255)
          gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                            <  2:103329/149‑1 (MQ=255)
          gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                            <  2:1534285/149‑1 (MQ=255)
          gcgcCAGGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATAACTGCCGCTGAAGGCGCTCCAGCGTc                                                                                                                                            <  1:1034527/149‑1 (MQ=255)
                 gggTTGGAAGAATGTTTCCGTGTATGCCGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAgtaatt                                                                                                                                     <  1:1499867/149‑6 (MQ=255)
                 gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGACTGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGTCCCGGGCAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAgtaatt                                                                                                                                     <  1:1165807/149‑6 (MQ=255)
                   gTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa                                                                                                                                   <  2:1158113/149‑1 (MQ=255)
                    ttGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAAAGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGCACAGTAGCAAAGAACATAGCAATTACTGACGCTGGCGCGACTTCAGCGTCAGATACCCAc                                                                                                                                  >  2:616505/1‑149 (MQ=255)
                       taagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTcc                                                                                                                               <  1:1173808/148‑1 (MQ=255)
                       gaagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTcc                                                                                                                               >  1:888385/1‑149 (MQ=255)
                          gaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcaca                                                                                                                            <  2:1544395/149‑1 (MQ=255)
                            aTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCAcccc                                                                                                                          >  2:1581903/1‑149 (MQ=255)
                                ttCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcgcg                                                                                                                      <  1:1323920/149‑1 (MQ=255)
                                  ccGTGTAGGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTa                                                                                                                    <  2:673933/149‑1 (MQ=255)
                                   cGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaa                                                                                                                   <  2:1566150/149‑1 (MQ=255)
                                    gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTaaa                                                                                                                  >  2:1145922/1‑149 (MQ=255)
                                      gtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAgg                                                                                                                <  1:349733/149‑1 (MQ=255)
                                        aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                                              <  2:1503632/149‑1 (MQ=255)
                                         tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                                             >  2:537632/1‑149 (MQ=255)
                                         tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                                             <  1:1278503/149‑1 (MQ=255)
                                          gACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcgt                                                                                                            <  2:1029881/149‑1 (MQ=255)
                                              gcggcTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATc                                                                                                        <  1:1115343/149‑1 (MQ=255)
                                                  cTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGc                                                                                                    >  1:282007/1‑149 (MQ=255)
                                                   tGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGca                                                                                                   >  1:1286657/1‑149 (MQ=255)
                                                      tACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGACCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGACAGATACCCACTCCACACCCGCGTAAAGGACGTCATCCAGCACaa                                                                                                >  2:1070926/1‑149 (MQ=255)
                                                          aaTGGAGCAGCGATCCAATAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAATTTACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACgga                                                                                            <  2:745748/149‑3 (MQ=255)
                                                            tGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGCCGCTGGAGCGACTTCAGCGTCAGATACCCACTCCCCACCCGCGTAAAGGTCGTCATCCAGCACCACGcccc                                                                                          >  1:1238076/1‑149 (MQ=255)
                                                            tGGAGCAGCGATCCAAAAAATCAGGTAGCAACCGGTGAACGCGGCCCGTACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACAGCCGCGTAAAGGTCGTCATCCAGCACAACGCAcc                                                                                          <  1:1070926/149‑1 (MQ=255)
                                                                    cGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcc                                                                                  <  2:888385/149‑1 (MQ=255)
                                                                     taTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATcca                                                                                 <  2:331385/148‑1 (MQ=255)
                                                                      aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATACAGCACAACGCACCGAAGATccac                                                                                <  1:250073/149‑1 (MQ=255)
                                                                         cAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgca                                                                             >  1:1459581/1‑149 (MQ=255)
                                                                          aaaaaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGCCTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcac                                                                            >  2:1310562/1‑148 (MQ=255)
                                                                            aaaaTCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGGAAATAACATTGCAATTACTGACGCTGGAGCTACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGGCATCCAGCACAACGCACCGAAGATccacgcaaac                                                                          <  1:774577/149‑1 (MQ=255)
                                                                                    gTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacgcaaa                                                                  <  1:6146/149‑1 (MQ=255)
                                                                                     atgCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATGCCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCcacgcaaaccacgcaaac                                                                 <  2:794883/147‑1 (MQ=255)
                                                                                              gTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAggg                                                        >  1:646397/1‑149 (MQ=255)
                                                                                                 aaCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAAt                                                     >  1:421770/1‑149 (MQ=255)
                                                                                                   cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACACGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATaa                                                   <  1:1212040/149‑1 (MQ=255)
                                                                                                    gcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATaaa                                                  >  1:1169391/1‑149 (MQ=255)
                                                                                                        cccGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACagg                                              >  1:1482976/1‑147 (MQ=255)
                                                                                                         ccGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATcc                                             <  2:1481670/149‑1 (MQ=255)
                                                                                                          cGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATccc                                            <  2:347234/149‑1 (MQ=255)
                                                                                                           ggACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCa                                           >  2:606072/1‑149 (MQ=255)
                                                                                                                 tAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACaa                                     >  2:1306429/1‑149 (MQ=255)
                                                                                                                 tAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTAAGATACACACTACAAACACGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGCAATAAACATCCCAGCACaa                                     >  1:94433/1‑149 (MQ=255)
                                                                                                                   gCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACCTCCCAGCACAAAc                                   >  1:1188846/1‑149 (MQ=255)
                                                                                                                      aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCCCGCAAACCATGCAGGGAATAAACATCCCAGAACAAACAGt                                >  2:1352112/1‑149 (MQ=255)
                                                                                                                      aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGt                                <  1:537632/149‑1 (MQ=255)
                                                                                                                      aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGt                                >  2:369921/1‑149 (MQ=255)
                                                                                                                        aaaCATTGCAATTACTGACGCTGGTGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTaa                              <  2:784975/148‑1 (MQ=255)
                                                                                                                           cATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGAt                           >  1:909607/1‑149 (MQ=255)
                                                                                                                            aTTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATc                          >  1:1546917/1‑149 (MQ=255)
                                                                                                                                cAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGttt                      <  2:651458/149‑1 (MQ=255)
                                                                                                                                 ggtatcTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCAAGCACAACGCACCGAAGATCAACGCAACCCACGCAAACCATGCAGGGATTAAACATCCAAGCACAAACAGTAAGAGCGttta                     >  2:318812/6‑148 (MQ=255)
                                                                                                                                 aaTTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCAAAACGCACCAAAGATCCACGCAAACCACGCAAACCATCAAGGCAATAAACCTCCCAGCACAAACAGTAAGATCgagac                     >  2:62327/1‑145 (MQ=255)
                                                                                                                                  aTTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCt                    <  1:1469443/149‑1 (MQ=255)
                                                                                                                                   tatcTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtg                   >  2:647927/4‑149 (MQ=255)
                                                                                                                                      cTGACGCTGGAGCGACTTCAGCGTCAGATACCGACTGCACACCCGCGTAAAGGTCGTTATACAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCTCAAACAGTAAGATCGTTTCTGTGc                <  2:300141/149‑1 (MQ=255)
                                                                                                                                        gACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCt              >  2:1476855/1‑149 (MQ=255)
                                                                                                                                        gACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCt              >  2:678459/1‑149 (MQ=255)
                                                                                                                                         aCGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCtt             >  2:433778/1‑149 (MQ=255)
                                                                                                                                         aCGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCAGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCtt             <  1:177051/149‑1 (MQ=255)
                                                                                                                                           gCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCg           >  1:687276/1‑149 (MQ=255)
                                                                                                                                            cTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCgg          <  1:1476855/149‑1 (MQ=255)
                                                                                                                                            cTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCgg          <  1:678459/149‑1 (MQ=255)
                                                                                                                                                aGCGAGTTCAGCGTCAGATACCGACTACACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAATCCACGCAAACCATGCAGGGATTAATCATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTaa      <  2:1459581/149‑1 (MQ=255)
                                                                                                                                                aGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATACACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTaa      >  1:732579/1‑149 (MQ=255)
                                                                                                                                                    aCTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATccc  >  2:953391/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCC  >  NC_000913/1840652‑1840948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: