Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,377,173 A→C F83C (TTC→TGC)  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,1730AC98.5% 225.5 / ‑3.4 68F83C (TTC→TGC) 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 (1/0);  new base C (29/38);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

AGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTT  >  NC_000913/3377026‑3377319
                                                                                                                                                   |                                                                                                                                                  
aGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAAc                                                                                                                                                 <  1:332776/151‑1 (MQ=255)
   tGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGc                                                                                                                                               <  2:994008/150‑1 (MQ=255)
       tACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCa                                                                                                                                            <  1:901857/149‑1 (MQ=255)
       tACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCAt                                                                                                                                           <  1:99399/150‑1 (MQ=255)
            gCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGat                                                                                                                                     <  1:355452/151‑1 (MQ=255)
                     gCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCa                                                                                                                               <  1:92222/148‑1 (MQ=255)
                     gCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATa                                                                                                                             >  1:86890/1‑150 (MQ=255)
                       ttCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAat                                                                                                                           >  1:144562/1‑150 (MQ=255)
                        tCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATaag                                                                                                                           <  2:387231/149‑2 (MQ=255)
                         cAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAatgat                                                                                                                        <  2:738002/151‑1 (MQ=255)
                         cAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAatg                                                                                                                          <  2:1028552/149‑1 (MQ=255)
                         cAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAat                                                                                                                           >  1:614963/1‑148 (MQ=255)
                            aaGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCg                                                                                                                     <  2:642349/151‑1 (MQ=255)
                            aaGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCg                                                                                                                     >  1:785337/1‑151 (MQ=255)
                            aaGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCg                                                                                                                     >  1:518971/1‑151 (MQ=255)
                             aGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGc                                                                                                                      >  1:161776/1‑149 (MQ=255)
                              gCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCgag                                                                                                                   <  1:251867/151‑1 (MQ=255)
                               cTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCCACGCTCATCCAGATATTCCATAatgatg                                                                                                                       >  1:1121566/1‑146 (MQ=255)
                               cTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAAAGCTCATCCAGATATTCCATAATGATGCg                                                                                                                     >  1:771581/1‑148 (MQ=255)
                                 gAACCGTTGATGATGGTGTTCATCAGCGTTTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGAt                                                                                                                 <  1:656597/150‑1 (MQ=255)
                                 gAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCgaga                                                                                                                  <  1:1113149/149‑1 (MQ=255)
                                  aaCCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTTGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTc                                                                                                               <  1:68751/151‑1 (MQ=255)
                                      gTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCcac                                                                                                           >  1:435812/1‑151 (MQ=255)
                                        tgatgaTGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCcaca                                                                                                          <  2:183953/150‑1 (MQ=255)
                                         gatgatGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAgg                                                                                                        >  1:665633/1‑151 (MQ=255)
                                         gatgatGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAgg                                                                                                        >  2:218829/1‑151 (MQ=255)
                                               ggTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCCACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGc                                                                                                  >  2:306672/1‑151 (MQ=255)
                                                gtgtTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCa                                                                                                    >  1:170913/1‑148 (MQ=255)
                                                gtgtTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCt                                                                                                 >  1:829560/1‑151 (MQ=255)
                                                    tcatcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACg                                                                                             >  1:838741/1‑151 (MQ=255)
                                                     catcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACg                                                                                             <  2:476976/150‑1 (MQ=255)
                                                     catcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGa                                                                                            <  2:293019/151‑1 (MQ=255)
                                                          gCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATcca                                                                                        <  2:374720/150‑1 (MQ=255)
                                                                aCCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAggg                                                                                  >  1:1059793/1‑150 (MQ=255)
                                                                 ccAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTc                                                                                <  1:967774/151‑1 (MQ=255)
                                                                 ccAGTCTTTTTCGATGCGATGCATGTACAGACCGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGAATCGGGCAACGCTCATCCAGATATTCCATAATGCTGCGAGATCCCCACAGGGTCAGCGCACGATCCACCAgg                                                                                   >  2:414582/1‑148 (MQ=255)
                                                                    gTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGa                                                                             <  2:264256/151‑1 (MQ=255)
                                                                     tCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCgg                                                                              >  2:916465/1‑149 (MQ=255)
                                                                           tCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGctc                                                                       >  1:878420/1‑150 (MQ=255)
                                                                                gcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGAtt                                                                  <  2:727471/150‑1 (MQ=255)
                                                                                gcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTc                                                                 <  2:728263/151‑1 (MQ=255)
                                                                                  gatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCgg                                                               >  2:301837/1‑151 (MQ=255)
                                                                                        tGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCAGCAGGGTCGGAACGCTCTGATTCGGGTTg                                                           >  1:643912/1‑149 (MQ=255)
                                                                                               aCGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATc                                                  <  2:831360/151‑1 (MQ=255)
                                                                                               aCGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATc                                                  <  2:219960/151‑1 (MQ=255)
                                                                                                 ggCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAg                                                <  2:303824/151‑1 (MQ=255)
                                                                                                       cACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCt                                           <  1:157812/150‑1 (MQ=255)
                                                                                                       cACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                          <  2:198222/151‑1 (MQ=255)
                                                                                                        aCCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGa                                         <  1:141541/151‑1 (MQ=255)
                                                                                                        aCCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGa                                         <  1:881984/151‑1 (MQ=255)
                                                                                                         ccGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGGTGAGGTCAATCAGATCCTg                                          >  2:137995/1‑149 (MQ=255)
                                                                                                           gcgAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGc                                      <  2:207829/151‑1 (MQ=255)
                                                                                                             gAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                    >  1:932171/1‑151 (MQ=255)
                                                                                                              aGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGa                                   <  1:741473/151‑1 (MQ=255)
                                                                                                              aGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGa                                   <  1:252106/151‑1 (MQ=255)
                                                                                                               gCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTAATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                    >  2:330630/1‑149 (MQ=255)
                                                                                                                       gTAAACAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCttt                          >  1:304177/1‑151 (MQ=255)
                                                                                                                       gTAAACAGGCATCAGTGGCGGATGAGGGCACCGCTCATTCAGAAATTCCATAATGATGGGAGATTCCCACAGGGTCAGCTAACGAACCACCAGGGTCGGAACGCTCCGATTCGGGTTGAGGTCAATCAGAACCTGAGGCGGATTGTCCttt                          >  2:775966/1‑151 (MQ=255)
                                                                                                                          aaCAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTcc                       <  2:1025689/151‑1 (MQ=255)
                                                                                                                          aaCAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTcc                       <  1:1090532/151‑1 (MQ=255)
                                                                                                                          aaCAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGGTGAGGTCAATCAGATCCTGAGGGGGAGTGTCCtttt                         >  2:54529/1‑149 (MQ=255)
                                                                                                                            cAGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCAc                     <  2:276115/151‑1 (MQ=255)
                                                                                                                             aGGCATCAGTGGCGGATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACg                    <  2:1149401/151‑1 (MQ=255)
                                                                                                                                      tGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGAtc           >  1:278325/1‑151 (MQ=255)
                                                                                                                                          ggATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCg        >  1:762785/1‑150 (MQ=255)
                                                                                                                                           gATGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGaa      <  1:279792/151‑1 (MQ=255)
                                                                                                                                            aTGCGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGaa      >  1:714554/1‑150 (MQ=255)
                                                                                                                                               cGGGCAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt  <  1:931319/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTT  >  NC_000913/3377026‑3377319

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