Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,185,057 C→T D710D (GAC→GAT)  rpoC → RNA polymerase beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,185,0570CT99.2% 1320.8 / ‑3.3 367D710D (GAC→GAT) rpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (194/170);  total (196/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 375 total aligned reads displayed.

CAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCT  >  NC_000913/4184967‑4185146
                                                                                          |                                                                                         
caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATgg                                                                                         <  2:4655460/93‑1 (MQ=255)
caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATgg                                                                                         <  2:5006484/93‑1 (MQ=255)
 aaaGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGt                                                                                        >  1:144280/1‑93 (MQ=255)
 aaaGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGt                                                                                        >  2:1704183/1‑93 (MQ=255)
  aaGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTc                                                                                       <  2:2250406/93‑1 (MQ=255)
  aaGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTc                                                                                       >  1:2331319/1‑93 (MQ=255)
  aaGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTc                                                                                       >  2:7726169/1‑93 (MQ=255)
   aGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTcg                                                                                      >  2:2191838/1‑92 (MQ=255)
   aGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCa                                                                                      <  1:7195342/93‑1 (MQ=255)
     ttATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAgg                                                                                    >  1:5208894/1‑93 (MQ=255)
       aTCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGaa                                                                                  >  2:5721371/1‑93 (MQ=255)
        tCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAag                                                                                 <  1:505345/93‑1 (MQ=255)
        tCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAag                                                                                 >  1:6210861/1‑93 (MQ=255)
        tCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAag                                                                                 >  1:688050/1‑93 (MQ=255)
         cGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAaga                                                                                <  2:3177568/93‑1 (MQ=255)
          gATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAagag                                                                               <  2:7370505/93‑1 (MQ=255)
           atatCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAgaga                                                                              <  1:4139210/93‑1 (MQ=255)
           atatCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAgaga                                                                              >  1:2868588/1‑93 (MQ=255)
           atatCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAgaga                                                                              <  1:3081491/93‑1 (MQ=255)
            tatCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGaa                                                                             >  2:5914340/1‑93 (MQ=255)
            tatCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGaa                                                                             >  1:4277686/1‑93 (MQ=255)
            tatCTGGGCGGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAACCCGTGATTACCCGTGATGGTCAGGAAGAGaa                                                                             <  1:5525289/93‑1 (MQ=255)
             atCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAg                                                                            <  1:6378361/93‑1 (MQ=255)
             atCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAg                                                                            >  1:6064334/1‑93 (MQ=255)
                 gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                        >  2:1914688/1‑93 (MQ=255)
                 gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                        <  2:776383/93‑1 (MQ=255)
                 gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                        >  1:1574716/1‑93 (MQ=255)
                 gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCggg                                                                        >  1:4675803/1‑93 (MQ=255)
                   gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                      <  2:898152/93‑1 (MQ=255)
                   gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                      <  1:2030951/93‑1 (MQ=255)
                   gCTGCGGAGAGCGATCGTGTAGCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTTATTAACCGTGATGGTCAGGAAGAGAAGCGGGtt                                                                      <  2:2378339/93‑1 (MQ=255)
                    cTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGttt                                                                     <  1:2182021/93‑1 (MQ=255)
                      gcggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTcc                                                                   >  2:2625328/1‑93 (MQ=255)
                       cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                  <  2:2009538/93‑1 (MQ=255)
                       cggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCt                                                                  <  1:3613917/93‑1 (MQ=255)
                        ggcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAATAGAAGCGGGTTTCCtt                                                                 >  1:4648577/1‑93 (MQ=255)
                          cgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCa                                                               <  1:4793691/93‑1 (MQ=255)
                            aaCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGAGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAAc                                                             <  1:7027217/93‑1 (MQ=255)
                             aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                            >  1:5289394/1‑93 (MQ=255)
                             aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                            >  1:4100990/1‑93 (MQ=255)
                             aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACa                                                            <  2:1590607/93‑1 (MQ=255)
                             aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCTGGTTTCCTTCAACa                                                            >  2:5550492/1‑93 (MQ=255)
                                aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                         >  1:7690362/1‑93 (MQ=255)
                                aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                         <  2:7645994/93‑1 (MQ=255)
                                aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCa                                                         <  2:6493594/93‑1 (MQ=255)
                                aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCa                                                         >  2:2946770/1‑93 (MQ=255)
                                 tCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCAt                                                        >  2:3856317/1‑93 (MQ=255)
                                 tCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCAt                                                        >  2:4297335/1‑93 (MQ=255)
                                    tgtATCCAAAGCTATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                     <  2:4726473/93‑1 (MQ=255)
                                    tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                     <  2:2052882/93‑1 (MQ=255)
                                    tgtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTa                                                     >  2:4169556/1‑93 (MQ=255)
                                     gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTAc                                                    >  2:7594114/1‑93 (MQ=255)
                                      tATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACa                                                   >  1:4597466/1‑93 (MQ=255)
                                        tCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatg                                                 >  1:6493524/1‑93 (MQ=255)
                                         ccACAGCGAGGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                <  1:2149481/93‑1 (MQ=255)
                                         ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatga                                                >  2:6261945/1‑93 (MQ=255)
                                          cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                               <  2:7104501/93‑1 (MQ=255)
                                          cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                               <  1:6689878/93‑1 (MQ=255)
                                          cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgat                                               >  1:1374438/1‑93 (MQ=255)
                                           aaaGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACatgatg                                              >  1:3792116/1‑93 (MQ=255)
                                             aGCTATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                            <  1:5697646/93‑1 (MQ=255)
                                             aGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                            <  1:25648/93‑1 (MQ=255)
                                             aGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGc                                            >  2:5091228/1‑93 (MQ=255)
                                              gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                           <  2:233331/93‑1 (MQ=255)
                                              gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                           <  2:5430508/93‑1 (MQ=255)
                                              gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGcc                                           >  2:7879535/1‑93 (MQ=255)
                                                 atgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGAc                                        <  1:4063278/93‑1 (MQ=255)
                                                  tgatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACt                                       <  2:3643382/93‑1 (MQ=255)
                                                   gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                      >  2:4605615/1‑93 (MQ=255)
                                                   gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTc                                      >  2:5827394/1‑93 (MQ=255)
                                                     tgGATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCg                                    >  2:1212139/1‑93 (MQ=255)
                                                      ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                   >  1:1070276/1‑93 (MQ=255)
                                                      ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                   <  2:5189822/93‑1 (MQ=255)
                                                      ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                   <  1:1823506/93‑1 (MQ=255)
                                                      ggATAACCTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCgg                                   <  1:3932557/93‑1 (MQ=255)
                                                             cTGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGt                            <  2:6816296/93‑1 (MQ=255)
                                                              tGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTg                           <  2:3429024/93‑1 (MQ=255)
                                                              tGCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTg                           >  2:1441645/1‑93 (MQ=255)
                                                               gCAAACTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTgg                          >  2:6929780/1‑93 (MQ=255)
                                                                  aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                       <  2:2519889/93‑1 (MQ=255)
                                                                  aaCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTc                       <  2:7344713/93‑1 (MQ=255)
                                                                   aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                      <  1:1031834/93‑1 (MQ=255)
                                                                   aCTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGGTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                      >  1:1222695/1‑93 (MQ=255)
                                                                   aCTGAAACCGTGAGTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCt                      >  1:3244989/1‑93 (MQ=255)
                                                                    cTGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATAATGGCCGACTCCGGTGCGCGTGGTTCTg                     <  1:3370839/93‑1 (MQ=255)
                                                                     tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGc                    >  1:7438942/1‑93 (MQ=255)
                                                                     tGAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGACCGCGGGGTTCTGc                    >  2:6102907/1‑93 (MQ=255)
                                                                      gAAACCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCg                   >  2:5649948/1‑93 (MQ=255)
                                                                       aaaCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCgg                  >  1:7608902/1‑93 (MQ=255)
                                                                         aCCGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGca                >  2:6418515/1‑93 (MQ=255)
                                                                          ccGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcac               >  2:6588787/1‑93 (MQ=255)
                                                                          ccGTGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGcac               >  1:2495498/1‑93 (MQ=255)
                                                                             tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa            >  1:4731560/1‑93 (MQ=255)
                                                                             tGATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGa            <  1:7267763/93‑1 (MQ=255)
                                                                              gATTAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAt           >  1:3272875/1‑93 (MQ=255)
                                                                                ttAACCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTc         <  1:5999570/93‑1 (MQ=255)
                                                                                   aCCGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTc      >  2:391877/1‑93 (MQ=255)
                                                                                    ccGTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCa     <  1:834717/93‑1 (MQ=255)
                                                                                      gTGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGc   >  1:5763269/1‑93 (MQ=255)
                                                                                       tGATGGTCAGGAAGAGAAGCGGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCt  >  1:2802413/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCT  >  NC_000913/4184967‑4185146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑