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,2150AG99.1% 373.9 / ‑5.3 109L69P (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 (1/0);  new base G (53/55);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 113 total aligned reads displayed.

GATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTAT  >  NC_000913/3377067‑3377359
                                                                                                                                                    |                                                                                                                                                
gatgatGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACggg                                                                                                                                                >  3:56169/1‑151 (MQ=255)
 atgatgGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACgggg                                                                                                                                               <  3:205663/151‑1 (MQ=255)
  tgatgGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACgggg                                                                                                                                               >  4:1647480/1‑150 (MQ=255)
     tgGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTc                                                                                                                                             >  4:1297532/1‑149 (MQ=255)
        tgtTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTc                                                                                                                                        >  2:1085473/1‑151 (MQ=255)
         gtTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCCTCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTc                                                                                                                                        >  2:607149/1‑150 (MQ=255)
          ttCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCAc                                                                                                                                      >  4:1098908/1‑151 (MQ=255)
           tcatcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACg                                                                                                                                     >  2:653733/1‑151 (MQ=255)
            catcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGa                                                                                                                                    <  4:1153884/151‑1 (MQ=255)
             atcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGa                                                                                                                                    >  2:468478/1‑150 (MQ=255)
              tcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGGAACGCGCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGa                                                                                                                                    >  4:1298017/1‑149 (MQ=255)
                aGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATcca                                                                                                                                <  1:88438/151‑1 (MQ=255)
                  cGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATccacc                                                                                                                              <  4:1089780/151‑1 (MQ=255)
                  cGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATccacc                                                                                                                              >  4:2880/1‑151 (MQ=255)
                     gtACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAg                                                                                                                            <  2:1577045/150‑1 (MQ=255)
                     gtACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAg                                                                                                                            <  4:297029/150‑1 (MQ=255)
                      tACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAggg                                                                                                                          >  2:211129/1‑151 (MQ=255)
                       aCCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAggg                                                                                                                          <  4:732209/150‑1 (MQ=255)
                       aCCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAggg                                                                                                                          <  4:531561/150‑1 (MQ=255)
                       aCCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGt                                                                                                                         <  2:755612/151‑1 (MQ=255)
                         cAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCg                                                                                                                       >  4:311987/1‑151 (MQ=255)
                          aGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCgg                                                                                                                      <  2:1373323/151‑1 (MQ=255)
                             cTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGaa                                                                                                                    >  1:279015/1‑150 (MQ=255)
                             cTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGaa                                                                                                                    >  2:298338/1‑150 (MQ=255)
                              tttttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCCGCTCACGATCCACCAGGGTCGGAACg                                                                                                                  >  4:793315/1‑151 (MQ=255)
                               ttttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGaa                                                                                                                    <  2:230320/148‑1 (MQ=255)
                                tttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGc                                                                                                                 >  4:1186981/1‑150 (MQ=255)
                                 ttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGct                                                                                                                <  4:315157/150‑1 (MQ=255)
                                       gcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGAtt                                                                                                          >  2:203088/1‑150 (MQ=255)
                                       gcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGAtt                                                                                                          >  2:789149/1‑150 (MQ=255)
                                        cgatgcATGTACAGACGGCTTTCACAGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCAACCAGATATTCCATAATGATGGGGGATTCCCACGGGGTCAACTCACGATCCACCAGGGTCGGAACGCTCTGATTCg                                                                                                        >  2:797317/1‑151 (MQ=255)
                                         gatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCCGGGTCGGAACGCTCTGATTCgg                                                                                                       >  1:135222/1‑151 (MQ=255)
                                              aTGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGa                                                                                                  <  2:380192/151‑1 (MQ=255)
                                                   cAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCGCGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCa                                                                                             <  2:1229017/151‑1 (MQ=255)
                                                   cAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCa                                                                                             <  4:466859/151‑1 (MQ=255)
                                                   cAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACAATCCACCAGGGTCGGGACGCTCTGATTCGGGTTTAGGTc                                                                                              >  2:1334708/1‑150 (MQ=255)
                                                    aGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAgg                                                                                                >  2:1111422/1‑147 (MQ=255)
                                                       cGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCa                                                                                         <  3:404942/151‑1 (MQ=255)
                                                       cGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCa                                                                                         <  4:1487290/151‑1 (MQ=255)
                                                        ggCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATc                                                                                          >  1:111699/1‑149 (MQ=255)
                                                         gCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCCCAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGGCAATc                                                                                          >  2:1014005/1‑148 (MQ=255)
                                                            ttCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATcc                                                                                    >  3:439921/1‑151 (MQ=255)
                                                              cACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                                                                  <  3:440679/151‑1 (MQ=255)
                                                              cACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                                                                  <  2:1379151/151‑1 (MQ=255)
                                                               aCCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                                                                  >  4:1075958/1‑150 (MQ=255)
                                                                ccGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAg                                                                                <  2:382638/151‑1 (MQ=255)
                                                                 cgcgAGCTACCGGGTAAACAGGCATCAGTGGCGGAGGCGGGAAACGCTCATCCCGGGATTCCATAATGATGCGAGATTACCACGGGGTCAGCTCACGATCCACCCGGGTCGGAACGCTCCGATTCGGGGTGAGGGCAATCAGATCCTGAgg                                                                               >  1:247333/1‑151 (MQ=255)
                                                                  gcgAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAgg                                                                               >  4:1287675/1‑150 (MQ=255)
                                                                  gcgAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGc                                                                              <  2:1097513/151‑1 (MQ=255)
                                                                   cgAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCg                                                                             <  2:514906/151‑1 (MQ=255)
                                                                    gAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                                                            <  2:311516/151‑1 (MQ=255)
                                                                    gAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                                                            <  2:521812/151‑1 (MQ=255)
                                                                     aGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGa                                                                           <  2:1295897/151‑1 (MQ=255)
                                                                      gCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAt                                                                          <  4:922801/151‑1 (MQ=255)
                                                                        tACCGGGTAAACAGGCATCAGTGGCGGATTCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGGTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                         >  4:1492737/1‑150 (MQ=255)
                                                                        tACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                         >  2:586839/1‑150 (MQ=255)
                                                                         aCCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                         <  4:1178033/149‑1 (MQ=255)
                                                                         aCCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTg                                                                        <  2:568848/150‑1 (MQ=255)
                                                                                 aaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTcc                                                               <  3:418171/151‑1 (MQ=255)
                                                                                 aaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTcc                                                               >  4:615671/1‑151 (MQ=255)
                                                                                   cAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCAc                                                             <  4:211469/151‑1 (MQ=255)
                                                                                    aGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACg                                                            <  1:46603/151‑1 (MQ=255)
                                                                                      gCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgtg                                                          <  4:1086852/151‑1 (MQ=255)
                                                                                       cATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACgtgt                                                         <  4:257344/151‑1 (MQ=255)
                                                                                          cAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTc                                                       <  2:560789/150‑1 (MQ=255)
                                                                                          cAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTc                                                       <  2:301284/150‑1 (MQ=255)
                                                                                           aGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGa                                                     <  4:1447283/151‑1 (MQ=255)
                                                                                            gTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGGGTTCGAt                                                    >  2:343916/1‑151 (MQ=255)
                                                                                               gCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGAtctc                                                 <  4:1514179/151‑1 (MQ=255)
                                                                                                cGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCg                                                >  4:215338/1‑151 (MQ=255)
                                                                                                 ggATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCg                                                <  3:161311/150‑1 (MQ=255)
                                                                                                 ggATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCg                                                >  1:161297/1‑150 (MQ=255)
                                                                                                 ggATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCg                                                >  2:438901/1‑150 (MQ=255)
                                                                                                     gCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACt                                           <  4:547976/151‑1 (MQ=255)
                                                                                                      cGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt                                          <  2:1635167/151‑1 (MQ=255)
                                                                                                       gggAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACtt                                          >  2:970999/1‑150 (MQ=255)
                                                                                                           aaCGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACAcc                                     <  3:111713/151‑1 (MQ=255)
                                                                                                            aCGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCt                                    >  1:160051/1‑151 (MQ=255)
                                                                                                               cTCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACGCCttt                                  <  2:476067/150‑1 (MQ=255)
                                                                                                                tCATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtct                                >  1:172588/1‑151 (MQ=255)
                                                                                                                 cATCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtctc                               >  2:111648/1‑151 (MQ=255)
                                                                                                                  aTCCAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTtctc                               >  4:927157/1‑150 (MQ=255)
                                                                                                                    ccAGATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagc                            <  4:469131/151‑1 (MQ=255)
                                                                                                                       gATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagccag                         <  1:94943/151‑1 (MQ=255)
                                                                                                                       gATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagcc                           >  4:1561692/1‑149 (MQ=255)
                                                                                                                       gATATTCCATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCGCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagccag                         <  2:1570596/151‑1 (MQ=255)
                                                                                                                             ccATAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAAt                   <  2:865914/151‑1 (MQ=255)
                                                                                                                               aTAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGc                 <  2:721534/151‑1 (MQ=255)
                                                                                                                               aTAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGc                 >  4:142161/1‑151 (MQ=255)
                                                                                                                               aTAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGc                 <  4:867439/151‑1 (MQ=255)
                                                                                                                                tAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGGTCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCg                >  4:160728/1‑151 (MQ=255)
                                                                                                                                tAATGATGCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCg                >  2:685053/1‑151 (MQ=255)
                                                                                                                                  atgatgCGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGa              <  4:1630703/151‑1 (MQ=255)
                                                                                                                                        cGAGATTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGAt        <  3:263108/151‑1 (MQ=255)
                                                                                                                                         gagaTTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATg       >  4:1428923/1‑151 (MQ=255)
                                                                                                                                         gagaTTCCCACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATg       >  1:313868/1‑151 (MQ=255)
                                                                                                                                               cccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGACACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat  >  4:653338/1‑150 (MQ=255)
                                                                                                                                               cccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat  <  3:51771/150‑1 (MQ=255)
                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat  >  2:1328058/1‑149 (MQ=255)
                                                                                                                                                ccACGGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCtat  >  2:990368/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTAT  >  NC_000913/3377067‑3377359

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: