Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,803 T→C 7.2% K201E (AAG→GAG)  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,8030TC 7.2% 198.1 / 5.4 69K201E (AAG→GAG) ynjFCDP‑alcohol phosphatidyltransferase family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base C (2/3);  total (32/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 = 2.43e-01

CACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCC  >  NC_000913/1840659‑1840948
                                                                                                                                                |                                                                                                                                                 
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                                                                                                             <  1:1278503/149‑1 (MQ=255)
                                  tGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                                             >  2:537632/1‑149 (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)
                                                                                                                                                |                                                                                                                                                 
CACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCC  >  NC_000913/1840659‑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: ‑