Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,962,662 A→G M195T (ATG→ACG)  rhlB ← putative ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,962,6620AG99.7% 1064.8 / ‑5.5 342M195T (ATG→ACG) rhlBputative ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (161/180);  total (161/181)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 353 total aligned reads displayed.

GTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCG  >  NC_000913/3962573‑3962751
                                                                                         |                                                                                         
   ccGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGa                                                                                     <  1:3543090/93‑1 (MQ=255)
   ccGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGa                                                                                     <  1:1341243/93‑1 (MQ=255)
    cGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAg                                                                                    >  1:6533746/1‑93 (MQ=255)
     ggTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGACCAGCGTGTTGAgg                                                                                   >  1:2729094/1‑93 (MQ=255)
     ggTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAgg                                                                                   <  1:327459/93‑1 (MQ=255)
     ggTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAgg                                                                                   >  1:5869620/1‑93 (MQ=255)
     ggTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAgg                                                                                   >  2:3785361/1‑93 (MQ=255)
      gTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGgc                                                                                  >  2:7384831/1‑93 (MQ=255)
      gTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGgc                                                                                  >  2:2088660/1‑93 (MQ=255)
       ttCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGGACACGGTACGAAAGCGTGGCGGAGAACAGCGAGTTGAGgcg                                                                                 >  2:6227391/1‑93 (MQ=255)
        tCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGgcgc                                                                                <  2:4519664/93‑1 (MQ=255)
        atttcTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTTCTCGGTACGAAAGCGTGGCGGTGTATAGCGTGGTTAGgcgc                                                                                >  2:7037617/5‑93 (MQ=255)
           aCTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTgg                                                                             <  2:5737533/93‑1 (MQ=255)
             ttCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGtt                                                                           >  1:2447145/1‑93 (MQ=255)
             ttCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGtt                                                                           >  2:462642/1‑93 (MQ=255)
                aaTATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGc                                                                        >  2:5107269/1‑93 (MQ=255)
                aaTATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGc                                                                        <  1:857916/93‑1 (MQ=255)
                 atatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCa                                                                       >  2:3861509/1‑93 (MQ=255)
                 atatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCa                                                                       <  2:1885440/93‑1 (MQ=255)
                 atatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCa                                                                       <  1:1128326/93‑1 (MQ=255)
                 atatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCa                                                                       <  2:75234/93‑1 (MQ=255)
                 atatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCa                                                                       >  2:6681964/1‑93 (MQ=255)
                  tatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAg                                                                      >  2:275676/1‑93 (MQ=255)
                  tatatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAg                                                                      >  1:7390978/1‑93 (MQ=255)
                   atatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAgg                                                                     <  1:4930457/93‑1 (MQ=255)
                   atatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAgg                                                                     <  1:3458263/93‑1 (MQ=255)
                    tatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAggc                                                                    <  1:7656191/93‑1 (MQ=255)
                    tatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAggc                                                                    >  2:4722518/1‑93 (MQ=255)
                      ttCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAggcgg                                                                  >  2:1428204/1‑93 (MQ=255)
                       tCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAggcggc                                                                 <  2:6599742/93‑1 (MQ=255)
                         ggCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCAt                                                               >  1:5016260/1‑93 (MQ=255)
                           cATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGTAGGCGGCATAc                                                             <  2:7667566/93‑1 (MQ=255)
                           cATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAc                                                             <  2:4958241/93‑1 (MQ=255)
                            aTTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcg                                                            >  2:4444304/1‑93 (MQ=255)
                            aTTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcg                                                            <  2:7294842/93‑1 (MQ=255)
                             ttGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcgg                                                           >  1:7760584/1‑93 (MQ=255)
                              tGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGGAGGCGGGCTAcggg                                                          >  2:4284582/1‑92 (MQ=255)
                              tGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcggc                                                          >  1:5083397/1‑93 (MQ=255)
                               gTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcg                                                         <  2:1657801/93‑1 (MQ=255)
                               gTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcg                                                         >  2:7261774/1‑93 (MQ=255)
                                ttCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcgg                                                        >  1:5911093/1‑93 (MQ=255)
                                  cATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGaa                                                      >  1:3933914/1‑93 (MQ=255)
                                   aTCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAc                                                     <  1:570798/93‑1 (MQ=255)
                                    tCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAca                                                    <  1:7146385/93‑1 (MQ=255)
                                    tCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATAAGGCGGAAca                                                    >  2:5933925/1‑93 (MQ=255)
                                    tCTGCTCGAACGCCACTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAca                                                    >  1:2846535/1‑93 (MQ=255)
                                      tGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGCAcagc                                                  >  2:3006732/1‑93 (MQ=255)
                                      tGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagc                                                  >  2:3380748/1‑93 (MQ=255)
                                      tGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGGTTGCAGGCGGCATACGGCGGAAcagc                                                  >  1:3836785/1‑93 (MQ=255)
                                       gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                                 <  2:875267/93‑1 (MQ=255)
                                       gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                                 <  1:2548184/93‑1 (MQ=255)
                                       gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                                 >  1:2794034/1‑93 (MQ=255)
                                       gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                                 <  1:5406143/93‑1 (MQ=255)
                                       gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                                 <  1:152493/93‑1 (MQ=255)
                                       gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagcc                                                 <  1:6387594/93‑1 (MQ=255)
                                           gAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCg                                             <  1:2050608/93‑1 (MQ=255)
                                           gAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCg                                             <  1:6041274/93‑1 (MQ=255)
                                             aCGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGa                                           >  2:2218689/1‑93 (MQ=255)
                                             aCGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGa                                           <  1:1937841/93‑1 (MQ=255)
                                              cGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGat                                          <  2:4251334/93‑1 (MQ=255)
                                               gCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGata                                         <  2:6487548/93‑1 (MQ=255)
                                               gCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGata                                         <  2:3142259/93‑1 (MQ=255)
                                               gCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGata                                         >  1:4080057/1‑93 (MQ=255)
                                                 cAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATc                                       <  1:7536483/93‑1 (MQ=255)
                                                 cAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATc                                       <  1:3125513/93‑1 (MQ=255)
                                                   gTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATAtttt                                     >  2:1031133/1‑93 (MQ=255)
                                                   gTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                                     >  1:4268235/1‑93 (MQ=255)
                                                   gTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                                     >  1:4469683/1‑93 (MQ=255)
                                                    ttCACGTACCTGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCttt                                    >  1:7829679/1‑93 (MQ=255)
                                                    ttCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCttt                                    >  1:4807917/1‑93 (MQ=255)
                                                              cGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGccc                          <  1:6737861/93‑1 (MQ=255)
                                                              cGGTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGacc                          >  2:3512523/1‑93 (MQ=255)
                                                               ggTACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCa                         <  1:6966680/93‑1 (MQ=255)
                                                                 tACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGa                       >  1:552199/1‑93 (MQ=255)
                                                                 tACGAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGa                       <  2:4851240/93‑1 (MQ=255)
                                                                    gAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCg                    <  1:4555482/93‑1 (MQ=255)
                                                                    gAAAGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCg                    >  2:6872261/1‑93 (MQ=255)
                                                                      aaGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTa                  <  1:129637/93‑1 (MQ=255)
                                                                      aaGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTa                  <  2:3208323/93‑1 (MQ=255)
                                                                       aGCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTAc                 <  1:1570283/93‑1 (MQ=255)
                                                                        gCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                <  2:5282536/93‑1 (MQ=255)
                                                                        gCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                <  1:3687121/93‑1 (MQ=255)
                                                                        gCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                <  2:5754761/93‑1 (MQ=255)
                                                                        gCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                <  2:5985829/93‑1 (MQ=255)
                                                                        gCGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                <  1:6461079/93‑1 (MQ=255)
                                                                        gCGTGGCGAAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                <  2:7798260/93‑1 (MQ=255)
                                                                         cGTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACAt               <  1:2209363/93‑1 (MQ=255)
                                                                          gTGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATg              >  1:6259401/1‑93 (MQ=255)
                                                                           tGGCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGc             >  2:2708175/1‑93 (MQ=255)
                                                                            ggCGGAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCg            >  1:3738265/1‑93 (MQ=255)
                                                                               ggAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATc         <  2:1166962/93‑1 (MQ=255)
                                                                               ggAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATc         >  1:3358185/1‑93 (MQ=255)
                                                                               ggAGAACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATc         <  2:3667673/93‑1 (MQ=255)
                                                                                gagaACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCg        <  1:102021/93‑1 (MQ=255)
                                                                                 agaACAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCgg       <  2:1780042/93‑1 (MQ=255)
                                                                                   aaCAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCggtt     >  1:726080/1‑91 (MQ=255)
                                                                                   aaCAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCt     >  2:6023830/1‑93 (MQ=255)
                                                                                   aaCAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCt     >  2:5406547/1‑93 (MQ=255)
                                                                                   aaCAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCt     >  1:5616542/1‑93 (MQ=255)
                                                                                    aCAGCGTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCtt    <  1:7930387/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCG  >  NC_000913/3962573‑3962751

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