Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,840,800 A→G 7.7% S202P (TCG→CCG) ‡ ynjF ← putative phosphatidyl transferase, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,840,8000AG 7.7% 188.6 / 5.7 67S202P (TCG→CCG) ‡ynjFputative phosphatidyl transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/36);  new base G (2/3);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCC  >  NC_000913/1840652‑1840949
                                                                                                                                                    |                                                                                                                                                     
aTGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                                                                                                                                                      <  1:226189/150‑1 (MQ=255)
aTGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGAc                                                                                                                                                      <  2:670819/150‑1 (MQ=255)
 tGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACt                                                                                                                                                     <  1:1158461/150‑1 (MQ=255)
 tGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAAGGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCACATAACATTGCAATTACTGACGCTGGAGCGACt                                                                                                                                                     <  2:602987/150‑1 (MQ=255)
  gTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACtt                                                                                                                                                    >  2:1426127/1‑150 (MQ=255)
   tACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTc                                                                                                                                                   <  2:846466/150‑1 (MQ=255)
      gcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCGAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                                                                                                                                                <  1:1387801/150‑1 (MQ=255)
      gcacgcGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGc                                                                                                                                                <  1:418554/150‑1 (MQ=255)
         cgcgCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                             >  1:235248/1‑150 (MQ=255)
         cgcgCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTc                                                                                                                                             <  1:194779/150‑1 (MQ=255)
          gcgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCa                                                                                                                                            <  2:128378/150‑1 (MQ=255)
           cgcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGAAGTCGCTCCAGCGTCAg                                                                                                                                           <  1:501077/150‑1 (MQ=255)
            gcCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGa                                                                                                                                          <  2:235248/150‑1 (MQ=255)
             ccATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAATAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGAt                                                                                                                                         <  1:929060/150‑1 (MQ=255)
                 gggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATAccc                                                                                                                                     <  2:857775/150‑1 (MQ=255)
                  ggTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCa                                                                                                                                    >  2:142605/1‑150 (MQ=255)
                        aagaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCcac                                                                                                                              <  2:501922/150‑1 (MQ=255)
                         agaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                             >  1:180209/1‑134 (MQ=255)
                         agaaTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTc                                                                                                                                             <  2:180209/134‑1 (MQ=255)
                            aTGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCTGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCCCTCCACcccc                                                                                                                          >  2:1430728/1‑150 (MQ=255)
                             tGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCcg                                                                                                                         >  1:287843/1‑150 (MQ=255)
                                    gtgtATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAg                                                                                                                  <  1:1181922/150‑1 (MQ=255)
                                       tATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtc                                                                                                               >  2:247030/1‑150 (MQ=255)
                                        aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                                              <  2:287843/150‑1 (MQ=255)
                                        aTGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGgtcg                                                                                                              <  1:1046307/150‑1 (MQ=255)
                                           aCGGCGGCTGGTACGAATCGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCa                                                                                                           <  1:927463/150‑1 (MQ=255)
                                              gcggcTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATcc                                                                                                        <  2:814185/150‑1 (MQ=255)
                                              gcggcTGGTAAGAATGGAGCAGCGAGCCAAAAAATCCGGTAGCACCCGGTGAACGCGGCCCGGACAGGAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTGCAGCGTCAGATACCCACTCCACACCCTCGTAAAGGTCGTCATcc                                                                                                        <  2:1057874/150‑1 (MQ=255)
                                                        cGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATc                                                                                                         <  2:829812/139‑1 (MQ=255)
                                                        cGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATc                                                                                                         >  1:829812/1‑139 (MQ=255)
                                                        cGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGccatcgat                                                                                              >  1:629293/1‑143 (MQ=255)
                                                             ggAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGa                                                                                         >  1:1197717/1‑150 (MQ=255)
                                                                 cagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGAt                                                                                     <  2:1137367/150‑1 (MQ=255)
                                                                 cagcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGAt                                                                                     >  2:610264/1‑150 (MQ=255)
                                                                   gcGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGGCGTCATCCAGCCCAACGCACCGAAGATcc                                                                                   <  2:501881/150‑1 (MQ=255)
                                                                      aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacg                                                                                <  2:1265298/150‑1 (MQ=255)
                                                                      aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacg                                                                                >  2:486361/1‑150 (MQ=255)
                                                                      aTCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGAGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATACAGCACAACGCACCGAAGATccacg                                                                                >  1:496350/1‑150 (MQ=255)
                                                                         cAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCCCACCCGCGTAAAGGACGACATCCAGCACAACCCACCGAAGATccacgcaa                                                                             >  2:1257295/1‑150 (MQ=255)
                                                                                tCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacgc                                                                      >  2:552484/1‑150 (MQ=255)
                                                                                tCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATccacgcaaaccacgc                                                                      <  1:864183/150‑1 (MQ=255)
                                                                                    gTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCcacgcaaaccacgcaaac                                                                  >  2:178104/1‑150 (MQ=255)
                                                                                        cAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATg                                                              >  2:450644/1‑150 (MQ=255)
                                                                                           ccGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTGAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAg                                                           <  1:1106198/150‑1 (MQ=255)
                                                                                             ggTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAggg                                                         >  1:234527/1‑150 (MQ=255)
                                                                                               tGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGaa                                                       <  1:142605/150‑1 (MQ=255)
                                                                                                gAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAAt                                                      <  1:1426127/150‑1 (MQ=255)
                                                                                                  aCGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATaa                                                    <  2:486596/150‑1 (MQ=255)
                                                                                                   cgcgGCCCGGACAGTAGCAAATAACATTGCAATTACTGACCCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATaaa                                                   >  2:971472/1‑150 (MQ=255)
                                                                                                      ggCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACAt                                                >  2:389294/1‑150 (MQ=255)
                                                                                                         ccGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATccc                                             >  2:1018980/1‑150 (MQ=255)
                                                                                                          cGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCa                                            <  2:142478/150‑1 (MQ=255)
                                                                                                             aCAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGca                                         <  1:1275036/150‑1 (MQ=255)
                                                                                                               aGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGcaca                                       <  1:486361/150‑1 (MQ=255)
                                                                                                                gTAGCAAATAACATTGCAATTACTAACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACaa                                      >  2:1279189/1‑150 (MQ=255)
                                                                                                                      aaTAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTa                                <  1:555215/150‑1 (MQ=255)
                                                                                                                            aTTGCAATTACTGACGCTGGAGCGCCTTCAGCGTCAGATACCCACTCCACCCCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCg                          >  2:1356356/1‑150 (MQ=255)
                                                                                                                             ttGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGt                         >  1:479467/1‑150 (MQ=255)
                                                                                                                                    tACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtgtg                  <  1:450644/150‑1 (MQ=255)
                                                                                                                                    tACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCtgtg                  >  2:83776/1‑150 (MQ=255)
                                                                                                                                     aCTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGc                 <  1:1257295/150‑1 (MQ=255)
                                                                                                                                      cTGACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGcc                >  1:1231148/1‑150 (MQ=255)
                                                                                                                                        gACGCTGAAGTCGCTCCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGCTCGTTTCTGTGCCtt              >  1:421532/1‑150 (MQ=255)
                                                                                                                                         aCGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTc             <  1:1279189/150‑1 (MQ=255)
                                                                                                                                             tGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGtt         <  1:552484/150‑1 (MQ=255)
                                                                                                                                                aGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAAt      >  2:487234/1‑150 (MQ=255)
                                                                                                                                                aGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAAt      <  2:479467/150‑1 (MQ=255)
                                                                                                                                                    aCTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATcccc  <  1:247030/150‑1 (MQ=255)
                                                                                                                                                    aCTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACAGAAGATCCACGCAACCCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATcccc  >  2:50160/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGTACGCACGCGCCATGGGTTGGAAGAATGTTTCCGTGTATGACGGCGGCTGGTACGAATGGAGCAGCGATCCAAAAAATCCGGTAGCAACCGGTGAACGCGGCCCGGACAGTAGCAAATAACATTGCAATTACTGACGCTGGAGCGACTTCAGCGTCAGATACCCACTCCACACCCGCGTAAAGGTCGTCATCCAGCACAACGCACCGAAGATCCACGCAAACCACGCAAACCATGCAGGGAATAAACATCCCAGCACAAACAGTAAGATCGTTTCTGTGCCTTCGGTTAATCCCC  >  NC_000913/1840652‑1840949

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