Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,928,799 A→C 6.2% L14R (CTT→CGT)  ptrB ← oligopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,928,7990AC 6.2% 281.5 / 7.7 94L14R (CTT→CGT) ptrBoligopeptidase B
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base C (2/4);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGAT  >  NC_000913/1928650‑1928948
                                                                                                                                                     |                                                                                                                                                     
ccTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGaa                                                                                                                                                       >  2:588739/1‑150 (MQ=255)
 cTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                      >  1:861178/1‑150 (MQ=255)
 cTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                      >  2:907779/1‑150 (MQ=255)
       gATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCAt                                                                                                                                                >  2:1047062/1‑150 (MQ=255)
         tGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATgg                                                                                                                                              >  1:1460265/1‑150 (MQ=255)
          gCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGc                                                                                                                                             <  1:595390/150‑1 (MQ=255)
           cGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGGCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCg                                                                                                                                            <  2:772892/150‑1 (MQ=255)
             aTCCTGCAAGGCTTCTTGTGAGGACATCACCGGATGACCTTAACTATTTTCTTTGTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCACGCTCCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTGCATGGCgat                                                                                                                                          <  2:1233268/150‑3 (MQ=255)
              tCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgata                                                                                                                                         <  1:498316/150‑4 (MQ=255)
                 tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgatacgc                                                                                                                                      <  1:1401405/150‑7 (MQ=255)
                     aGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgc                                                                                                                                  >  1:368707/1‑150 (MQ=255)
                          tgttgtGAGGCCATCACCCGACGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcg                                                                                                                             >  2:108351/1‑150 (MQ=255)
                                 aGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTg                                                                                                                      <  2:694639/150‑1 (MQ=255)
                                 aGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTg                                                                                                                      <  1:1568751/150‑1 (MQ=255)
                                    ccATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTa                                                                                                                   >  2:491935/1‑150 (MQ=255)
                                     cATCACCCGATGACCGTAACTATTGTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCCTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAg                                                                                                                  >  2:984853/1‑150 (MQ=255)
                                         aCCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATg                                                                                                              <  2:1612949/150‑1 (MQ=255)
                                         aCCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATg                                                                                                              >  2:1447013/1‑150 (MQ=255)
                                          cccGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGt                                                                                                             <  1:1319308/150‑1 (MQ=255)
                                          cccGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGt                                                                                                             <  1:1365407/150‑1 (MQ=255)
                                            cGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTa                                                                                                           <  1:1442667/150‑1 (MQ=255)
                                             gATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTgggg                                                                                                                                       >  2:1321353/1‑121 (MQ=255)
                                             gATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTgggg                                                                                                                                       <  1:1321353/121‑1 (MQ=255)
                                               tGACCGTAACTATTTTCTTGTTGCAGGTATTCCTGGACTTCTGGCTGAGAACGAGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCTCGTATCTCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTg                                                                                                        <  2:76857/150‑1 (MQ=255)
                                                taCCGTAACTATTTTCTTGTTTCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGt                                                                                                       <  1:704010/149‑1 (MQ=255)
                                                    gTAACTATTTTCTTGTTGCAGGTAGTCAAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGttctt                                                                                                   >  1:1564341/1‑150 (MQ=255)
                                                       agtATTTTCTTGTTTCAGGTAGTCGAGGACTTCTGTCTGTGAACGCGTTTGGTCCGGGAGCGAGTAGTAATTATCGATGCGCGTATCTCCATGAAGCGTCATGTCGTGGTGAATTCGGGGGTCTTTTGGTAGCATGTTATTGttctttct                                                                                                <  2:45289/148‑1 (MQ=255)
                                                        cTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                               <  2:1457635/150‑1 (MQ=255)
                                                         tATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTgg                                                                                              <  1:1365805/150‑1 (MQ=255)
                                                            tttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTg                                                                                           <  1:911005/150‑1 (MQ=255)
                                                             ttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGa                                                                                          >  1:1092297/1‑150 (MQ=255)
                                                                 tgttgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACa                                                                                      <  2:1449867/150‑1 (MQ=255)
                                                                  gttgCAGGTAGTCCAGGACTTCTCGCTGAGAACGCGTATCGTCCCCCAGCCAGTAGTAATTATCGATGCGCGTAACGCCATGAACCGTCATGGCGTGGCGAATGCGGGCGCCATTTGGTATCATGTTATTCTTCTTTCTGGTTGAAACaa                                                                                     >  2:1417750/1‑149 (MQ=255)
                                                                    tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                   >  1:202695/1‑150 (MQ=255)
                                                                     gCAGGTAGTCCAGGACTTCTGGCAGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCtt                                                                                  >  2:959105/1‑150 (MQ=255)
                                                                      cAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTta                                                                                 <  2:368707/150‑1 (MQ=255)
                                                                        ggTAGTCAAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATTAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                               >  1:1291083/1‑150 (MQ=255)
                                                                          tAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAg                                                                             <  1:1200302/150‑1 (MQ=255)
                                                                          tAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAg                                                                             <  1:907779/150‑1 (MQ=255)
                                                                           aGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAgg                                                                            >  2:923626/1‑150 (MQ=255)
                                                                            gTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAggg                                                                           <  2:607038/150‑1 (MQ=255)
                                                                             tCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGACAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGt                                                                          >  2:215212/1‑150 (MQ=255)
                                                                              ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                         >  2:274182/1‑150 (MQ=255)
                                                                               cAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTgg                                                                        <  2:996237/150‑1 (MQ=255)
                                                                                 ggACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCa                                                                      <  1:1502295/150‑1 (MQ=255)
                                                                                    cTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCaaaa                                                                   >  1:370785/1‑150 (MQ=255)
                                                                                     ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAAc                                                                  <  1:1319271/150‑1 (MQ=255)
                                                                                     ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAAc                                                                  <  1:491935/150‑1 (MQ=255)
                                                                                      tCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACt                                                                 >  2:953211/1‑150 (MQ=255)
                                                                                      tCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACt                                                                 >  2:1556544/1‑150 (MQ=255)
                                                                                       ctgactgaGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTc                                                                >  2:697200/1‑150 (MQ=255)
                                                                                       cTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTc                                                                >  2:589505/1‑150 (MQ=255)
                                                                                       cTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTc                                                                >  1:250007/1‑150 (MQ=255)
                                                                                        tGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGTGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCa                                                               >  1:873078/1‑150 (MQ=255)
                                                                                        tGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTcc                                                               >  2:642896/1‑149 (MQ=255)
                                                                                         ggCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCAc                                                              >  2:235667/1‑150 (MQ=255)
                                                                                           cTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCg                                                            >  2:459139/1‑150 (MQ=255)
                                                                                           cTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCg                                                            >  2:676971/1‑150 (MQ=255)
                                                                                            tGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgg                                                           <  2:1353571/150‑1 (MQ=255)
                                                                                             gagaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                               <  2:910473/113‑1 (MQ=255)
                                                                                             gagaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                               >  1:910473/1‑113 (MQ=255)
                                                                                               gaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGa                                                        >  2:244725/1‑150 (MQ=255)
                                                                                                 aCGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATg                                                      <  2:1460265/150‑1 (MQ=255)
                                                                                                   gcgTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCa                                                    <  1:697200/150‑1 (MQ=255)
                                                                                                     gTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAg                                                  <  2:1564341/150‑1 (MQ=255)
                                                                                                     gTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAg                                                  >  1:1104634/1‑150 (MQ=255)
                                                                                                      tATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGTAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGc                                                 <  2:1148065/150‑1 (MQ=255)
                                                                                                        tCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGa                                               >  1:202913/1‑150 (MQ=255)
                                                                                                        tCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGa                                               <  1:1530163/150‑1 (MQ=255)
                                                                                                         cGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGaa                                              >  2:7755/1‑150 (MQ=255)
                                                                                                          gTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGaaa                                             >  2:1574090/1‑150 (MQ=255)
                                                                                                             ccGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAg                                          >  2:201579/1‑150 (MQ=255)
                                                                                                             ccGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAg                                          >  1:662299/1‑150 (MQ=255)
                                                                                                              cGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgg                                         >  2:836141/1‑150 (MQ=255)
                                                                                                                   ccAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGTTGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAAAAGGGGAGt                                    <  2:1291083/150‑1 (MQ=255)
                                                                                                                     agtagtAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCCGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGAAAGCGAAACAGGGGAGTCa                                  >  2:1537827/1‑150 (MQ=255)
                                                                                                                        agtaATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGTGGAGTCATTg                               >  1:1541225/1‑150 (MQ=255)
                                                                                                                        agtaATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTg                               <  2:1104634/150‑1 (MQ=255)
                                                                                                                         gtaATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGc                              <  2:202695/150‑1 (MQ=255)
                                                                                                                           aaTTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCtt                            >  2:51122/1‑150 (MQ=255)
                                                                                                                             ttATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAg                          <  2:202913/150‑1 (MQ=255)
                                                                                                                              tATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAga                         >  2:819826/1‑150 (MQ=255)
                                                                                                                                tCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAgatg                       <  1:108351/150‑1 (MQ=255)
                                                                                                                                   aTGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGatgatg                    <  1:929955/150‑1 (MQ=255)
                                                                                                                                      cgcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACa                 >  2:1600846/1‑150 (MQ=255)
                                                                                                                                      cgcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACa                 >  1:784961/1‑150 (MQ=255)
                                                                                                                                         gTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGt              <  2:370785/150‑1 (MQ=255)
                                                                                                                                          tATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGTGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGTAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTa             <  2:873078/150‑1 (MQ=255)
                                                                                                                                           cacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTaa            >  2:1414551/3‑150 (MQ=255)
                                                                                                                                              gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATgg         >  1:305288/1‑150 (MQ=255)
                                                                                                                                              gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATgg         <  1:1556544/150‑1 (MQ=255)
                                                                                                                                               ccATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGgc        >  1:521477/1‑150 (MQ=255)
                                                                                                                                                cATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGgcg       >  2:61361/1‑150 (MQ=255)
                                                                                                                                                 aTGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGgcgc      >  1:101613/1‑150 (MQ=255)
                                                                                                                                                   gAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCgg    <  1:1616548/150‑1 (MQ=255)
                                                                                                                                                    aaGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGa   <  1:959105/150‑1 (MQ=255)
                                                                                                                                                     aGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGat  <  1:244725/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGAT  >  NC_000913/1928650‑1928948

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