Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,655,150 T→C 5.6% K403K (AAA→AAG pepB ← aminopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,655,1500TC 5.6% 216.1 / 2.9 72K403K (AAA→AAGpepBaminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base C (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

TTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCT  >  NC_000913/2655005‑2655298
                                                                                                                                                 |                                                                                                                                                    
ttCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                                                                                                                   <  1:1025743/149‑1 (MQ=255)
ttCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                                                                                                                   <  2:1128951/149‑1 (MQ=255)
  cGCAATTTGCTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTGCAACCGGCGCGTTACg                                                                                                                                                 <  1:821967/149‑1 (MQ=255)
       tttGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAg                                                                                                                                            <  2:487520/149‑1 (MQ=255)
          gttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCg                                                                                                                                         >  2:1577948/1‑149 (MQ=255)
             gAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCg                                                                                                                                      >  1:371548/1‑149 (MQ=255)
              aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCGTATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                                                                                                                     >  1:189024/1‑149 (MQ=255)
              aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAGCCGGCGCTTTACGGTAAGTCGCCGa                                                                                                                                     >  1:1106325/1‑149 (MQ=255)
                 ttGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCa                                                                                                                                  <  2:1245601/149‑1 (MQ=255)
                   gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                                                                                                                                >  2:1275518/1‑149 (MQ=255)
                   gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                                                                                                                                <  1:867308/149‑1 (MQ=255)
                   gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGACTTACCGTAAAGCGCCggttgaa                                                                                                                                <  1:531311/149‑7 (MQ=38)
                     aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTATGTCGCCGAGCAGTCg                                                                                                                              >  2:307627/1‑149 (MQ=255)
                      gAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGCTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGa                                                                                                                             >  1:550470/1‑149 (MQ=255)
                       aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGat                                                                                                                            <  1:1470589/149‑1 (MQ=255)
                            gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                                                                                                                       <  2:1380910/149‑1 (MQ=255)
                            gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                                                                                                                       <  1:364948/149‑1 (MQ=255)
                               tgCCGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGGCGATATGcagcc                                                                                                                    >  1:37033/2‑149 (MQ=255)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcc                                                                                                                    >  1:741629/1‑149 (MQ=255)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcc                                                                                                                    >  1:733762/1‑149 (MQ=255)
                                gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                                                                                                                   >  2:800136/1‑149 (MQ=255)
                                    gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCt                                                                                                               >  2:849021/1‑149 (MQ=255)
                                       aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc                                                                                                            <  1:839315/149‑1 (MQ=255)
                                        aGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCTCACCAAGTCCCGTAGCGCCCGCAGACGACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt                                                                                                           <  1:595325/149‑1 (MQ=255)
                                          gCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGa                                                                                                         >  2:996287/1‑149 (MQ=255)
                                            gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa                                                                                                       <  1:1546332/149‑1 (MQ=255)
                                            gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa                                                                                                       <  1:1376122/149‑1 (MQ=255)
                                            gTTTACGACGCATCAGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa                                                                                                       <  2:447303/149‑1 (MQ=255)
                                                aCGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTtc                                                                                                   <  2:427972/149‑1 (MQ=255)
                                                   ccgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTca                                                                                                <  2:1243531/149‑1 (MQ=255)
                                                    cgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGCCCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTcaa                                                                                               >  1:1302744/1‑149 (MQ=255)
                                                           ggCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTg                                                                                        >  2:652632/1‑149 (MQ=255)
                                                           ggCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTg                                                                                        >  1:424929/1‑149 (MQ=255)
                                                            gCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGc                                                                                       <  2:53971/149‑1 (MQ=255)
                                                            gCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGc                                                                                       <  2:495990/149‑1 (MQ=255)
                                                                atTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                                                                                                                                 >  2:1052692/1‑87 (MQ=255)
                                                                atTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                                                                                                                                 <  1:1052692/87‑1 (MQ=255)
                                                                 ttAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAg                                                                                  <  1:1571432/149‑1 (MQ=255)
                                                                 ttAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTGCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAg                                                                                  <  1:267943/149‑1 (MQ=255)
                                                                    gTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGaa                                                                               <  1:643745/149‑1 (MQ=255)
                                                                     tttACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAg                                                                              <  2:1210740/149‑1 (MQ=255)
                                                                           ccGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCccgcc                                                                        >  1:826291/1‑149 (MQ=255)
                                                                              ttAACAGCTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCccgccgcc                                                                     >  1:73887/1‑149 (MQ=255)
                                                                              ttAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCGTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACCGGAAGCccgccgcc                                                                     >  2:1406353/1‑149 (MQ=255)
                                                                                  cAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGc                                                                 >  1:847576/1‑149 (MQ=255)
                                                                                     aTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAGGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTcg                                                              >  1:5657/1‑149 (MQ=255)
                                                                                     aTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTcg                                                              >  2:687771/1‑149 (MQ=255)
                                                                                        aGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCgcgc                                                           <  2:937517/149‑1 (MQ=255)
                                                                                          cTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCt                                                         <  1:540268/149‑1 (MQ=255)
                                                                                          cTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCt                                                         <  1:425245/149‑1 (MQ=255)
                                                                                           tATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTg                                                        >  2:425446/1‑149 (MQ=255)
                                                                                             tCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTTAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGcc                                                      >  2:1455737/1‑149 (MQ=255)
                                                                                             tCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGcc                                                      >  2:1449263/1‑149 (MQ=255)
                                                                                                       acCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATAcgccg                                            >  1:41804/1‑149 (MQ=255)
                                                                                                       acCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATAcgccg                                            >  2:1550084/1‑149 (MQ=255)
                                                                                                             tCCCGTAGCGCCAGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTc                                      <  1:1070147/149‑1 (MQ=255)
                                                                                                              cccGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTcc                                     <  1:651707/149‑1 (MQ=255)
                                                                                                               ccGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTGGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCg                                    >  2:1288944/1‑149 (MQ=255)
                                                                                                                  tAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTa                                 <  1:652632/149‑1 (MQ=255)
                                                                                                                  tAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTa                                 <  1:1566736/149‑1 (MQ=255)
                                                                                                                        ccGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTc                           <  1:474528/149‑1 (MQ=255)
                                                                                                                        ccGCAGACCACTGTTCAACCCGCGCTTTACGGTAAGTCGCCGAGCAGTCGATCTGCAGCCAGCCTTGCTGATAGGTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTc                           <  1:1468420/149‑1 (MQ=255)
                                                                                                                         cGCAGACCACTGTTCAGCCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCa                          <  2:1106325/149‑1 (MQ=255)
                                                                                                                          gCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAg                         <  1:849021/149‑1 (MQ=255)
                                                                                                                                cACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGc                   >  1:189127/1‑149 (MQ=255)
                                                                                                                                   tGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCaaa                <  1:1577948/149‑1 (MQ=255)
                                                                                                                                   tGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCaaa                <  1:1550084/149‑1 (MQ=255)
                                                                                                                                        aaCCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAg           >  2:506064/1‑149 (MQ=255)
                                                                                                                                         ctcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGa          >  1:958803/3‑149 (MQ=255)
                                                                                                                                          tcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc         >  2:541115/2‑149 (MQ=255)
                                                                                                                                          tcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc         <  2:958803/148‑1 (MQ=255)
                                                                                                                                          ccGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc         <  2:733762/149‑1 (MQ=255)
                                                                                                                                                 tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCgg                    >  1:260148/1‑131 (MQ=255)
                                                                                                                                                 tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCgg                    <  2:260148/131‑1 (MQ=255)
                                                                                                                                                 tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGCCGGCAgct  >  1:54287/1‑149 (MQ=255)
                                                                                                                                                 tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAgct  >  2:658051/1‑149 (MQ=255)
                                                                                                                                                 tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAgct  >  1:192302/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCT  >  NC_000913/2655005‑2655298

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: