Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,799 G→C 7.3% S202W (TCG→TGG) ‡ 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,7990GC 7.3% 211.5 / 8.4 69S202W (TCG→TGG) ‡ynjFCDP‑alcohol phosphatidyltransferase family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/34);  new base C (2/3);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAA  >  NC_000913/1840652‑1840944
                                                                                                                                                   |                                                                                                                                                 
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:678459/1‑149 (MQ=255)
                                                                                                                                        gACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCt          >  2:1476855/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)
                                                                                                                                                   |                                                                                                                                                 
ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAA  >  NC_000913/1840652‑1840944

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: ‑