Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,187,090 G→C 7.1% intergenic (+20/+29) pepT → / ← roxA peptidase T/ribosomal protein‑arginine oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,187,0900GC 7.1% 287.5 / 10.1 85intergenic (+20/+29)pepT/roxApeptidase T/ribosomal protein‑arginine oxygenase
Reads supporting (aligned to +/- strand):  ref base G (43/36);  new base C (4/2);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

ACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAG  >  NC_000913/1186941‑1187235
                                                                                                                                                     |                                                                                                                                                 
aCCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtg                                                                                                                                                   >  2:1343556/1‑150 (MQ=255)
 ccGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                                                                  <  2:601405/150‑1 (MQ=255)
 ccGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCAtgt                                                                                                                                                  >  1:557024/1‑150 (MQ=255)
   gTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGa                                                                                                                                                >  1:171825/1‑150 (MQ=255)
     gCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATg                                                                                                                                              >  2:279175/1‑150 (MQ=255)
      cccGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGc                                                                                                                                             >  1:670935/1‑150 (MQ=255)
         gAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGc                                                                                                                                          >  2:441182/1‑150 (MQ=255)
          aaCCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCa                                                                                                                                         <  2:1183862/150‑1 (MQ=255)
           aCCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCAt                                                                                                                                        >  2:1350758/1‑150 (MQ=255)
            ccTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATc                                                                                                                                       <  2:568708/150‑1 (MQ=255)
                ttCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGc                                                                                                                                   <  1:910974/150‑1 (MQ=255)
                 tCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCt                                                                                                                                  <  1:1124114/150‑1 (MQ=255)
                 tCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCt                                                                                                                                  <  1:1118837/150‑1 (MQ=255)
                     tGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTaaa                                                                                                                              >  2:1260468/1‑150 (MQ=255)
                      ggCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                                                             <  2:456437/150‑1 (MQ=255)
                      ggCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAAt                                                                                                                             <  1:115120/150‑1 (MQ=255)
                        cGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATcc                                                                                                                           >  1:356075/1‑150 (MQ=255)
                              cAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                                                                      >  1:1240027/1‑133 (MQ=255)
                              cAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATCACATGCCGCATcc                                                                                                                                      <  2:1240027/133‑1 (MQ=255)
                                 cTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAcc                                                                                                                  >  1:1150447/1‑150 (MQ=255)
                                  tATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTAccc                                                                                                                 >  2:1291232/1‑150 (MQ=255)
                                   aTCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCt                                                                                                                >  1:1292362/1‑150 (MQ=255)
                                    tCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCtt                                                                                                               >  2:1328487/1‑150 (MQ=255)
                                     cATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                                                              >  1:1257099/1‑150 (MQ=255)
                                     cATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTc                                                                                                              <  1:421109/150‑1 (MQ=255)
                                           aaGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaagaa                                                                                                        <  1:1260468/150‑1 (MQ=255)
                                            aGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAc                                                                                                       <  1:964280/150‑1 (MQ=255)
                                             gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                      <  2:912023/150‑1 (MQ=255)
                                             gCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAAcc                                                                                                      <  2:1291421/150‑1 (MQ=255)
                                                tGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAt                                                                                                   <  1:380190/150‑1 (MQ=255)
                                                   gTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATAcc                                                                                                >  1:1296780/1‑150 (MQ=255)
                                                       gtgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCt                                                                                            <  2:1257099/150‑1 (MQ=255)
                                                       gtgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCt                                                                                            >  1:564147/1‑150 (MQ=255)
                                                        tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTa                                                                                           <  1:368381/150‑1 (MQ=255)
                                                        tgACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTa                                                                                           >  2:883304/1‑150 (MQ=255)
                                                              tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAcc                                                                                     <  1:771746/150‑1 (MQ=255)
                                                              tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGAcc                                                                                     <  1:968549/150‑1 (MQ=255)
                                                              tGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAAAACCCGCTATTGAcc                                                                                     >  1:1139944/1‑150 (MQ=255)
                                                                    gTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGcgcc                                                                               >  1:1293941/1‑150 (MQ=255)
                                                                       tGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCcgcg                                                                            <  1:279175/150‑1 (MQ=255)
                                                                        ggAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                                                                           >  2:1053294/1‑150 (MQ=255)
                                                                          aaaaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGc                                                                         >  2:935609/1‑150 (MQ=255)
                                                                            aaaGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCAt                                                                       <  1:1343556/150‑1 (MQ=255)
                                                                                cGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgcg                                                                   <  2:1296780/150‑1 (MQ=255)
                                                                                cGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATcgcg                                                                   <  2:1150447/150‑1 (MQ=255)
                                                                                 ggTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGa                                                                  >  2:238524/1‑150 (MQ=255)
                                                                                   tGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAgg                                                                <  1:1053294/150‑1 (MQ=255)
                                                                                   ggCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAgg                                                                <  1:906929/149‑1 (MQ=255)
                                                                                    gCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGa                                                               >  1:1292871/1‑150 (MQ=255)
                                                                                       ggTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg                                                            >  1:824892/1‑150 (MQ=255)
                                                                                       ggTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATg                                                            <  2:1292871/150‑1 (MQ=255)
                                                                                           aTCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACgg                                                        <  2:1292362/150‑1 (MQ=255)
                                                                                            tCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                                                                           >  2:1215647/1‑130 (MQ=255)
                                                                                            tCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGa                                                                           <  1:1215647/130‑1 (MQ=255)
                                                                                             cGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGAt                                                      <  2:670935/150‑1 (MQ=255)
                                                                                             cGTCCGTATTGCCGAGCTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGAt                                                      >  2:650781/1‑150 (MQ=255)
                                                                                               tCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCt                                                    <  1:663452/150‑1 (MQ=255)
                                                                                                  gTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTcc                                                 >  1:1336266/1‑150 (MQ=255)
                                                                                                     ttGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgc                                              <  2:824892/150‑1 (MQ=255)
                                                                                                       gCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAgcgc                                            >  2:1209034/1‑150 (MQ=255)
                                                                                                        ccGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCa                                           >  1:993750/1‑150 (MQ=255)
                                                                                                              tAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCa                                     >  2:19003/1‑150 (MQ=255)
                                                                                                              tAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCa                                     >  1:1097540/1‑150 (MQ=255)
                                                                                                               aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaa                                    >  1:1143934/1‑150 (MQ=255)
                                                                                                               aaCGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCAGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaa                                    >  1:225122/1‑150 (MQ=255)
                                                                                                                 cGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCaaaa                                  >  1:690772/1‑150 (MQ=255)
                                                                                                                  ggCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAAt                                 >  1:1236950/1‑150 (MQ=255)
                                                                                                                        aCGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCTAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCg                           <  2:898622/150‑1 (MQ=255)
                                                                                                                             aGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaag                                                                                                          >  2:1039519/1‑66 (MQ=255)
                                                                                                                             aGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCgaag                                                                                                          <  1:1039519/66‑1 (MQ=255)
                                                                                                                             aGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGCAGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTc                      >  2:710871/1‑150 (MQ=255)
                                                                                                                             aGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTc                      <  1:1350758/150‑1 (MQ=255)
                                                                                                                                  gCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAgcgc                 <  2:1143934/150‑1 (MQ=255)
                                                                                                                                   cGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCa                <  2:1293941/150‑1 (MQ=255)
                                                                                                                                    gAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCaa               >  2:789409/1‑150 (MQ=255)
                                                                                                                                     aaaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAAt              <  1:1328487/150‑1 (MQ=255)
                                                                                                                                      aaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGCCCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATg             >  2:1231080/1‑150 (MQ=255)
                                                                                                                                      aaaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATg             >  1:1207256/1‑150 (MQ=255)
                                                                                                                                       gccggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt            >  2:1383538/4‑150 (MQ=255)
                                                                                                                                       aaGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGt            <  1:264726/150‑1 (MQ=255)
                                                                                                                                          ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGc         >  2:667902/1‑150 (MQ=255)
                                                                                                                                          ggATGCGGCATCACATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGc         >  2:1422276/1‑150 (MQ=255)
                                                                                                                                           gATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt        >  2:56028/1‑150 (MQ=255)
                                                                                                                                           gATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt        >  2:50751/1‑150 (MQ=255)
                                                                                                                                           gATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCt        >  2:217210/1‑150 (MQ=255)
                                                                                                                                             tGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTgg      >  1:911578/1‑150 (MQ=255)
                                                                                                                                                ggCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGa   >  1:335153/1‑150 (MQ=255)
                                                                                                                                                 gCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAg  <  1:1291232/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCATGAGTTTGTGACTCTGGAAGGTATGGAAAAAGCGGTGCAGGTGATCGTCCGTATTGCCGAGTTAACGGCGCAACGGAAGTAAGCGAAAAGGGATGCGGCATGTGATGCCGCATCCGGCTTAAATCCAAACTTACCCTTCGAAGAACCAATACCCGCTATTGACCAGCGCCGCGAGCATCGCGAGGAATGACGGATCTTCCAGCGCATCGCCAAAATTCTCCGCAGTCAGCGCAATGTTGCTGGCGAG  >  NC_000913/1186941‑1187235

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