Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,798 C→T 6.5% S202S (TCG→TCA) ‡ 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,7980CT 6.5% 212.4 / 2.7 69S202S (TCG→TCA) ‡ynjFCDP‑alcohol phosphatidyltransferase family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/2);  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 = 2.89e-01

TTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAA  >  NC_000913/1840650‑1840944
                                                                                                                                                    |                                                                                                                                                  
ttATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGc                                                                                                                                                    <  2:1020215/149‑1 (MQ=255)
  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)
                                                                                                                                                    |                                                                                                                                                  
TTATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAA  >  NC_000913/1840650‑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: ‑