Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,377,215 A→G L69P (CTG→CCG)  sspA ← stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,377,2150AG98.7% 262.3 / NA 74L69P (CTG→CCG) sspAstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/41);  minor base C (1/0);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATA  >  NC_000913/3377067‑3377360
                                                                                                                                                    |                                                                                                                                                 
gatgatGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACggg                                                                                                                                                 <  1:909189/151‑1 (MQ=255)
    atgGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGGAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTc                                                                                                                                              >  2:1369820/1‑150 (MQ=255)
      ggTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGc                                                                                                                                           <  1:723131/151‑1 (MQ=255)
      ggTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGATCATCCAGATATTCCATACTGATGCGAGATTCCCACgggggc                                                                                                                                              >  2:848588/1‑146 (MQ=255)
       gtgtTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCt                                                                                                                                          <  2:716294/151‑1 (MQ=255)
           tcatcaGCGTGTACCAGACTTTTTCGATGCGATGCGTGTACAGACGGCTTTCAACGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGTAACGCTCAACCAAATATTCCGTAATGATACGAAATTCCCCCCGGGACAGCGCAc                                                                                                                                       >  2:889580/1‑150 (MQ=255)
            catcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGGCAGCTCAc                                                                                                                                       >  2:325137/1‑149 (MQ=255)
             atcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGAt                                                                                                                                    >  1:1110577/1‑151 (MQ=255)
                   gtgtACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATccacca                                                                                                                              <  1:493115/151‑1 (MQ=255)
                    tgtACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATccacc                                                                                                                               >  1:1352511/1‑149 (MQ=255)
                     gtACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAgg                                                                                                                            <  2:881857/151‑1 (MQ=255)
                       aCCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGt                                                                                                                          <  2:1442201/151‑1 (MQ=255)
                         cAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGGCg                                                                                                                        >  1:1203455/1‑151 (MQ=255)
                               ttttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGc                                                                                                                  >  1:1270671/1‑151 (MQ=255)
                                tttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGc                                                                                                                  <  2:348990/150‑1 (MQ=255)
                                 ttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGGCAGCTCACGATCCACCAGGGTCGGAACGctc                                                                                                                >  1:1067116/1‑151 (MQ=255)
                                     atgcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGAt                                                                                                            >  2:799210/1‑151 (MQ=255)
                                      tgcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGAtt                                                                                                           >  2:1135687/1‑151 (MQ=255)
                                       gcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTc                                                                                                          <  1:605837/151‑1 (MQ=255)
                                           tgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCggg                                                                                                       >  1:1268717/1‑150 (MQ=255)
                                            gcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGtt                                                                                                     <  2:281577/151‑1 (MQ=255)
                                            gcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGt                                                                                                      <  2:49095/150‑1 (MQ=255)
                                                   cAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTc                                                                                               >  2:757181/1‑150 (MQ=255)
                                                       cGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCa                                                                                          <  2:1326913/151‑1 (MQ=255)
                                                         gCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGa                                                                                        <  2:668415/151‑1 (MQ=255)
                                                             tCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATcc                                                                                     <  2:934537/150‑1 (MQ=255)
                                                              cACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                                                                   >  2:754241/1‑151 (MQ=255)
                                                               aCCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGa                                                                                  <  2:308329/151‑1 (MQ=255)
                                                                   cgAGCTACCGGGTAAACAGGCATCAGAGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGCTCCTGAgg                                                                                >  2:1096519/1‑149 (MQ=255)
                                                                    gAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                                                             <  1:771368/151‑1 (MQ=255)
                                                                        tACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                          >  1:1303489/1‑150 (MQ=255)
                                                                        tACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                          >  1:798746/1‑150 (MQ=255)
                                                                          ccGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGt                                                                        >  1:134878/1‑150 (MQ=255)
                                                                          ccGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTc                                                                       <  1:872912/151‑1 (MQ=255)
                                                                              gTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCttt                                                                   <  2:635568/151‑1 (MQ=255)
                                                                              gTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGACCttt                                                                   >  1:1174965/1‑151 (MQ=255)
                                                                                aaaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTc                                                                 <  2:1393428/151‑1 (MQ=255)
                                                                                  aCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCa                                                               <  1:518919/151‑1 (MQ=255)
                                                                                  aCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCa                                                               <  2:699748/151‑1 (MQ=255)
                                                                                     ggCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgt                                                            <  2:795210/151‑1 (MQ=255)
                                                                                      gCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgtg                                                           <  1:1369924/151‑1 (MQ=255)
                                                                                       cATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgtgt                                                          >  1:853609/1‑151 (MQ=255)
                                                                                         tCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGtt                                                         >  1:495859/1‑150 (MQ=255)
                                                                                           aGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGa                                                      <  1:1339127/151‑1 (MQ=255)
                                                                                                      cGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt                                           >  2:1090296/1‑151 (MQ=255)
                                                                                                      cGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt                                           <  1:544942/151‑1 (MQ=255)
                                                                                                       gggAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt                                           >  1:82243/1‑150 (MQ=255)
                                                                                                       gggAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTa                                          <  2:1411867/151‑1 (MQ=255)
                                                                                                          aaaCGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTacac                                       <  2:1194137/151‑1 (MQ=255)
                                                                                                          aaaCGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTacac                                       >  1:982626/1‑151 (MQ=255)
                                                                                                             cGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCtt                                    >  1:1260170/1‑151 (MQ=255)
                                                                                                              gCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCttt                                   >  1:498696/1‑151 (MQ=255)
                                                                                                              gCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCttt                                   <  2:427669/151‑1 (MQ=255)
                                                                                                               cTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtc                                  >  1:1341450/1‑151 (MQ=255)
                                                                                                               cTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtc                                  >  1:85637/1‑151 (MQ=255)
                                                                                                                 cATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtc                                  <  1:153331/149‑1 (MQ=255)
                                                                                                                     cAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagcc                            <  2:209812/151‑1 (MQ=255)
                                                                                                                      aGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagcca                           <  2:262359/151‑1 (MQ=255)
                                                                                                                       gATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagccag                          <  1:984994/151‑1 (MQ=255)
                                                                                                                         tatTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagccagc                         >  1:842753/1‑150 (MQ=255)
                                                                                                                          atTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGcac                       <  2:273740/151‑1 (MQ=255)
                                                                                                                          atTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGca                        <  1:485874/150‑1 (MQ=255)
                                                                                                                             ccATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAAt                    <  1:886062/151‑1 (MQ=255)
                                                                                                                              cATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATg                   <  2:1424598/151‑1 (MQ=255)
                                                                                                                                  atgatgCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGa               <  2:915252/151‑1 (MQ=255)
                                                                                                                                         gagaTTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATg        <  1:848674/151‑1 (MQ=255)
                                                                                                                                         gagaTTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATg        <  2:288853/151‑1 (MQ=255)
                                                                                                                                             ttCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGCTGAGGTCAATCAGCTCCTGAGGCGGAATGTCCTTTTCCACGTGTTCGGTCACGAGACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGAttgat     >  2:349145/1‑146 (MQ=255)
                                                                                                                                              tCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCta    <  1:889668/150‑1 (MQ=255)
                                                                                                                                               cccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtata  <  1:97482/151‑1 (MQ=255)
                                                                                                                                               cccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtata  <  2:798668/151‑1 (MQ=255)
                                                                                                                                               cccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat   >  1:20518/1‑150 (MQ=255)
                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGAc              <  1:397897/138‑1 (MQ=255)
                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGAc              >  2:397841/1‑138 (MQ=255)
                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat   >  2:248823/1‑149 (MQ=255)
                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat   >  1:817206/1‑149 (MQ=255)
                                                                                                                                                 cacgGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat   >  2:697037/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATA  >  NC_000913/3377067‑3377360

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