Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,878 T→A 6.0% intergenic (+20/+29) slyX → / ← slyD phi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,8780TA 6.0% 266.2 / 3.6 96intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (37/53);  new base A (3/3);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTG  >  NC_000913/3477729‑3478026
                                                                                                                                                     |                                                                                                                                                    
ggTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                                      <  2:1588396/150‑1 (MQ=255)
      cGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGATATAAGCCAGCAAGCCGTCGAACATAGCGTCGCACGCTGAAGAAACGCCACCGCCACATTATTGAGGGGTAAAAAAAGCGGGGATTCCCCg                                                                                                                                                >  2:207454/1‑150 (MQ=255)
          cATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttt                                                                                                                                            >  2:919154/1‑150 (MQ=255)
           aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttt                                                                                                                                           <  1:127677/150‑1 (MQ=255)
           aTGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGGCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtttt                                                                                                                                           <  1:1066490/150‑1 (MQ=255)
            tcaaaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACAGAGAAGTTATTAGGCAGCCAGCGGTCGAAGTTCGCGTCGCAGGCTGAATATACTCCACCGCAACATTATTTAGGCTTAAAAAAATCGGGGATTCCCCGCttttt                                                                                                                                          <  2:324688/148‑1 (MQ=255)
            tGATATGGAGTTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTTAAAGCCTGGCTGCCGTCGAACATCGCGGCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCttttt                                                                                                                                          <  2:1426757/150‑1 (MQ=255)
            tGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCttttt                                                                                                                                          <  2:1336202/150‑1 (MQ=255)
             gAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCtttttt                                                                                                                                         >  1:1247227/1‑150 (MQ=255)
              aaaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTg                                                                                                                                        <  2:603733/150‑1 (MQ=255)
               aaTGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGt                                                                                                                                       >  1:752716/1‑150 (MQ=255)
                 tGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGACCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAATCGGGGATTCCCCGCTTTTTTGTCa                                                                                                                                     <  1:346543/150‑1 (MQ=255)
                 tGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCa                                                                                                                                     <  2:785027/150‑1 (MQ=255)
                   gagaTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACt                                                                                                                                   >  1:1578954/1‑150 (MQ=255)
                    agaTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGTGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACtt                                                                                                                                  >  1:968923/1‑150 (MQ=255)
                     gaTGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCGTTTTTGTCACttt                                                                                                                                 >  1:1569911/1‑150 (MQ=255)
                        ggCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTc                                                                                                                              <  1:1198400/150‑1 (MQ=255)
                         gCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCg                                                                                                                             <  1:977683/150‑1 (MQ=255)
                          cGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTGAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCgg                                                                                                                            <  1:1533108/150‑1 (MQ=255)
                             aaCTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGttt                                                                                                                         >  1:806289/1‑148 (MQ=255)
                                tGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAg                                                                                                                      <  2:1309073/150‑1 (MQ=255)
                                  cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTg                                                                                                                    <  1:1420673/150‑1 (MQ=255)
                                  cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTg                                                                                                                    <  1:1666764/150‑1 (MQ=255)
                                  cgcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTg                                                                                                                    <  2:362210/150‑1 (MQ=255)
                                   gcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTgg                                                                                                                   <  2:1717130/150‑1 (MQ=255)
                                   gcgATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTgg                                                                                                                   <  1:472402/150‑1 (MQ=255)
                                   gcgATCATCTGCGGCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTgg                                                                                                                   <  2:1320538/150‑1 (MQ=255)
                                     gATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgca                                                                                                                 <  1:919154/150‑1 (MQ=255)
                                      atcatcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaa                                                                                                                >  2:1641664/1‑150 (MQ=255)
                                        catcTGCGGCTTCTGACCGAGAATGTAAAAGCCAGCCAGCCTTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAATAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGgcaacc                                                                                                              <  2:69024/150‑1 (MQ=255)
                                          tcTGCGTCTGCTGACCGAGAAGTTAAAAGCCAAACAGACGTCGAACATCGCGTCGCAGGCTGAAGCAACGCCACCGCCACATTATTGCGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGAGgcaaccgc                                                                                                            >  1:1224892/1‑150 (MQ=255)
                                                 ctgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgcc                                                                                                     <  2:1179262/150‑1 (MQ=255)
                                                  tgctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGGCGCAGGCTGAAGAAACTCCACCGCCAGATTATTTAGGCGTAATAAAAGCGGGGATTCCGCGCTTTTTTGTCACTTTTTCGTTATTAGTGGCAACCGCATccgccg                                                                                                    <  1:1312213/150‑1 (MQ=255)
                                                  tgcagACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGGCTCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAAccgccg                                                                                                    <  1:123671/150‑1 (MQ=255)
                                                   gctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGt                                                                                                   >  1:55890/1‑150 (MQ=255)
                                                    ctgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGtt                                                                                                  >  2:1437663/1‑150 (MQ=255)
                                                     tgACCGAGAAGTTAACAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTg                                                                                                 <  1:676254/150‑1 (MQ=255)
                                                     tgACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTg                                                                                                 <  1:1633511/150‑1 (MQ=255)
                                                     tgACCGAGAAGTTAAAAGCCAGCCAGACGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACGGCCACATTATTGAGGCGTAAATAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTg                                                                                                 <  1:271283/150‑1 (MQ=255)
                                                      gACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGc                                                                                                >  2:552549/1‑150 (MQ=255)
                                                       aCCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGcc                                                                                               <  1:1057638/150‑1 (MQ=255)
                                                        ccGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCt                                                                                              <  2:112695/150‑1 (MQ=255)
                                                             aaGTTAAAAGCCAGCCAGCCGTCGAACCTCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAcc                                                                                         <  1:1332790/150‑1 (MQ=255)
                                                             aaGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTAcc                                                                                         <  2:754205/150‑1 (MQ=255)
                                                                ttAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTTAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCTCCGTTGCCTTTACCGcc                                                                                      <  2:1224892/150‑1 (MQ=255)
                                                                 tAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCGTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCca                                                                                     <  1:206292/150‑1 (MQ=255)
                                                                  aaaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcac                                                                                    <  2:995138/150‑1 (MQ=255)
                                                                   aaaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCcaca                                                                                   <  1:166454/150‑1 (MQ=255)
                                                                    aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcag                                                                                  >  2:1385321/1‑150 (MQ=255)
                                                                    aaGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCAcag                                                                                  <  2:278307/150‑1 (MQ=255)
                                                                      gccagccagccGTCGAACATCGCGTCACAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGTCAACCGCCACCGCCGTTGCCTTTACCGCCAcagca                                                                                >  1:958946/1‑150 (MQ=255)
                                                                        cagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc                                                                              <  2:520859/150‑1 (MQ=255)
                                                                        cagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc                                                                              <  2:996912/150‑1 (MQ=255)
                                                                        cagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACagcagc                                                                              >  1:1292199/1‑150 (MQ=255)
                                                                               ccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGcc                                                                       <  1:738770/150‑1 (MQ=255)
                                                                                 gTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCAc                                                                     >  2:207972/1‑150 (MQ=255)
                                                                                   cGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCGCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCg                                                                   >  1:867449/1‑150 (MQ=255)
                                                                                    gAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                  <  2:30648/150‑1 (MQ=255)
                                                                                    gAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                  <  1:635006/150‑1 (MQ=255)
                                                                                    gAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTATTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgt                                                                  >  2:1654408/1‑150 (MQ=255)
                                                                                      aCATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCgtgt                                                                >  1:1571002/1‑150 (MQ=255)
                                                                                         tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCCCAGCAGCCTTCGCCACCGTGTTCa                                                             >  1:41016/1‑150 (MQ=255)
                                                                                         tCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCa                                                             <  1:792221/150‑1 (MQ=255)
                                                                                             gTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGAc                                                         >  1:936433/1‑150 (MQ=255)
                                                                                                 cAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTg                                                     <  2:7608/150‑1 (MQ=255)
                                                                                                   ggCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGAt                                                   <  2:1247227/150‑1 (MQ=255)
                                                                                                   ggCTGAAGAAACGCCACCGCCACATTAATGAGGCGTAAAAAAAGAGGGGATTCACCGCTATTTTGACACTTTTTCGGTATTAGTGGCAACCACAACCGCCGTTGCCTTTACCGCCCCAGCAGCCTTCGCCACCGTGTTCATGACCGTGAt                                                   >  2:736532/1‑150 (MQ=255)
                                                                                                    gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATc                                                  >  1:567130/1‑150 (MQ=255)
                                                                                                    gCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATc                                                  <  2:914940/150‑1 (MQ=255)
                                                                                                     cTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCg                                                 >  1:861408/1‑150 (MQ=255)
                                                                                                       gaagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTGTATTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGATGCCTTTACCGCCACAGCAGCCTTCCCCACCGTGTTCAAGACCGAGATCGTg                                               >  1:203088/1‑150 (MQ=255)
                                                                                                        aagaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTgg                                              <  2:752716/150‑1 (MQ=255)
                                                                                                         agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCGACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                             >  1:1236889/1‑150 (MQ=255)
                                                                                                         agaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACGGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGc                                             <  1:1379300/150‑1 (MQ=255)
                                                                                                          gaaACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGcc                                            <  2:1007397/150‑1 (MQ=255)
                                                                                                            aaCGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCCCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCaa                                          >  2:1230227/1‑149 (MQ=255)
                                                                                                              cgccaccgccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGa                                        >  1:1094436/1‑150 (MQ=255)
                                                                                                                 cacagccacATTATTGAGGCGTAAAACAAGCGGTGATACCCCGCTTTTTTGTCACTTTTTCGGTATTATTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGGCCACGATCACGGTCATGaaca                                     <  1:154633/150‑4 (MQ=255)
                                                                                                                   ccgccacATTATTGAGGCGTAAAAAAAGCGGGGAATCCCAGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGGCCACGATCACGGTCATGaacacg                                   <  1:954444/150‑6 (MQ=255)
                                                                                                                     gccacATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCATACGAAACCGCCTTTGCGTTTACATCCACAGCAGACTTCGCCACCGTGTTAATGAACGTGATCGTGGCCATGACggccgc                                 >  1:309211/1‑150 (MQ=255)
                                                                                                                         caTTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGTCAACCGCAACCGCCGTTGCCTTTACCGCCAAAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCgcagca                             >  1:805261/1‑150 (MQ=255)
                                                                                                                            tattGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTGTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCGTTACCGCCCGAGCCGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCCTGACCGCCGCAGCAGcc                          >  1:62797/1‑150 (MQ=255)
                                                                                                                                    tgacaaaaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt                  >  2:995304/5‑150 (MQ=255)
                                                                                                                                    cGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCCACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCCACCgtcggggg                  >  1:603521/1‑149 (MQ=255)
                                                                                                                                        aaaaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgtg              >  1:803339/1‑150 (MQ=255)
                                                                                                                                        aaaaaaGCGGGGAATCCACGCTTTTATGTCACTTTTTCGGTATTAGTGGAAACCGCAACCGCCGTTGCCTTTAACTCAACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggtcgtg              >  1:1620128/1‑150 (MQ=255)
                                                                                                                                          aaaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGAt            <  2:1716093/150‑1 (MQ=255)
                                                                                                                                          aaaaGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGAt            >  2:1467827/1‑150 (MQ=255)
                                                                                                                                           aaaGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCACCCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATc           >  1:1249662/1‑150 (MQ=255)
                                                                                                                                            aaGCGGGGATTCGGCGCTTTTTTGTCACTTTTTCGGTATTCGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGGAGCCTTCGCCTGCGTGTTGATGAGCGTGTTAGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCg          <  2:922337/150‑1 (MQ=255)
                                                                                                                                             aGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCTCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACAGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGGAACCGTCGTGGTGATGATCgg         >  2:160901/1‑149 (MQ=255)
                                                                                                                                               cGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtgg       <  1:839165/150‑1 (MQ=255)
                                                                                                                                               cGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtgg       >  2:1238869/1‑150 (MQ=255)
                                                                                                                                               cGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtgg       <  1:592123/150‑1 (MQ=255)
                                                                                                                                                ggggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggt      <  1:1476245/150‑1 (MQ=255)
                                                                                                                                                  ggtttCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTTGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtgg    <  2:1578954/147‑1 (MQ=255)
                                                                                                                                                    aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggtg  <  2:1571002/150‑1 (MQ=255)
                                                                                                                                                    aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCgtggtggtg  <  2:41016/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTGACCGCTCATGAAATGGAGATGGCGAAACTGCGCGATCATCTGCGTCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGTCGTGATCGTGGTGGTG  >  NC_000913/3477729‑3478026

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: