Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,067,323 A→T 7.1% intergenic (‑150/+17) ygfI ← / ← yggE putative DNA‑binding transcriptional regulator/oxidative stress defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,067,3230AT 7.1% 279.3 / 7.6 102intergenic (‑150/+17)ygfI/yggEputative DNA‑binding transcriptional regulator/oxidative stress defense protein
Reads supporting (aligned to +/- strand):  ref base A (46/48);  new base T (2/6);  total (48/54)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 103 total aligned reads displayed.

CCCCCTCCCCCGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGG  >  NC_000913/3067174‑3067470
                                                                                                                                                     |                                                                                                                                                   
ccccctcccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATa                                                                                                                                                     >  2:489209/1‑150 (MQ=255)
  ccctcgcccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAACGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAACAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACg                                                                                                                                                   <  2:1213297/150‑1 (MQ=255)
   cctcccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                                                                                                  >  1:700056/1‑150 (MQ=255)
    ctcccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTAGCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAATAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                                                                                                 <  2:63883/150‑1 (MQ=255)
     tcccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCg                                                                                                                                                <  1:1079118/150‑1 (MQ=255)
      cccccGCATGCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                                                                                               <  2:883488/150‑1 (MQ=255)
       ccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                              <  1:43313/150‑1 (MQ=255)
       ccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                                                                                                              <  1:730410/150‑1 (MQ=255)
       ccccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCggg                                                                                                                                              <  2:894698/150‑1 (MQ=255)
        cccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGc                                                                                                                                             <  1:1121583/150‑1 (MQ=255)
        cccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCATTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGc                                                                                                                                             <  2:1021618/150‑1 (MQ=255)
         ccGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCt                                                                                                                                            <  2:20277/150‑1 (MQ=255)
         cagCATCCTTTTCCGCTTATCCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTTATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGATGCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCGGGCt                                                                                                                                            <  2:1037407/148‑1 (MQ=255)
            cATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACa                                                                                                                                         <  2:1348878/150‑1 (MQ=255)
              tCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTTATTACGTCATGGATTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACaaa                                                                                                                                       <  1:1588826/150‑1 (MQ=255)
              tCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACaag                                                                                                                                       <  2:773622/150‑2 (MQ=255)
                cTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACAGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGTAGGCAATATCGCGTAAAGGCATTAGCCCGACGAAAACGTCGTGCTACAAata                                                                                                                                     >  2:273644/1‑150 (MQ=255)
                   gTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtatta                                                                                                                                  >  1:1583343/1‑150 (MQ=255)
                      cGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                                                                                               >  1:1711430/1‑150 (MQ=255)
                       gCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgt                                                                                                                              >  2:1276085/1‑150 (MQ=255)
                            aCCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgctgc                                                                                                                         <  1:605383/150‑1 (MQ=255)
                              cGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAg                                                                                                                       >  2:328743/1‑150 (MQ=255)
                              cGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAg                                                                                                                       >  2:208213/1‑150 (MQ=255)
                                   cagtcagCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                                                                                                  >  2:1199299/1‑150 (MQ=255)
                                     gtcagCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                                                                                                >  2:65549/1‑150 (MQ=255)
                                     gtcagCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                                                                                                >  1:1050414/1‑150 (MQ=255)
                                         gcatATCCTTGATGAAACCGAGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAAAATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                                                                                                            >  1:712368/4‑150 (MQ=255)
                                           cTATCCTTGATGAAACCGCGAGCAAAGATCGGTGATTACGTCATGGTGTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                                                                          <  1:434906/150‑1 (MQ=255)
                                           cTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                                                                                                          >  2:34395/1‑150 (MQ=255)
                                             aTCCTTGATGAATCCGCGAGCAAAGTTAGGTGATTACGTCATGGTTTTACAGAAAATTAGAGAATAAGGAGTCAATATCGGGTAAAGGCATTAGCCCGACGAATTCGTCGGGCTACAAATATAATTGTGCTGCAGGTGTTTTAGCGGgtt                                                                                                        <  2:712368/150‑1 (MQ=255)
                                               ccTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                                                                                                      >  2:61296/1‑150 (MQ=255)
                                                cTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                                                                     <  1:34395/150‑1 (MQ=255)
                                                 ttGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                                                                    >  1:365110/1‑150 (MQ=255)
                                                 ttGATGAAACCGCGAGCAAAGATAGGTGAGTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGTCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGttgttg                                                                                                    <  2:1129375/150‑1 (MQ=255)
                                                   gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                                                                     <  1:1166306/147‑1 (MQ=255)
                                                   gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgtt                                                                                                     >  2:1166306/1‑147 (MQ=255)
                                                   gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                                                                  <  1:1018157/150‑1 (MQ=255)
                                                   gatgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                                                                                                  >  2:888899/1‑150 (MQ=255)
                                                     tgaAACCGCGAGCAAAGTTAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAATCAGGTGGCTATATCGGGTAAAGGCATTAGCCCGACGAATACGTCTGGCTACAAATTTTCTTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                                                                                <  1:87934/150‑1 (MQ=255)
                                                     tgaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                                                                                                >  2:328226/1‑150 (MQ=255)
                                                      gaAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCca                                                                                               <  1:495719/150‑1 (MQ=255)
                                                        aaCCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                                                                             <  2:1471307/150‑1 (MQ=255)
                                                        aaCCGCGAGCAAAGATAGGTGATTACGTAATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcaca                                                                                             <  1:1202483/150‑1 (MQ=255)
                                                          ccGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg                                                                                           >  1:262974/1‑150 (MQ=255)
                                                           cgcgAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGt                                                                                          >  2:829626/1‑150 (MQ=255)
                                                               aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTa                                                                                      >  1:1337763/1‑150 (MQ=255)
                                                               aGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTa                                                                                      <  2:1282236/150‑1 (MQ=255)
                                                                   aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                                                                      <  1:554399/98‑1 (MQ=255)
                                                                   aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAat                                                                                                                                      >  2:554399/1‑98 (MQ=255)
                                                                   aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                                                                  >  2:1687745/1‑150 (MQ=255)
                                                                   aaGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTg                                                                                  <  2:165532/150‑1 (MQ=255)
                                                                     gATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGa                                                                                <  2:1583343/150‑1 (MQ=255)
                                                                      aTAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAATGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGaa                                                                               <  2:1050414/150‑1 (MQ=255)
                                                                        aGGTGATTACGTCATGGTTTTACAGAATATTACAGAAAAAGGAGGCATTATCGGGTAAAGGCATTAGACCGACGAATTCGTCGGGCTACAATTATTATTGTGCTGCAGGTGTTTTTGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGa                                                                             <  1:273644/150‑1 (MQ=255)
                                                                            gATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcac                                                                         <  1:592555/150‑1 (MQ=255)
                                                                             aTTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGCCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACcaca                                                                        >  1:318428/1‑150 (MQ=255)
                                                                                aCGTCATGGGTTTACAGAAAATTACCGCAAAACGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCg                                                                     <  2:1495191/150‑1 (MQ=255)
                                                                                 cGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCgc                                                                    >  1:1302172/1‑149 (MQ=255)
                                                                                  gTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAc                                                                   <  2:731396/150‑1 (MQ=255)
                                                                                   tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                                                                  >  1:122333/1‑150 (MQ=255)
                                                                                   tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                                                                  >  2:386965/1‑150 (MQ=255)
                                                                                   tCATGGTATTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGAcc                                                                  <  1:26055/150‑1 (MQ=255)
                                                                                     aTGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTg                                                                <  1:1199299/150‑1 (MQ=255)
                                                                                      tGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGa                                                               >  2:17560/1‑150 (MQ=255)
                                                                                        gTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGatc                                                             >  2:217313/1‑150 (MQ=255)
                                                                                            tacagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGAtcatca                                                         <  1:888899/150‑1 (MQ=255)
                                                                                              cagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCaaa                                                       >  2:44010/1‑150 (MQ=255)
                                                                                               agaaagttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCACTAAGAATACCTCGTGATACAAATAAAATTGTGCAGCAGGTGTTTAAGCGGGTTGTTGATCCACAGGTTCAAACTGGAAGACAACATAGACCTGAtcatcataa                                                      >  1:1564421/1‑147 (MQ=255)
                                                                                                  aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGa                                                   <  1:372125/150‑1 (MQ=255)
                                                                                                    attacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAg                                                                              >  1:989541/1‑121 (MQ=255)
                                                                                                      tacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAg                                               >  2:563638/1‑150 (MQ=255)
                                                                                                      tacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAg                                               <  2:537437/150‑1 (MQ=255)
                                                                                                          gaaaaAGGAGGCAATATCGGGTAAAGGCATTATCCCGACGAATACGTCGGGCTACAAATATTATTGTGCGGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGATGACCACATCGACCTGATCATCAAACTGAATAGCGGc                                           >  1:219347/1‑150 (MQ=255)
                                                                                                          gaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACATGATCATCAAACTGACTAGCGGc                                           >  2:1546096/1‑150 (MQ=255)
                                                                                                              aaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGAGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCCAATGGACCTGATAAGCAAACTGAATAGCGGCCTGc                                       >  2:285132/1‑150 (MQ=255)
                                                                                                               aggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCt                                      <  2:365110/150‑1 (MQ=255)
                                                                                                                ggaggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTc                                     >  2:852045/1‑150 (MQ=255)
                                                                                                                   ggCAATATCGGGTAAAGGCATTAGCCCGAAGAAAACGTCGTTCAACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCAACAAACTGAATAGCCGCCTGCTCGTa                                  >  1:171645/1‑150 (MQ=255)
                                                                                                                     cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAg                                <  1:328743/150‑1 (MQ=255)
                                                                                                                     cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAg                                <  2:817869/150‑1 (MQ=255)
                                                                                                                     cAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAg                                <  1:208213/150‑1 (MQ=255)
                                                                                                                           cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCt                          >  1:1657511/1‑150 (MQ=255)
                                                                                                                              gTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCAAAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTggg                       >  1:971272/1‑150 (MQ=255)
                                                                                                                               tAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTTGTAAGTTTCCTGGGc                      <  2:318428/150‑1 (MQ=255)
                                                                                                                               tAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGc                      >  1:53052/1‑150 (MQ=255)
                                                                                                                               tAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGc                      <  1:829626/150‑1 (MQ=255)
                                                                                                                               tAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGc                      <  2:374328/150‑1 (MQ=255)
                                                                                                                                     cATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGacac                >  2:451023/1‑150 (MQ=255)
                                                                                                                                      aTTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACAcc               <  1:17560/150‑1 (MQ=255)
                                                                                                                                         tgCCCGACGAATACTTCGGGATACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGTTTCTAACTGTAAGAACACATCGACCTGATCATCAAACTGAATAGATGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc            <  2:262974/149‑1 (MQ=255)
                                                                                                                                         aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgc            <  2:1337763/150‑1 (MQ=255)
                                                                                                                                          tcccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTTGAAGACCACATATAACTGATCATCAAACTGAATAGAGGCCTGATAGTAAGTTTCCTGGGCGGACAACGtcg           >  1:220425/2‑150 (MQ=255)
                                                                                                                                          gCCCGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACAGgcg           >  1:472381/1‑150 (MQ=255)
                                                                                                                                          gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgcg           >  2:812351/1‑150 (MQ=255)
                                                                                                                                          gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgcg           >  2:1208598/1‑150 (MQ=255)
                                                                                                                                          gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGgcg           >  1:395393/1‑150 (MQ=255)
                                                                                                                                            ccGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGcgcg         <  1:386965/150‑1 (MQ=255)
                                                                                                                                             cgacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCgg        >  1:369968/1‑150 (MQ=255)
                                                                                                                                              gacgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGc       >  2:1457251/1‑150 (MQ=255)
                                                                                                                                                   aTACGTCGGGCTACAAATATTATTGTGCTGCAGTTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCgg  <  1:44010/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCCCTCCCCCGCATCCTTGTCCGCTTAACCGTTTCAGTCAGCCTATCCTTGATGAAACCGCGAGCAAAGATAGGTGATTACGTCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGGAAGACCACATCGACCTGATCATCAAACTGAATAGCGGCCTGCTCGTAAGTTTCCTGGGCGGACACCGGCGCGGCATCGG  >  NC_000913/3067174‑3067470

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