Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,998,290 +G coding (379/486 nt) b4566 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,2891.G99.5% 571.5 / ‑4.1 182coding (380/486 nt)b4566hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (84/97);  total (85/97)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376
                                                                                             |                                                                                       
gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                          >  2:3962381/1‑93 (MQ=255)
gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                          <  2:3019748/93‑1 (MQ=255)
 aCTGCACCGCAAACTGATTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAg                                                                                        <  1:311260/93‑1 (MQ=255)
 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAg                                                                                        <  1:120936/93‑1 (MQ=255)
 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAg                                                                                        <  2:3547498/93‑1 (MQ=255)
  cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAgg                                                                                       >  2:2261890/1‑92 (MQ=255)
    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAggac                                                                                     <  1:145480/93‑4 (MQ=255)
       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGCCTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGa                                                                                   <  2:712814/93‑1 (MQ=255)
        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGAt                                                                                  <  1:3601633/93‑1 (MQ=255)
          cAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaa                                                                                <  2:631804/93‑1 (MQ=255)
           aaaCTGTATATCTTTACTGGTTTTGAACGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaaa                                                                               <  1:2283961/93‑1 (MQ=255)
           aaaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATaaa                                                                               <  1:4009176/93‑1 (MQ=255)
            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAg                                                                              >  1:3851356/1‑93 (MQ=255)
                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAtt                                                                         <  1:4090541/93‑1 (MQ=255)
                 aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAAtt                                                                         >  2:2757600/1‑93 (MQ=255)
                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGATGGATAAAGAAATTc                                                                        <  1:2635254/93‑1 (MQ=255)
                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTc                                                                        <  2:2197279/93‑1 (MQ=255)
                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTc                                                                        >  2:2054102/1‑93 (MQ=255)
                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTc                                                                        >  1:4424702/1‑93 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     >  2:891495/1‑93 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     <  1:1008998/93‑1 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     >  1:3685902/1‑93 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     <  1:3607737/93‑1 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     >  2:4667870/1‑93 (MQ=255)
                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTGTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCt                                                                     <  2:2154123/93‑1 (MQ=255)
                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCtt                                                                    <  1:318208/93‑1 (MQ=255)
                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   <  2:2575952/93‑1 (MQ=255)
                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTc                                                                   >  2:3496212/1‑93 (MQ=255)
                        ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCa                                                                  <  1:2964490/93‑1 (MQ=255)
                         tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCAc                                                                 >  2:1812522/1‑93 (MQ=255)
                          aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCAcc                                                                <  2:156954/93‑1 (MQ=255)
                          aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCAcc                                                                >  1:4341066/1‑93 (MQ=255)
                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCt                                                               >  1:52046/1‑93 (MQ=255)
                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                                              >  1:2980817/1‑93 (MQ=255)
                             ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGa                                                             <  2:4383264/93‑1 (MQ=255)
                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            <  2:3378239/93‑1 (MQ=255)
                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                                            <  1:119494/93‑1 (MQ=255)
                               ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                                           <  2:548797/93‑1 (MQ=255)
                               ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                                           >  2:13988/1‑93 (MQ=255)
                               ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                                           >  1:221343/1‑93 (MQ=255)
                                   gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgag                                                       >  1:3823759/1‑93 (MQ=255)
                                    aaGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGATGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                                      <  2:2096111/93‑1 (MQ=255)
                                     aGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAt                                                     <  2:2694147/93‑1 (MQ=255)
                                     aGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAt                                                     >  1:3495865/1‑93 (MQ=255)
                                        cATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTa                                                  <  1:1422379/93‑1 (MQ=255)
                                         aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                                                 <  1:1026056/93‑1 (MQ=255)
                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGTCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:4415150/1‑93 (MQ=255)
                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:2802970/1‑93 (MQ=255)
                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                <  2:1496087/93‑1 (MQ=255)
                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:992378/1‑93 (MQ=255)
                                            cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGc                                              <  1:2417345/93‑1 (MQ=255)
                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                             >  1:825826/1‑93 (MQ=255)
                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  2:3866408/1‑93 (MQ=255)
                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  2:1582452/1‑93 (MQ=255)
                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  1:1152396/1‑93 (MQ=255)
                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         <  1:1908116/93‑1 (MQ=255)
                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        <  1:2038876/93‑1 (MQ=255)
                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        <  2:1843474/93‑1 (MQ=255)
                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        <  1:3201564/93‑1 (MQ=255)
                                                     gaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctg                                     >  1:2003583/1‑93 (MQ=255)
                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    >  1:1469284/1‑93 (MQ=255)
                                                      aTGCCGACTTTAGCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    <  2:3887901/93‑1 (MQ=255)
                                                        gCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgctg                                  <  1:1951284/93‑1 (MQ=255)
                                                         ccGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGa                                 <  1:4399938/93‑1 (MQ=255)
                                                          cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc                                >  2:3117587/1‑93 (MQ=255)
                                                          cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc                                >  2:1668224/1‑93 (MQ=255)
                                                          cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc                                >  2:2867970/1‑93 (MQ=255)
                                                           gACTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAcc                               >  1:1732328/1‑93 (MQ=255)
                                                             cTTTATCTCCCGCGACATAAAACGCGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt                             <  2:2362538/93‑1 (MQ=255)
                                                             cTTTATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt                             >  1:899667/1‑93 (MQ=255)
                                                              tttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            >  2:2274180/1‑93 (MQ=255)
                                                              tttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            >  2:4559558/1‑93 (MQ=255)
                                                              tttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            <  2:3158794/93‑1 (MQ=255)
                                                              tttATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            >  1:3007190/1‑93 (MQ=255)
                                                                tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:2675736/93‑1 (MQ=255)
                                                                tATCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          >  2:4150561/1‑93 (MQ=255)
                                                                 aTCTCCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         >  1:4246150/1‑93 (MQ=255)
                                                                  tctcCCGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCg                        >  2:1285853/1‑93 (MQ=255)
                                                                    ccccGCGACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  2:4030186/92‑1 (MQ=255)
                                                                      ccGCGACATAAAACACGCCGCCAGGACGAATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    <  2:171055/93‑1 (MQ=255)
                                                                        gcgACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  >  2:3200212/1‑93 (MQ=255)
                                                                          gACATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                >  1:2381859/1‑93 (MQ=255)
                                                                           aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               <  1:1943488/93‑1 (MQ=255)
                                                                           aCATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  2:1136939/1‑93 (MQ=255)
                                                                            cATAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              <  2:804655/93‑1 (MQ=255)
                                                                             aTAAAACACGCCGCCAGGACGGATACAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             >  2:2989587/1‑93 (MQ=255)
                                                                             aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             >  1:1999539/1‑93 (MQ=255)
                                                                             aTAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAcc             >  1:563266/1‑93 (MQ=255)
                                                                              tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  1:4467280/1‑93 (MQ=255)
                                                                              tAAAACACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  2:51409/1‑93 (MQ=255)
                                                                               aaaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGa           <  1:4597255/93‑1 (MQ=255)
                                                                                aaaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAg          <  2:2646421/93‑1 (MQ=255)
                                                                                 aaCACGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc         <  2:3680760/93‑1 (MQ=255)
                                                                                    acGCCGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      >  2:3913892/1‑93 (MQ=255)
                                                                                     cgccgcCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggcg     >  2:3657177/1‑93 (MQ=255)
                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   <  2:1401434/93‑1 (MQ=255)
                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   <  1:3757821/93‑1 (MQ=255)
                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   >  2:3221859/1‑93 (MQ=255)
                                                                                       ccgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   <  2:2462652/93‑1 (MQ=255)
                                                                                        cgccAGGACGGATAAAGAAATTTCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt  >  2:2187290/1‑93 (MQ=255)
                                                                                             |                                                                                       
GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

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