Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,966 T→A 7.3% intergenic (‑39/+18) ralR ← / ← recT Rac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,9660TA 7.3% 324.6 / 11.6 110intergenic (‑39/+18)ralR/recTRac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base T (38/64);  new base A (2/6);  total (40/70)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 112 total aligned reads displayed.

GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATATTT  >  NC_000913/1413818‑1414115
                                                                                                                                                    |                                                                                                                                                     
gTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGTTGCAtt                                                                                                                                                      >  2:1143458/1‑150 (MQ=255)
      gATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCg                                                                                                                                                >  1:1080774/1‑150 (MQ=255)
          cGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaa                                                                                                                                            <  1:683447/150‑2 (MQ=255)
          cGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATAATTCAAGTTGGCGGTGCATTACACCGCCAg                                                                                                                                            >  1:1109723/1‑150 (MQ=255)
           gCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaac                                                                                                                                           <  1:1655847/150‑3 (MQ=255)
            cATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                                                                          <  1:1367182/150‑4 (MQ=255)
             aTCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaactt                                                                                                                                         <  2:147278/150‑5 (MQ=255)
               ccGTTACACTTCGCTCGGTAATATCCTGTAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCAAGGCTGa                                                                                                                                       >  2:1342366/1‑150 (MQ=255)
                cGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTTGCGGTGCATTACACCGCCAGGCTGaa                                                                                                                                      >  1:1269115/1‑150 (MQ=255)
                cGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGaa                                                                                                                                      >  1:805001/1‑150 (MQ=255)
                 tttACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                                                                     <  1:1217270/149‑1 (MQ=255)
                    acacTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattat                                                                                                                                  <  1:1394987/150‑1 (MQ=255)
                          cGCTCGGTAATATCCTGAAATGGCTTTCACCATTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctct                                                                                                                            <  2:1310009/150‑1 (MQ=255)
                             tatttAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATTGACATGGTATAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGaa                                                                                                                         <  2:187784/146‑1 (MQ=255)
                               ggTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAtt                                                                                                                       <  2:160395/150‑1 (MQ=255)
                                  aaTATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTTGCGGTGCATTACACCGCCAGGCTGAATTAATCCTCTGAATTATc                                                                                                                    <  2:1623509/150‑1 (MQ=255)
                                   ttataCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACTCCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                                                                   <  1:1071815/149‑1 (MQ=255)
                                   atatCCTGAAATGGCTTTCACCGTTACTGATGGAAAACCACAAAATGGACATGGTTTAAGATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                                                                   <  1:1052715/150‑1 (MQ=255)
                                     atCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAt                                                                                                                 >  2:1817245/1‑150 (MQ=255)
                                      tCCTGAAATGGCTTTCACCGTTGCTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAtt                                                                                                                >  1:519835/1‑150 (MQ=255)
                                        cTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTaa                                                                                                              >  2:1115559/1‑149 (MQ=255)
                                           aaaTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACt                                                                                                           <  1:810818/150‑1 (MQ=255)
                                             aTGGCTTTCACCGTTACTGATTGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGt                                                                                                         <  2:1043096/150‑1 (MQ=255)
                                              tGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTa                                                                                                        <  2:1532486/150‑1 (MQ=255)
                                               ggCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAt                                                                                                       <  2:1755260/150‑1 (MQ=255)
                                               ggCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAt                                                                                                       >  2:1177914/1‑150 (MQ=255)
                                                  tttCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcc                                                                                                    >  1:433583/1‑150 (MQ=255)
                                                     cACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCg                                                                                                 >  2:415448/1‑150 (MQ=255)
                                                        cGTTACTGATGGAAAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                                                                              <  2:653807/150‑1 (MQ=255)
                                                         gTTACTGATGGACAACCACAAAATGGACATGGTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTaa                                                                                            <  1:1002524/150‑1 (MQ=255)
                                                           tACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTAAt                                                                                           <  1:1170466/150‑1 (MQ=255)
                                                           tACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAt                                                                                           >  1:1434107/1‑150 (MQ=255)
                                                            aCTGATGGACAACCACAAAATGTACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATa                                                                                          <  2:1261981/150‑1 (MQ=255)
                                                            aCTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATa                                                                                          <  2:214183/150‑1 (MQ=255)
                                                             cTGATGGACAACCACAAAATGGACATGGTTTAATATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAc                                                                                         <  2:511300/150‑1 (MQ=255)
                                                             cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAc                                                                                         <  1:1170465/150‑1 (MQ=255)
                                                             cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAc                                                                                         <  2:582132/150‑1 (MQ=255)
                                                             cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAc                                                                                         <  2:164764/150‑1 (MQ=255)
                                                              tGATGGACAACCACAAAATGGACATGGTTTATCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACa                                                                                        <  1:1809783/150‑1 (MQ=255)
                                                               gATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACag                                                                                       >  2:42112/1‑150 (MQ=255)
                                                               gATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACag                                                                                       >  2:1473824/1‑150 (MQ=255)
                                                                 tGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                                                                     <  1:595321/150‑1 (MQ=255)
                                                                  ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                                                                    >  2:1107314/1‑150 (MQ=255)
                                                                  ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                                                                    >  2:43510/1‑150 (MQ=255)
                                                                  ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                                                                    <  1:614658/150‑1 (MQ=255)
                                                                     cAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAAt                                                                                 >  2:701797/1‑150 (MQ=255)
                                                                      aaCCACAAAATGGACATTGTTTAACATTGTCATATCTATTTATTTTTCTCATAACAAATATTTCAAGTTGTCGTTGCATTACACCGCCAGGCTGACTTATTCCTCTGAATTATCGATTACACTGTATTAACCGGGGAATACAGTGTAATc                                                                                <  2:242511/150‑1 (MQ=255)
                                                                       aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATAGATTACACTGTATTCCCCGGTTAATACAGAGGAATCt                                                                               >  1:646285/1‑150 (MQ=255)
                                                                         cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAACGCCAGGCTGAATTATTCCTATGAATTATAGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                                                                             >  2:1239631/1‑150 (MQ=255)
                                                                               aTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                                                                       <  1:1610132/150‑1 (MQ=255)
                                                                               aTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                                                                       <  1:628947/150‑1 (MQ=255)
                                                                               aTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCATCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                                                                       >  1:354581/1‑150 (MQ=255)
                                                                                 ggACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGACAGGCTGAATTATTCATCTGAATTATCGATTACACTGTATTCACCGGTTAATACAGAGTAATCTGCAGGATCGa                                                                     >  1:192902/1‑150 (MQ=255)
                                                                                  gACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTTTCGATTACACTGTATTCACCGGGTAATACAGAGGAATCTGCAGGATCGAt                                                                    >  2:1230328/1‑150 (MQ=255)
                                                                                    cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAACACAGAGGAATCTGCAGGATCGATTg                                                                  <  2:1288196/150‑1 (MQ=255)
                                                                                    cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGCTAATACAGAGGAATCTGCAGGATCGATTg                                                                  <  1:112735/150‑1 (MQ=255)
                                                                                      tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCTAGTTGGCGGTGCATTACACCGCAAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc                                                                <  1:1097955/150‑1 (MQ=255)
                                                                                       ggTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACAGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAATGAGGAATCTGCAGGATCGATTGTCa                                                               >  2:44178/1‑150 (MQ=255)
                                                                                        gTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCTCCAGGCTGAATTATTACTCTGAATTATCGATTACACTGTATTAACCGGTTAATACAGAGGAATCTGCAGGATCGATTGTAAg                                                              >  1:1356064/1‑150 (MQ=255)
                                                                                             aCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCTGTGCATTACACCGCCTGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGtt                                                         <  1:1570797/150‑1 (MQ=255)
                                                                                               aTTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTcc                                                       <  1:1352985/150‑1 (MQ=255)
                                                                                                  ttCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt                                                    <  2:1191866/149‑1 (MQ=255)
                                                                                                  gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt                                                    <  1:410593/150‑1 (MQ=255)
                                                                                                  gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGGTAATACAGAGGAATCTGCAGTATCGATTGTCAGTGGTTCCttt                                                    >  1:1261967/1‑150 (MQ=255)
                                                                                                   tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtttt                                                   <  2:1080774/150‑1 (MQ=255)
                                                                                                   tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtttt                                                   <  1:1350833/150‑1 (MQ=255)
                                                                                                   tCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtttt                                                   <  1:1021287/150‑1 (MQ=255)
                                                                                                    cATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTc                                                  <  1:1238890/150‑1 (MQ=255)
                                                                                                    cATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTc                                                  <  2:1654007/150‑1 (MQ=255)
                                                                                                    cATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTc                                                  <  1:221238/150‑1 (MQ=255)
                                                                                                     atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCa                                                 <  2:81606/150‑1 (MQ=255)
                                                                                                     atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCa                                                 <  1:368759/150‑1 (MQ=255)
                                                                                                     atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCa                                                 <  1:225467/150‑1 (MQ=255)
                                                                                                      tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCAGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCAt                                                >  1:441498/1‑150 (MQ=255)
                                                                                                       atCTAATAATTTTTCTCATAAAAAATATTTAAAGTTGGAGGTGCATTACAACTCCATGCTGAATTAATAATCTGAATTATCGATTACACTGTATTCACCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTAATTTTAATc                                               >  2:489649/1‑150 (MQ=255)
                                                                                                        tctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATcc                                              <  1:991573/150‑1 (MQ=255)
                                                                                                            aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTg                                          <  1:1308610/150‑1 (MQ=255)
                                                                                                             tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGa                                         <  2:449893/150‑1 (MQ=255)
                                                                                                              aaTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGAt                                        >  2:79/1‑150 (MQ=255)
                                                                                                               aTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTTAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATa                                       <  2:67455/150‑1 (MQ=255)
                                                                                                               aTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATa                                       >  2:1021729/1‑150 (MQ=255)
                                                                                                                    tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCa                                  <  2:754098/150‑1 (MQ=255)
                                                                                                                        aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCCTCCATTGATACTGCACGCt                              >  1:25780/1‑150 (MQ=255)
                                                                                                                        aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCt                              >  2:709651/1‑150 (MQ=255)
                                                                                                                        aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTACCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCt                              >  1:1599876/1‑150 (MQ=255)
                                                                                                                         tAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTg                             >  2:1478279/1‑150 (MQ=255)
                                                                                                                            aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAt                          <  2:1427369/150‑1 (MQ=255)
                                                                                                                             aaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGATGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGGTCCTTTTCATCCATTGATACTGCACGCTGGAtc                         <  1:352406/150‑1 (MQ=255)
                                                                                                                               atatTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAtctc                       <  1:1817245/150‑1 (MQ=255)
                                                                                                                                tatTTCAAGTTGGCGTTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGTAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCa                      <  1:836412/150‑1 (MQ=255)
                                                                                                                                  tttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAAt                    >  2:1068532/1‑150 (MQ=255)
                                                                                                                                   ttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAAtt                   <  1:415448/150‑1 (MQ=255)
                                                                                                                                    tCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGTAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTg                  <  1:1342366/150‑1 (MQ=255)
                                                                                                                                    tCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTg                  <  2:147243/150‑1 (MQ=255)
                                                                                                                                        cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATAc              >  1:1751733/3‑150 (MQ=255)
                                                                                                                                        cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATAc              >  1:493593/3‑150 (MQ=255)
                                                                                                                                           ggCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATCCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTGCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACggg           <  1:1818689/150‑1 (MQ=255)
                                                                                                                                            gCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGc          >  1:1046054/1‑150 (MQ=255)
                                                                                                                                                 gCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAata     <  1:1017269/150‑1 (MQ=255)
                                                                                                                                                    ttACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATAttt  <  1:1368006/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGCAAATATTT  >  NC_000913/1413818‑1414115

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: