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 ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,928,7990AC 6.2% 270.1 / 6.6 95L14R (CTT→CGT) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base A (42/46);  new base C (2/4);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCG  >  NC_000913/1928651‑1928945
                                                                                                                                                    |                                                                                                                                                  
cTTTAAGATGCGATCCTGGAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                   <  1:2024096/150‑1 (MQ=255)
cTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                   >  2:1060511/1‑150 (MQ=255)
cTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAg                                                                                                                                                   <  1:190331/150‑1 (MQ=255)
 tttAAGATGCTATCCTGCACGGCGTGTTGTGAGGGCATCACCCGATGACCGTAAGTATTTTCTTGTTTCAGTTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGc                                                                                                                                                  <  1:128785/150‑1 (MQ=255)
    aaGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTTTCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTc                                                                                                                                               >  1:457613/1‑150 (MQ=255)
       aTGCGATCCTGCAAGGCTTGTTGTGCGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGGAGGCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATg                                                                                                                                            <  2:106032/150‑1 (MQ=255)
          ggATCCTGCAAGGCTTGTTGTGAGGCCATCACGCGATGAAAGTATCTTTTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGAAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCg                                                                                                                                         <  1:473266/149‑1 (MQ=255)
           gATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTAGCGTCCCGCAGCCAGTAGTAATTTTCGATGCGCGTATCGCCATGAAGCGTCATGGCGt                                                                                                                                        >  2:1961464/1‑150 (MQ=255)
            aTCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgat                                                                                                                                       <  1:573732/150‑3 (MQ=255)
            aTCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgat                                                                                                                                       <  2:527669/150‑3 (MQ=255)
             cccTGCAAGGCTGGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCGGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgata                                                                                                                                      <  2:1044354/149‑4 (MQ=255)
                tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCTTCATGGCgatacgc                                                                                                                                   <  1:32190/150‑7 (MQ=255)
                tGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGa                                                                                                                                   >  2:891252/1‑150 (MQ=255)
                tGCAAGGCATGTTGTGAGGCCATCACCCGATGACAGTAACAATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCAGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGCTGAGCGTATCGCCAAGAAGCGTCATGGCGTgggta                                                                                                                                   >  2:1336630/1‑148 (MQ=255)
                 gCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGaa                                                                                                                                  >  2:303648/1‑150 (MQ=255)
                       cTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAATCGTCATGGAGTGGGGAATgcggg                                                                                                                            >  2:59074/1‑150 (MQ=255)
                         tgttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcg                                                                                                                          >  1:272655/1‑150 (MQ=255)
                         tgttgtGAGGCCATCACCCGATGACCGTAACTATTTTCTCGTTGCAGGCAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATgcgggcg                                                                                                                          <  1:142616/150‑1 (MQ=255)
                            tgtgAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCt                                                                                                                       >  1:1910802/1‑150 (MQ=255)
                               gAGGCCGTCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGAAGTAATTATCGATGCGCGTATCGCCATGCAGCGTCATGGCGTGTGGAATGCGGGCGGCtttt                                                                                                                    >  1:188588/1‑150 (MQ=255)
                               gAGGCCATCACCTGATGACCGTAACTATTTTCTTTTTGCAGGTAGTCCAGGACTTCTTGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGCTGCTCGTATCGCCATGAAGCGTCATGGCGTTGGGAATGCGGGCGGCtttt                                                                                                                    <  1:328875/150‑1 (MQ=255)
                               gAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCtttt                                                                                                                    >  2:1371136/1‑150 (MQ=255)
                                  gCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTAGCGTCCCGCAGCCAGTAGTAATTATCGATGCGCATATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGt                                                                                                                 <  2:1535363/150‑1 (MQ=255)
                                   ccATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAAAGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTa                                                                                                                >  2:1194365/1‑150 (MQ=255)
                                    cATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAg                                                                                                               >  1:1072277/1‑150 (MQ=255)
                                    cATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAg                                                                                                               >  2:1142525/1‑150 (MQ=255)
                                     aTCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGc                                                                                                              >  2:266807/1‑150 (MQ=255)
                                     aTCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGc                                                                                                              <  1:1234321/150‑1 (MQ=255)
                                      tCACCCGATGACCGTAACTATTTTCTTGTCGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCa                                                                                                             >  1:605459/1‑150 (MQ=255)
                                       cACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTCGCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCAt                                                                                                            <  1:1861422/150‑1 (MQ=255)
                                         cccGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGt                                                                                                          <  2:1314262/150‑1 (MQ=255)
                                         cccGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGt                                                                                                          >  2:559258/1‑150 (MQ=255)
                                          ccGATGACCGTAAATATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTCTCGGCCCGCAGCCAGTAGAAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGtt                                                                                                         <  2:120056/150‑1 (MQ=255)
                                           cGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTa                                                                                                        >  1:1230104/1‑150 (MQ=255)
                                            gATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTAt                                                                                                       <  1:303648/150‑1 (MQ=255)
                                                aCCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGtt                                                                                                   >  1:1426137/1‑150 (MQ=255)
                                                       cTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                            >  2:1080241/1‑150 (MQ=255)
                                                       cTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTg                                                                                            >  1:1281457/1‑150 (MQ=255)
                                                       cTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCAAGGCGTGGCGAATGCGGGCGGCTTTTGGTAGCATGTTATTGGTCTTTCTg                                                                                            >  2:1404433/1‑150 (MQ=255)
                                                          ttttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGtt                                                                                         >  1:435807/1‑150 (MQ=255)
                                                          ttttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATGGTTCTTTCTGGtt                                                                                         >  2:1514885/1‑150 (MQ=255)
                                                          ttttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTGTTGGTAGCAGGTTATTGTTGTTTCTGGtt                                                                                         >  2:1928880/1‑150 (MQ=255)
                                                           tttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGGGTCATGGCGTGTGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTg                                                                                        <  1:1514995/150‑1 (MQ=255)
                                                           tttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTATGGTTg                                                                                        >  2:52748/1‑150 (MQ=255)
                                                            ttCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGa                                                                                       <  2:2013184/150‑1 (MQ=255)
                                                              cTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaaa                                                                                     >  2:1253211/1‑150 (MQ=255)
                                                              cTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCACGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGACGCGTCAGGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaaa                                                                                     >  2:333564/1‑150 (MQ=255)
                                                              cTTGCTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGGCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGaaa                                                                                     <  2:350985/150‑1 (MQ=255)
                                                                   tgCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGGGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCt                                                                                <  1:1563845/150‑1 (MQ=255)
                                                                       ggTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTtata                                                                            <  2:689158/150‑1 (MQ=255)
                                                                         tAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAg                                                                          >  2:777228/1‑150 (MQ=255)
                                                                         tAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAg                                                                          >  2:769801/1‑150 (MQ=255)
                                                                          aGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAgg                                                                         >  1:947717/1‑150 (MQ=255)
                                                                            tCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGt                                                                       >  1:710723/1‑150 (MQ=255)
                                                                             ccAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTg                                                                      <  2:1644140/150‑1 (MQ=255)
                                                                                    ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCTTGGCGTGGGTAATTCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAAc                                                               <  2:2100794/150‑1 (MQ=255)
                                                                                    ttCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAAc                                                               <  1:1553456/150‑1 (MQ=255)
                                                                                        ggCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCAc                                                           >  1:962075/1‑150 (MQ=255)
                                                                                             agaACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                      <  2:1850453/150‑1 (MQ=255)
                                                                                             agaACGCGTATCGTCCCGCAGCCAGTAGTAAATATCGATGCGCGTATCGCCATGAAGGGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCgggg                                                      <  1:1978784/150‑1 (MQ=255)
                                                                                              caaCGCGTACCGTCCCGCAGCCAGTAGTAAATATCTATGCGCGTATCGCCATGAAGCGTCATGTCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGa                                                     <  2:1031648/149‑1 (MQ=255)
                                                                                                   cgTTTCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGATGCGTCATGGCGTGGGGAATGCGTGCGTCTTTTGGTAGCATGTTATTGTTCTTTCTGGGTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCaa                                                <  1:326705/150‑1 (MQ=255)
                                                                                                   cgTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCaa                                                >  1:603610/1‑150 (MQ=255)
                                                                                                    gTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAg                                               <  1:1860800/150‑1 (MQ=255)
                                                                                                     tATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGc                                              <  1:1142525/150‑1 (MQ=255)
                                                                                                        cGTCCCGCAGCCAGTTGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTTTTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACAGGGGATGCAAGCGaa                                           >  2:1193768/1‑150 (MQ=255)
                                                                                                         gTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGaaa                                          <  2:1281457/150‑1 (MQ=255)
                                                                                                          tCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAAc                                         >  1:1602062/1‑150 (MQ=255)
                                                                                                            ccGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAg                                       >  1:1464111/1‑150 (MQ=255)
                                                                                                            ccGCAGCCAGTAGTAATAATCGATGCGCGTAACGCCATGAAGCGTCATGGCGTGTGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCTGGGATGCTAGCGAAACAg                                       <  1:1714684/150‑1 (MQ=255)
                                                                                                               cAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAgggg                                    <  2:443106/150‑1 (MQ=255)
                                                                                                                 gCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGTGGTATGCTGGTGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAATCTTCTTATAAGGGTGGAAAAACTCACCGGGGATGCAAGCGAAACAGGGGAg                                  <  1:1336630/150‑1 (MQ=255)
                                                                                                                 gCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAg                                  >  2:405953/1‑150 (MQ=255)
                                                                                                                 gCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAg                                  <  2:1195446/150‑1 (MQ=255)
                                                                                                                  ccAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGt                                 <  1:1490901/150‑1 (MQ=255)
                                                                                                                  ccAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGt                                 <  2:1932323/150‑1 (MQ=255)
                                                                                                                   cAGTAGTAATTATCGATGCTCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTc                                <  2:89543/150‑1 (MQ=255)
                                                                                                                   cAGTAGTAATTATCGATGCGAGAATCGCCATGAAGCGTCTTGGCGGGGGGAATGCGGGGGGCTTTTGGTAGCATGTTATTGTTCGTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTc                                <  2:1595586/150‑1 (MQ=255)
                                                                                                                           aTTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTa                        >  1:1792463/1‑150 (MQ=255)
                                                                                                                             tataGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACTTCTTATACGGGTGGCAAAACTCACCGGGGATGAAAGCGAAACAGGGGAGTCATTGCTTAga                      >  2:1566677/1‑150 (MQ=255)
                                                                                                                             tATCGATGCGCGTATCGCCATGAAGCGGCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAga                      <  2:659275/150‑1 (MQ=255)
                                                                                                                               tCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCTGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAgatg                    <  1:686080/150‑1 (MQ=255)
                                                                                                                               tCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAgatg                    >  1:1459974/1‑150 (MQ=255)
                                                                                                                               tCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAgatg                    <  1:857970/150‑1 (MQ=255)
                                                                                                                                 gATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGATATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAgatgat                  <  2:1940969/150‑1 (MQ=255)
                                                                                                                                  ttGCGCGTATCGCCATGAAGCGGCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGatgatg                 <  2:347919/149‑1 (MQ=255)
                                                                                                                                   tGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGAtgatga                <  1:1650278/150‑1 (MQ=255)
                                                                                                                                      gcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAg             <  2:1874859/150‑1 (MQ=255)
                                                                                                                                      gcgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAg             <  1:1823686/150‑1 (MQ=255)
                                                                                                                                       cgTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAgg            <  2:947717/150‑1 (MQ=255)
                                                                                                                                        gTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAAAAGGGGAGTCATTGCTTAGATGATGACAGGt           >  1:1087470/1‑150 (MQ=255)
                                                                                                                                         ccacGCCATGACGCTTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTa          >  2:396646/4‑150 (MQ=255)
                                                                                                                                             gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATgg      <  1:1397558/150‑1 (MQ=255)
                                                                                                                                             gCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTAAAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATgg      <  1:1257467/150‑1 (MQ=255)
                                                                                                                                                 tGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGcgcg  >  2:454461/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTATCGCCATGAAGCGTCATGGCGTGGGGAATGCGGGCGGCTTTTGGTAGCATGTTATTGTTCTTTCTGGTTGAAACATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCG  >  NC_000913/1928651‑1928945

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