Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,508,438:1 +T 6.5% intergenic (‑196/+23) fryB ← / ← glk putative enzyme IIB component of PTS/glucokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,508,4381.T 6.5% 289.4 / 16.0 77intergenic (‑196/+23)fryB/glkputative enzyme IIB component of PTS/glucokinase
Reads supporting (aligned to +/- strand):  ref base . (31/41);  new base T (3/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

TGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAAT  >  NC_000913/2508289‑2508587
                                                                                                                                                      |                                                                                                                                                     
tGCAAAACAAATCGCCGTTACCAGACAAATAAGGGATTTTCTGGAGAATTTATCCCGCTGTACGGCAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTTAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAc                                                                                                                                                        <  2:54909/150‑1 (MQ=255)
 gCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAc‑c                                                                                                                                                      >  2:596126/1‑150 (MQ=255)
  cAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑Ct                                                                                                                                                     >  2:1297364/1‑150 (MQ=255)
  cAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑Ct                                                                                                                                                     >  1:401456/1‑150 (MQ=255)
   aaaaCAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCGTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGTAGAGTTAC‑CTc                                                                                                                                                    >  1:374956/1‑150 (MQ=255)
     aaCAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTccc                                                                                                                                                  <  1:623091/150‑1 (MQ=255)
     aaCAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTccc                                                                                                                                                  <  2:278645/150‑1 (MQ=255)
     aaCAAATCGCCGTTACCAGACAAATAATGCATTTAATGGAGAAATAATCCAGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAAATGAATAAACTGTATGAATTAAAAGATTATCGGGAGAGTTAC‑CTccc                                                                                                                                                  >  2:724409/1‑150 (MQ=255)
        aaaTCGCCTTAACCAGACAAATAATTCATTTACTGGAGAATTAATCCCGCTGTACGGAATTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGat                                                                                                                                               <  2:596403/150‑1 (MQ=255)
        aaaTCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAG‑GTAACTCTCCCGat                                                                                                                                               <  2:779886/150‑1 (MQ=255)
        aaaTCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAG‑GTAACTCTCCCGat                                                                                                                                               <  2:701120/150‑1 (MQ=255)
           tCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGAtata                                                                                                                                            >  1:476183/1‑150 (MQ=255)
             gCCGTTACCAGACAAATAAAGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATaaa                                                                                                                                          >  2:428045/1‑150 (MQ=255)
                   aCCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAaaggaag                                                                                                                                    <  2:779623/150‑1 (MQ=255)
                   aCCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAaaggaag                                                                                                                                    <  2:746914/150‑1 (MQ=255)
                     cAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAaggaagga                                                                                                                                  >  1:659781/1‑150 (MQ=255)
                      aGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGAt                                                                                                                                 >  2:714969/1‑150 (MQ=255)
                      aGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGAAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCTATATAAAAGGAAGGAt                                                                                                                                 >  1:1263827/1‑150 (MQ=255)
                       gACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACATAGAAAAATTGAATAAACTGTTTGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGTTATAAAAGGAAGGAtt                                                                                                                                >  2:759581/1‑150 (MQ=255)
                       gACAAAAAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTAGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGAtt                                                                                                                                >  2:445482/1‑150 (MQ=255)
                          aaaTAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTAc                                                                                                                             <  1:830517/150‑1 (MQ=255)
                            ataatGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAg                                                                                                                           <  1:242504/150‑1 (MQ=255)
                                tGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAAtg                                                                                                                       <  1:38733/150‑1 (MQ=255)
                                tGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAAtg                                                                                                                       <  2:1289610/150‑1 (MQ=255)
                                tGCATTTACTGGAGAATTAATCCAGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAAtg                                                                                                                       <  1:31935/150‑1 (MQ=255)
                                   aTTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGTGACACAAAGAAAAATTGAATAAACTGTATGTTTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGa                                                                                                                    <  1:1112492/150‑1 (MQ=255)
                                   aTTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGa                                                                                                                    <  2:1100690/150‑1 (MQ=255)
                                   aTTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGa                                                                                                                    <  2:759227/150‑1 (MQ=255)
                                   aTTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGGGAGGGGACACAAAGAAAAATTGAATACACTCTATGATTTAAAAGATTACCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGa                                                                                                                    <  1:1233902/150‑1 (MQ=255)
                                    tttACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGAc                                                                                                                   <  1:515687/150‑1 (MQ=255)
                                     atACTGGAGAATTAATCCCGCTTTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGAcc                                                                                                                  <  2:517487/149‑1 (MQ=255)
                                      tACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAATCTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCt                                                                                                                 >  1:94193/1‑150 (MQ=255)
                                      tACTGGAGAATTAATCCCGCTGTACCGAAGTGTTACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATAGAATATACTGTATGCTTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGGGACCt                                                                                                                 <  2:122363/150‑1 (MQ=255)
                                        cTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTAAAGAATGTGACCTaa                                                                                                               >  2:675406/1‑150 (MQ=255)
                                           gagaATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGt                                                                                                            <  1:1108715/150‑1 (MQ=255)
                                             taaTTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCt                                                                                                          <  2:829894/149‑1 (MQ=255)
                                                 tAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCg                                                                                                      <  1:834412/150‑1 (MQ=255)
                                                   aTCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTa                                                                                                    >  2:198396/1‑150 (MQ=255)
                                                   aTCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTa                                                                                                    <  1:1073538/150‑1 (MQ=255)
                                                   aTCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTATGGCGTa                                                                                                    >  1:408187/1‑150 (MQ=255)
                                                          tGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGc                                                                                             <  2:1060202/150‑1 (MQ=255)
                                                             aCGGAAGTGTGACGCAGGTCTACCTTGTTAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCGGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCAcc                                                                                          <  2:24939/150‑1 (MQ=255)
                                                              cGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAATCTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCAACg                                                                                         >  1:197602/1‑150 (MQ=255)
                                                              cGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCg                                                                                         >  1:1077783/1‑150 (MQ=255)
                                                                gAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGa                                                                                       >  2:498147/1‑150 (MQ=255)
                                                                    tgtgACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAATGTCTGGCGTAAATGTGCACCGGAACCg                                                                                   >  1:276659/1‑150 (MQ=255)
                                                                    tgtgACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCg                                                                                   <  2:1057512/150‑1 (MQ=255)
                                                                      tgACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCgag                                                                                 <  1:262204/150‑1 (MQ=255)
                                                                        actcAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATATAAAAGATTATCGGCAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGaa                                                                               <  2:506787/147‑1 (MQ=255)
                                                                         cGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAg                                                                              <  2:78271/150‑1 (MQ=255)
                                                                         cGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑ATCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAg                                                                              >  1:1140832/1‑150 (MQ=255)
                                                                          gCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAgg                                                                             >  2:473493/1‑150 (MQ=255)
                                                                             ggTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGccc                                                                          >  1:1075377/1‑150 (MQ=255)
                                                                              gTCGACATTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAATGGAAGGATTTACAGAATGTGACCTAAGGTGTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCg                                                                         >  2:890865/1‑150 (MQ=255)
                                                                                 gACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGAt                                                                      <  1:452603/150‑1 (MQ=255)
                                                                                 gACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGAt                                                                      >  2:682405/1‑150 (MQ=255)
                                                                                         gAGGTGACAGAAAGAAAAATTGAAATAACTGTTTGATATAAAAGATTATCGGGAGAGGTAC‑CTCCCGATATTAAATGAATGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATgg                                                              <  2:827249/150‑1 (MQ=255)
                                                                                           ggggACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGAc                                                            <  1:351608/150‑1 (MQ=255)
                                                                                                cacaAAGATATATTGAATAAACTGTATGATTTAAAAGATTAACGGGAGAGGTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGAAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGTCGATga                                                       <  2:157143/150‑1 (MQ=255)
                                                                                                cacaAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATga                                                       >  1:620318/1‑150 (MQ=255)
                                                                                                cacaAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATga                                                       <  2:620318/150‑1 (MQ=255)
                                                                                                  caAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATTTAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATgaga                                                     >  2:542474/1‑150 (MQ=255)
                                                                                                  caAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCACGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATTGACGATgaga                                                     >  2:957733/1‑150 (MQ=255)
                                                                                                   aaaGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGTAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGAt                                                    <  2:949381/150‑1 (MQ=255)
                                                                                                          aaTTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCgg                                             >  1:370402/1‑150 (MQ=255)
                                                                                                          aaTTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCgg                                             >  1:375925/1‑150 (MQ=255)
                                                                                                          aaTTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCgg                                             >  1:872131/1‑150 (MQ=255)
                                                                                                                       gTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACa                                <  2:683442/150‑1 (MQ=255)
                                                                                                                          tGATTTAAAATATTATCGGGAGATTTAC‑CTCCAGATATAAAAGGAATGATTTACAGAATGTGACCTAAGTTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTTTCATGGACGATGAGATACACCGGAATATCATGGACatat                             <  1:360314/150‑1 (MQ=255)
                                                                                                                               tAAAAGATTATCGGGATAGTTAC‑CTCCAGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCttt                        <  2:1263827/150‑1 (MQ=255)
                                                                                                                                 cttttatATCGGGAG‑GTAACTCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTaa                      >  1:357774/7‑150 (MQ=255)
                                                                                                                                 aaaGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAATTGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTaa                      >  2:1015425/1‑150 (MQ=255)
                                                                                                                                   aGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAg                    <  2:374956/150‑1 (MQ=255)
                                                                                                                                    caTTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAgc                   <  2:186132/149‑1 (MQ=255)
                                                                                                                                     tatATCGGGAG‑GTAACTCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAgcg                  >  2:884943/3‑150 (MQ=255)
                                                                                                                                       tATCGGGAG‑GTAACTCTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGcc                >  1:975534/1‑150 (MQ=255)
                                                                                                                                        aTCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGccc               <  2:415584/150‑1 (MQ=255)
                                                                                                                                        aTCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGGCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGccc               <  2:659781/150‑1 (MQ=255)
                                                                                                                                            ggAGAGTTAC‑CTCCAGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGCACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTa           <  2:476183/150‑1 (MQ=255)
                                                                                                                                                     c‑cTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAAt  <  1:1047479/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCAAAACAAATCGCCGTTACCAGACAAATAATGCATTTACTGGAGAATTAATCCCGCTGTACGGAAGTGTGACGCAGGTCGACCTTGTGAGGGGACACAAAGAAAAATTGAATAAACTGTATGATTTAAAAGATTATCGGGAGAGTTAC‑CTCCCGATATAAAAGGAAGGATTTACAGAATGTGACCTAAGGTCTGGCGTAAATGTGCACCGGAACCGAGAAGGCCCGGATTGTCATGGACGATGAGATACACCGGAATATCATGGACATATTCTTTAAAGCGCCCTTTATCTTCAAAT  >  NC_000913/2508289‑2508587

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