Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,810,386 C→T T261T (ACG→ACA b1730 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,810,3860CT84.1% 584.6 / 71.9 208T261T (ACG→ACAb1730hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base T (90/85);  total (106/102)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

AAACCGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTATC  >  NC_000913/1810294‑1810475
                                                                                            |                                                                                         
aaaCCGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGcc                                                                                           >  2:3358346/1‑93 (MQ=255)
  aCCGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGt                                                                                         >  2:1766356/1‑93 (MQ=255)
  aCCGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGt                                                                                         <  2:2745965/93‑1 (MQ=255)
   ccGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGtt                                                                                        <  2:5781751/93‑1 (MQ=255)
   ccGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGtt                                                                                        >  2:134475/1‑93 (MQ=255)
    cGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTg                                                                                       >  2:5476893/1‑93 (MQ=255)
     gATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGc                                                                                      >  1:4253257/1‑93 (MQ=255)
     gATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGc                                                                                      >  1:5561568/1‑93 (MQ=255)
     gATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGAGGCTGTTGc                                                                                      >  1:6602650/1‑93 (MQ=255)
       tAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCt                                                                                    <  1:6158672/93‑1 (MQ=255)
       tAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCt                                                                                    >  1:4419416/1‑93 (MQ=255)
         aaGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTa                                                                                  >  1:1910158/1‑93 (MQ=255)
          aGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTAc                                                                                 >  1:3101734/1‑93 (MQ=255)
          aGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTAc                                                                                 >  1:3680222/1‑93 (MQ=255)
           gCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACt                                                                                >  1:2866647/1‑93 (MQ=255)
                tCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTa                                                                           <  1:4887989/93‑1 (MQ=255)
                   ggATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAAt                                                                        <  2:1652775/93‑1 (MQ=255)
                   ggATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAAt                                                                        <  2:6373309/93‑1 (MQ=255)
                   ggATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAAt                                                                        <  1:4554855/93‑1 (MQ=255)
                    gATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATg                                                                       <  1:410489/93‑1 (MQ=255)
                     aTAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGt                                                                      >  2:4268433/1‑93 (MQ=255)
                       aGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTc                                                                    <  1:5545951/93‑1 (MQ=255)
                        gCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTcc                                                                   >  1:3764517/1‑93 (MQ=255)
                        gCGAAGACGACGCCGAACTGGCACACCATTAATCCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTcc                                                                   >  2:4301517/1‑93 (MQ=255)
                         cGAAGACGACGCCGAACTGGCACACCATTAATGCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCa                                                                  <  1:2641111/93‑1 (MQ=255)
                          gAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCAt                                                                 >  1:5511303/1‑93 (MQ=255)
                          gAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCAt                                                                 >  1:1596063/1‑93 (MQ=255)
                          gAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGATGTTGCCTTATTTTTTAAATGTTCCAt                                                                 >  2:4139395/1‑93 (MQ=255)
                           aaGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCTGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATc                                                                >  2:6435569/1‑93 (MQ=255)
                           aaGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATc                                                                >  1:677802/1‑93 (MQ=255)
                           aaGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATc                                                                >  2:4993999/1‑93 (MQ=255)
                           aaGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATc                                                                <  1:1378188/93‑1 (MQ=255)
                           aaGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATc                                                                >  2:4157024/1‑93 (MQ=255)
                             gacgacGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATcgc                                                              >  2:1291418/1‑93 (MQ=255)
                               cgacgCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCgcgg                                                            <  2:1414052/93‑1 (MQ=255)
                                 acgCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCgcggcg                                                          >  1:3880875/1‑93 (MQ=255)
                                  cgCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGa                                                         <  2:2548584/93‑1 (MQ=255)
                                    ccGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTCCTTTTTAAATGTTCCATCGCGGCGATg                                                       >  1:363160/1‑93 (MQ=255)
                                    ccGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATg                                                       <  2:399813/93‑1 (MQ=255)
                                    ccGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATg                                                       >  2:4488055/1‑93 (MQ=255)
                                    ccGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATg                                                       >  2:812365/1‑93 (MQ=255)
                                    ccGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATg                                                       >  1:1712544/1‑93 (MQ=255)
                                    ccGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATg                                                       >  1:3527669/1‑93 (MQ=255)
                                       aaCTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCAt                                                    <  2:5361810/93‑1 (MQ=255)
                                        aCTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATa                                                   >  2:6240497/1‑93 (MQ=255)
                                         cTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATAc                                                  >  1:2280103/1‑93 (MQ=255)
                                         cTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATAc                                                  <  1:4195216/93‑1 (MQ=255)
                                              acacCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACt                                             >  1:261243/1‑93 (MQ=255)
                                              acacCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACt                                             >  1:3033419/1‑93 (MQ=255)
                                               cacCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTc                                            >  1:6075627/1‑93 (MQ=255)
                                                  cATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCata                                         <  2:723422/93‑1 (MQ=255)
                                                  cATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCata                                         >  2:527737/1‑93 (MQ=255)
                                                  cATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCata                                         >  2:2351663/1‑93 (MQ=255)
                                                   aTTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCatat                                        >  2:2038229/1‑93 (MQ=255)
                                                      aaTTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTa                                     <  2:4647519/93‑1 (MQ=255)
                                                       aTTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTaa                                    <  2:4948197/93‑1 (MQ=255)
                                                       aTTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTaa                                    >  2:6608133/1‑93 (MQ=255)
                                                       aTTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTaa                                    >  2:4898519/1‑93 (MQ=255)
                                                        ttCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTaaa                                   >  1:5274021/1‑93 (MQ=255)
                                                          cTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACa                                 >  2:100593/1‑93 (MQ=255)
                                                          cTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACa                                 <  2:2982328/93‑1 (MQ=255)
                                                            tATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGa                               <  1:1041592/93‑1 (MQ=255)
                                                            tATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGa                               <  2:3669809/93‑1 (MQ=255)
                                                             aTGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGAt                              <  2:4525134/93‑1 (MQ=255)
                                                             aTGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGAt                              <  1:2513019/93‑1 (MQ=255)
                                                             aTGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGAt                              <  1:2577670/93‑1 (MQ=255)
                                                              tGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATc                             >  2:3616777/1‑93 (MQ=255)
                                                              tGCTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATc                             >  1:1987142/1‑93 (MQ=255)
                                                              tGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATc                             <  2:2369434/93‑1 (MQ=255)
                                                                cTGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGt                           <  2:5607924/93‑1 (MQ=255)
                                                                 tGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTc                          >  2:3501028/1‑93 (MQ=255)
                                                                 tGGCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTc                          >  1:2902471/1‑93 (MQ=255)
                                                                 tGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTc                          >  2:4180736/1‑93 (MQ=255)
                                                                   gCAACGGTGCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAg                        <  1:1875940/93‑1 (MQ=255)
                                                                   gCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAg                        <  1:2392750/93‑1 (MQ=255)
                                                                   gCAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAg                        >  2:6552466/1‑93 (MQ=255)
                                                                    cAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGc                       >  2:4751001/1‑93 (MQ=255)
                                                                    cAACGGTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGc                       <  2:3097709/93‑1 (MQ=255)
                                                                      aCGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGa                     >  2:5494798/1‑93 (MQ=255)
                                                                        ggTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGaaa                   <  1:871551/93‑1 (MQ=255)
                                                                         gTCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGaaaa                  <  2:4864866/93‑1 (MQ=255)
                                                                          tCCGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGaaaaa                 >  2:5190243/1‑93 (MQ=255)
                                                                           ccGTTTTGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGaaaaaa                >  1:2538208/1‑93 (MQ=255)
                                                                              ttttGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATg             <  2:4545774/93‑1 (MQ=255)
                                                                                ttGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTc           >  1:2128292/1‑93 (MQ=255)
                                                                                 tGTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTcc          <  2:505394/93‑1 (MQ=255)
                                                                                  gTATAGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCa         <  1:1018958/93‑1 (MQ=255)
                                                                                     taGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACt      <  2:3961049/93‑1 (MQ=255)
                                                                                     taGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACt      >  2:517399/1‑93 (MQ=255)
                                                                                      aGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGTTCGTCAGCGAAAAAAATGTCCAACtt     <  1:987241/93‑1 (MQ=255)
                                                                                      aGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACtt     >  2:1545652/1‑93 (MQ=255)
                                                                                      aGGGGCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACtt     >  2:734843/1‑93 (MQ=255)
                                                                                       ggggCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTa    <  1:3254817/93‑1 (MQ=255)
                                                                                       ggggCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTa    >  2:2024133/1‑93 (MQ=255)
                                                                                        gggCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTAt   <  2:1059760/93‑1 (MQ=255)
                                                                                        gggCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACGCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTAt   <  2:3565706/93‑1 (MQ=255)
                                                                                        gggCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTAt   >  2:5747739/1‑93 (MQ=255)
                                                                                         ggCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTATc  <  2:6041487/93‑1 (MQ=255)
                                                                                         ggCTGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTATc  <  1:851786/93‑1 (MQ=255)
                                                                                         ggCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTATc  <  1:5822997/93‑1 (MQ=255)
                                                                                            |                                                                                         
AAACCGATAAAGCCATTCTGGATAGCGAAGACGACGCCGAACTGGCACACCATTAATTCTTATGCTGGCAACGGTCCGTTTTGTATAGGGGCCGTTGCCTTACTTTTTAAATGTTCCATCGCGGCGATGCATACTGACTCATATTTAAACAGATCGTCAGCGAAAAAAATGTCCAACTTATC  >  NC_000913/1810294‑1810475

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