Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,304,479 C→G 5.6% P216A (CCG→GCG)  oppC → oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,304,4790CG 5.6% 689.3 / 25.0 198P216A (CCG→GCG) oppColigopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (78/107);  new base G (7/4);  total (86/112)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 203 total aligned reads displayed.

GATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGAT  >  NC_000913/1304331‑1304627
                                                                                                                                                    |                                                                                                                                                    
    ttCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGaa                                                                                                                                                  <  1:1004990/149‑1 (MQ=255)
     tCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAAt                                                                                                                                                 <  3:1302732/149‑1 (MQ=255)
        tGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGTAGCGCAAAGTGCCTATTGTGGCGGCACGAGTTGGCGGTGTATCCACATCGGGCGTTGTTATTCGCCAGAATGTGTCGAATGTa                                                                                                                                              >  2:763527/1‑149 (MQ=255)
         ggCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTAc                                                                                                                                             <  2:88966/149‑1 (MQ=255)
         ggCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTAc                                                                                                                                             >  4:637853/1‑149 (MQ=255)
              tAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGgt                                                                                                                                        <  3:637853/149‑1 (MQ=255)
                ggCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGgtgt                                                                                                                                      <  3:1121969/149‑1 (MQ=255)
                ggCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGgtgt                                                                                                                                      <  2:1111679/149‑1 (MQ=255)
                 gCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGtgtg                                                                                                                                     >  2:869371/1‑149 (MQ=255)
                 gCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGtgtg                                                                                                                                     >  2:688550/1‑149 (MQ=255)
                 gCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGtgtg                                                                                                                                     <  4:342628/149‑1 (MQ=255)
                   aTGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggt                                                                                                                                   >  2:1181042/1‑149 (MQ=255)
                    tGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGgtg                                                                                                                                       >  2:1462944/1‑144 (MQ=255)
                    tGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtg                                                                                                                                  <  3:46566/149‑1 (MQ=255)
                     ggTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtgg                                                                                                                                 >  1:74719/1‑149 (MQ=255)
                             ggctggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTAc                                                                                                                         <  1:688550/149‑1 (MQ=255)
                                tggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCa                                                                                                                      >  4:715300/1‑149 (MQ=255)
                                  gATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATc                                                                                                                    <  3:867303/149‑1 (MQ=255)
                                  gATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATc                                                                                                                    <  1:89998/149‑1 (MQ=255)
                                    tatGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCAc                                                                                                                  <  3:564158/149‑1 (MQ=255)
                                     atGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACt                                                                                                                 <  4:1180624/149‑1 (MQ=255)
                                        gCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCt                                                                                                              >  4:575878/1‑149 (MQ=255)
                                         cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTg                                                                                                             <  4:50830/149‑1 (MQ=255)
                                         cTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTg                                                                                                             <  1:444460/149‑1 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCt                                                                                                              >  1:129779/1‑147 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                            >  3:1229807/1‑149 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                            >  4:156528/1‑149 (MQ=255)
                                          tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                            >  1:421194/1‑149 (MQ=255)
                                           cGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGt                                                                                                           <  3:1395844/149‑1 (MQ=255)
                                               ttGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGccc                                                                                                       >  3:33939/1‑149 (MQ=255)
                                                tgtgCGTGGGTAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCa                                                                                                      <  1:637772/149‑1 (MQ=255)
                                                      tGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGa                                                                                                <  3:156528/149‑1 (MQ=255)
                                                          cAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCct                                                                                            <  4:1364634/149‑1 (MQ=255)
                                                           aaaCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCctc                                                                                           >  4:811804/1‑149 (MQ=255)
                                                             aCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                                         <  3:1189479/149‑1 (MQ=255)
                                                              cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCTACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCTTCttt                                                                                        <  1:456845/149‑1 (MQ=255)
                                                                   aGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATc                                                                                   <  1:955736/149‑1 (MQ=255)
                                                                    gTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                                         >  4:1486327/1‑142 (MQ=255)
                                                                    gTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                                         <  3:1486327/142‑1 (MQ=255)
                                                                       tGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTc                                                                               <  1:194201/149‑1 (MQ=255)
                                                                        gAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTcc                                                                              >  3:907194/1‑149 (MQ=255)
                                                                        gAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTcc                                                                              >  3:636506/1‑149 (MQ=255)
                                                                        gAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTcc                                                                              >  1:1464284/1‑149 (MQ=255)
                                                                          aGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCtt                                                                            <  2:1567733/149‑1 (MQ=255)
                                                                          aGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCtt                                                                            <  4:439306/149‑1 (MQ=255)
                                                                           gcgcAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTa                                                                           <  3:715300/149‑1 (MQ=255)
                                                                           gcgcAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTa                                                                           <  1:858156/149‑1 (MQ=255)
                                                                            cgcAAAGAGTTTATTGTGGCGGCACAAGTTGGCGGTGTATCGACATAGGGCATTGATATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTAAGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAg                                                                          <  4:937084/149‑1 (MQ=255)
                                                                             gcAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGc                                                                         <  4:128597/149‑1 (MQ=255)
                                                                                aaGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTc                                                                      >  2:742227/1‑149 (MQ=255)
                                                                                   agTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTg                                                                   <  4:329794/149‑1 (MQ=255)
                                                                                    gTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTgg                                                                  >  1:338788/1‑149 (MQ=255)
                                                                                    gTTTATTGGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTggg                                                                 <  2:87155/149‑1 (MQ=255)
                                                                                        aTTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGtt                                                              <  2:1193281/149‑1 (MQ=255)
                                                                                        aTTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGtt                                                              <  2:370798/149‑1 (MQ=255)
                                                                                        aTTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGtt                                                              >  2:489651/1‑149 (MQ=255)
                                                                                          tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                            <  3:811804/149‑1 (MQ=255)
                                                                                          tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                            <  4:1229807/149‑1 (MQ=255)
                                                                                              gcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATc                                                                                              >  3:522813/1‑111 (MQ=255)
                                                                                              gcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATc                                                                                              <  4:522813/111‑1 (MQ=255)
                                                                                                  cacaAGTTGGCGGTGTATCGACATCGGACATTGTTATTCGCCACAATGTGCCGAAGGTACTCGTTGTGGTGGAGGTCTAAGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCATCCTGGCGTTGGGTACGCaa                                                    >  4:762767/1‑149 (MQ=255)
                                                                                                    caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                                  <  3:969588/149‑1 (MQ=255)
                                                                                                          ggCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGtt                                            >  2:614106/1‑149 (MQ=255)
                                                                                                          ggCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGtt                                            <  2:1309955/149‑1 (MQ=255)
                                                                                                            cGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTaa                                          <  3:108068/149‑1 (MQ=255)
                                                                                                              gtgtATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACACAAGAGCCGTTAagc                                        >  3:1158518/1‑149 (MQ=255)
                                                                                                               tgtATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAagca                                       <  2:437961/149‑1 (MQ=255)
                                                                                                                 tATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAagcagc                                     >  1:1580150/1‑149 (MQ=255)
                                                                                                                   tCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTg                                   >  1:237031/1‑149 (MQ=255)
                                                                                                                     gACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTggg                                 >  1:79003/1‑149 (MQ=255)
                                                                                                                         tCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATc                                                                                              >  4:286097/1‑84 (MQ=255)
                                                                                                                         tCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGc                             <  3:950997/149‑1 (MQ=255)
                                                                                                                            ggCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCAtt                          >  1:530749/1‑149 (MQ=255)
                                                                                                                              cgagtacaTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGc                        >  2:1530869/8‑149 (MQ=255)
                                                                                                                               aTTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTcc                                                                              >  4:1105602/1‑94 (MQ=255)
                                                                                                                               aTTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTaa                                          >  4:1049141/1‑130 (MQ=255)
                                                                                                                               aTTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTaa                                          <  3:1049141/130‑1 (MQ=255)
                                                                                                                                 tGTTATTCGCTACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCGAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa                     <  4:1344121/149‑1 (MQ=255)
                                                                                                                                 tGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa                     <  4:94704/149‑1 (MQ=255)
                                                                                                                                 tGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGa                     <  2:713129/149‑1 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCt                                                                    >  3:419953/1‑100 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCt                       <  3:453757/145‑1 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCt                       >  4:453757/1‑145 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt                   <  1:1181042/149‑1 (MQ=255)
                                                                                                                                   ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt                   <  4:907194/149‑1 (MQ=255)
                                                                                                                                    tATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTg                  <  1:120925/149‑1 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATc                                                                                   <  4:1246182/83‑1 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATc                                                                                   >  3:1246182/1‑83 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                               >  3:942382/1‑135 (MQ=255)
                                                                                                                                     aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGa                 >  1:567751/1‑149 (MQ=255)
                                                                                                                                      ttCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGAt                <  2:203904/149‑1 (MQ=255)
                                                                                                                                          ccACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGgcg            >  4:1412394/1‑149 (MQ=255)
                                                                                                                                           cacaTTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGgcgc           >  1:1509832/1‑149 (MQ=255)
                                                                                                                                                tgtgGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACt      <  2:567751/149‑1 (MQ=255)
                                                                                                                                                 gtgGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTc     <  1:1530869/149‑1 (MQ=255)
                                                                                                                                                   gCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGa   <  1:1401649/149‑1 (MQ=255)
                                                                                                                                                   gCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGa   <  1:1381788/149‑1 (MQ=255)
                                                                                                                                                   gCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGa   >  3:1544945/1‑149 (MQ=255)
                                                                                                                                                    ccGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTgggg                                >  1:110481/1‑119 (MQ=255)
                                                                                                                                                    ccGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGAt  <  1:614106/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGCGCGAACTCGAT  >  NC_000913/1304331‑1304627

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