Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,178,668 A→G intergenic (+165/‑155) rplL → / → rpoB 50S ribosomal subunit protein L7/L12/RNA polymerase beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,178,6680AG98.7% 750.4 / ‑2.9 237intergenic (+165/‑155)rplL/rpoB50S ribosomal subunit protein L7/L12/RNA polymerase beta subunit
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (124/110);  total (127/110)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

TGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAAT  >  NC_000913/4178578‑4178757
                                                                                          |                                                                                         
tGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCg                                                                                         >  1:929156/1‑93 (MQ=255)
tGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCg                                                                                         >  1:1615209/1‑93 (MQ=255)
tGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCg                                                                                         >  2:2863617/1‑93 (MQ=255)
 gTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGt                                                                                        >  1:1909071/1‑93 (MQ=255)
 gTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGt                                                                                        >  2:2148179/1‑93 (MQ=255)
   aaGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTcc                                                                                      >  1:4015370/1‑93 (MQ=255)
   aaGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTcc                                                                                      <  2:4376668/93‑1 (MQ=255)
   aaGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTcc                                                                                      <  2:3911040/93‑1 (MQ=255)
   aaGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTcc                                                                                      >  2:4107329/1‑93 (MQ=255)
     ggCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGt                                                                                    <  1:468886/93‑1 (MQ=255)
     ggCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGt                                                                                    <  2:3679228/93‑1 (MQ=255)
        gcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCt                                                                                 >  1:3293130/1‑93 (MQ=255)
         ccAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtg                                                                                <  2:3459620/93‑1 (MQ=255)
         ccAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtg                                                                                >  1:3220792/1‑93 (MQ=255)
           aGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtgtg                                                                              >  1:4335896/1‑93 (MQ=255)
            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                             >  1:247567/1‑93 (MQ=255)
            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                             <  1:351073/93‑1 (MQ=255)
            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                             >  2:674935/1‑93 (MQ=255)
               gCGTTTCACGCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTaaa                                                                          >  2:741838/1‑93 (MQ=255)
               gCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTaaa                                                                          >  2:1750815/1‑93 (MQ=255)
                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATc                                                                        >  2:1293062/1‑93 (MQ=255)
                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATc                                                                        >  2:4518182/1‑93 (MQ=255)
                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCg                                                                       <  1:963980/93‑1 (MQ=255)
                   ttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGc                                                                      >  1:3968242/1‑93 (MQ=255)
                     cacaCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCaa                                                                    <  1:614015/93‑1 (MQ=255)
                       cacTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATg                                                                  >  1:4205745/1‑93 (MQ=255)
                        acTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGa                                                                 <  2:875263/93‑1 (MQ=255)
                          tGTTTGATTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGaaa                                                               <  1:4663653/93‑1 (MQ=255)
                          tGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGaaa                                                               >  1:22823/1‑93 (MQ=255)
                          tGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGaaa                                                               >  1:177792/1‑93 (MQ=255)
                           gTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATTGCAATGAAAt                                                              >  2:2685683/1‑93 (MQ=255)
                              tGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGt                                                           <  2:134577/93‑1 (MQ=255)
                                actactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGttt                                                         >  2:3934440/1‑93 (MQ=255)
                                 ctactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTa                                                        <  2:1839305/93‑1 (MQ=255)
                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTaa                                                       <  1:3735763/93‑1 (MQ=255)
                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTaa                                                       >  1:1604973/1‑93 (MQ=255)
                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTaa                                                       >  1:1388869/1‑93 (MQ=255)
                                   actGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAg                                                      <  1:3667088/93‑1 (MQ=255)
                                    ctgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGc                                                     >  2:1704489/1‑93 (MQ=255)
                                    ctgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGc                                                     <  2:2194182/93‑1 (MQ=255)
                                    ctgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGc                                                     <  2:3177896/93‑1 (MQ=255)
                                     cggtgtgCCTTTCAATGCTTGTTTCTATCGACTACTTAATATACTGCGACAGGGCGTCCTTTCTGTGTAAATCGCAATGAAATGGTTTAAGCg                                                    <  1:1277640/90‑1 (MQ=255)
                                      gctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGt                                                   >  2:3644400/1‑93 (MQ=255)
                                         gtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGAt                                                >  2:1560368/1‑93 (MQ=255)
                                          tgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATa                                               <  1:4607344/93‑1 (MQ=255)
                                          tgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATa                                               <  1:3173208/93‑1 (MQ=255)
                                          ggCCCTTCAATGCTTGTTTCGATCGACGACTTAATATACTGCGACAGGGCGTCCCGTCTGTGTAAATCGCAATTAAATGGTTTAAGCGTGATa                                               <  1:856654/92‑1 (MQ=255)
                                             cTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCa                                            >  1:3469993/1‑93 (MQ=255)
                                             cTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCa                                            >  2:4183894/1‑93 (MQ=255)
                                              tttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCaa                                           >  1:4035682/1‑93 (MQ=255)
                                              tttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCaa                                           >  2:504879/1‑93 (MQ=255)
                                                   aTGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGc                                      >  2:808177/1‑93 (MQ=255)
                                                    tGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCa                                     <  1:706610/93‑1 (MQ=255)
                                                    tGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCa                                     <  1:1948258/93‑1 (MQ=255)
                                                     gCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAt                                    >  2:593619/1‑93 (MQ=255)
                                                     gCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAt                                    <  1:3832617/93‑1 (MQ=255)
                                                      cTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAtt                                   <  1:876034/93‑1 (MQ=255)
                                                      cTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAtt                                   <  2:2268481/93‑1 (MQ=255)
                                                       ttGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTg                                  <  2:3395505/93‑1 (MQ=255)
                                                       ttGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTg                                  <  2:2556079/93‑1 (MQ=255)
                                                        tGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCTGGATAGCAACAGGCATTGc                                 <  2:3867394/93‑1 (MQ=255)
                                                         gTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCg                                <  1:2834183/93‑1 (MQ=255)
                                                         gTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCg                                >  2:1162514/1‑93 (MQ=255)
                                                         gTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCg                                >  1:1406528/1‑93 (MQ=255)
                                                          tttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCgg                               <  1:457649/93‑1 (MQ=255)
                                                          tttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCgg                               >  1:3627555/1‑93 (MQ=255)
                                                          tttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCgg                               >  2:3092845/1‑93 (MQ=255)
                                                           ttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGa                              >  2:1112218/1‑93 (MQ=255)
                                                             cTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCGATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaaa                            >  2:1437724/1‑93 (MQ=255)
                                                             cTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaaa                            <  2:1454575/93‑1 (MQ=255)
                                                              tATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAg                           >  2:1065572/1‑93 (MQ=255)
                                                              tATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAg                           >  1:2065589/1‑93 (MQ=255)
                                                               aTCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgt                          <  1:3780086/93‑1 (MQ=255)
                                                               aTCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgt                          <  2:3290645/93‑1 (MQ=255)
                                                                tCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgtg                         >  1:2389426/1‑93 (MQ=255)
                                                                tCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgtg                         >  1:3113979/1‑93 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                      >  1:2162982/1‑93 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                      <  1:4278940/93‑1 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                      <  2:2392249/93‑1 (MQ=255)
                                                                     gacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCa                    <  1:3146632/93‑1 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                   >  2:3045646/1‑93 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                   >  1:4162847/1‑93 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                   >  1:2643509/1‑93 (MQ=255)
                                                                        ttAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAttt                 <  1:4008722/93‑1 (MQ=255)
                                                                        ttAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAttt                 <  2:719160/93‑1 (MQ=255)
                                                                         tAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAtttt                <  1:1228049/93‑1 (MQ=255)
                                                                          aaTATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTc               >  2:3318971/1‑93 (MQ=255)
                                                                             ataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCgg            <  2:3496497/93‑1 (MQ=255)
                                                                               aCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTc          <  1:2782158/93‑1 (MQ=255)
                                                                                cTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTca         <  1:741838/93‑1 (MQ=255)
                                                                                cTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTca         <  1:4109606/93‑1 (MQ=255)
                                                                                 tGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaa        <  1:84597/93‑1 (MQ=255)
                                                                                 tGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaa        <  1:3915380/93‑1 (MQ=255)
                                                                                  gCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaac       <  2:1620400/93‑1 (MQ=255)
                                                                                  gCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaac       <  2:2486156/93‑1 (MQ=255)
                                                                                    gACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaacaa     >  1:1072969/1‑93 (MQ=255)
                                                                                    gACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaacaa     >  2:2605953/1‑93 (MQ=255)
                                                                                     aCAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaa    >  2:4655574/1‑93 (MQ=255)
                                                                                      cAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaaa   >  2:3124659/1‑93 (MQ=255)
                                                                                       aGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAAt  >  2:4286081/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAAT  >  NC_000913/4178578‑4178757

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: