Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,893,993 A→G A82A (GCT→GCC yieG ← putative membrane/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,893,9930AG96.7% 562.0 / 1.7 183A82A (GCT→GCCyieGputative membrane/transport protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (97/78);  total (101/82)
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.99e-01

Only 100 of 188 total aligned reads displayed.

GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAG  >  NC_000913/3893903‑3894082
                                                                                          |                                                                                         
gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACTCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         >  1:1741112/1‑93 (MQ=255)
gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         <  2:2149244/93‑1 (MQ=255)
gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         >  2:2811748/1‑93 (MQ=255)
gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         <  1:1828832/93‑1 (MQ=255)
   cTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCaaa                                                                                      >  1:346554/1‑93 (MQ=255)
    tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAAc                                                                                     >  2:159588/1‑93 (MQ=255)
      ccACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAg                                                                                   >  2:1814136/1‑93 (MQ=255)
         cGGCAAGCCGATCGCCTGTACAACGACAAAAGCGAAGAGCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAGCAGTcc                                                                                <  1:709372/93‑1 (MQ=255)
          ggCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTccc                                                                               <  2:3267031/93‑1 (MQ=255)
           gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCca                                                                              >  1:2081970/1‑93 (MQ=255)
             aaGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACCGTCCcatc                                                                            >  1:3427820/1‑93 (MQ=255)
             aaGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatc                                                                            >  2:4485933/1‑93 (MQ=255)
             aaGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatc                                                                            >  2:2705744/1‑93 (MQ=255)
              aGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatca                                                                           >  1:517771/1‑93 (MQ=255)
               gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatcat                                                                          >  2:3982384/1‑93 (MQ=255)
                 ccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataa                                                                        <  2:341370/93‑1 (MQ=255)
                 ccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataa                                                                        >  2:1437166/1‑93 (MQ=255)
                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                       >  1:3051729/1‑93 (MQ=255)
                   aTCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataata                                                                      >  2:1456927/1‑93 (MQ=255)
                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                     <  1:4175476/93‑1 (MQ=255)
                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                     >  1:2663705/1‑93 (MQ=255)
                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                     <  1:929551/93‑1 (MQ=255)
                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                     >  2:4088965/1‑93 (MQ=255)
                     cGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACt                                                                    <  1:1203625/93‑1 (MQ=255)
                      gCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTg                                                                   >  1:3457664/1‑93 (MQ=255)
                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                  <  2:3817562/93‑1 (MQ=255)
                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                  >  1:4501772/1‑93 (MQ=255)
                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                  >  2:3038898/1‑93 (MQ=255)
                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                  >  2:2255612/1‑93 (MQ=255)
                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGACCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                  >  1:2550347/1‑93 (MQ=255)
                         tGTACAACGACAAAAGCGAAGAACGCATTGAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCg                                                                >  2:1354340/1‑93 (MQ=255)
                          gTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCCACTGGCAGGCTGGCAAACAGTCCCATCATAATACTGCCGa                                                               >  1:2311922/1‑93 (MQ=255)
                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGaa                                                              >  2:2448483/1‑93 (MQ=255)
                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                             <  2:48797/93‑1 (MQ=255)
                             cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGAGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAAtc                                                            >  2:4246158/1‑92 (MQ=255)
                              aaCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                           <  1:1984277/93‑1 (MQ=255)
                              aaCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                           >  2:1113656/1‑93 (MQ=255)
                               aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCa                                                          >  1:3979643/1‑93 (MQ=255)
                               aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCa                                                          >  1:1528810/1‑93 (MQ=255)
                               aCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCa                                                          <  2:1940368/93‑1 (MQ=255)
                                cGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAg                                                         >  1:4305203/1‑93 (MQ=255)
                                 gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCNAACAGTCCCATCATAATACTGCCGAAtgctgc                                                        >  2:2421577/1‑93 (MQ=255)
                                 gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGc                                                        <  2:4368652/93‑1 (MQ=255)
                                 gACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGc                                                        <  2:3295732/93‑1 (MQ=255)
                                    aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGAt                                                     >  2:3444731/1‑93 (MQ=255)
                                    aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGAt                                                     >  1:2222858/1‑93 (MQ=255)
                                    aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGAt                                                     >  1:2936594/1‑93 (MQ=255)
                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                    <  1:589828/93‑1 (MQ=255)
                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                    >  1:2333410/1‑93 (MQ=255)
                                      aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGGAACTGGCCGGTTGGGAAACAGTCCCATCATCATACTGCCGAATGAAGCGATCa                                                   >  2:191339/1‑93 (MQ=255)
                                      aaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCa                                                   <  2:2747980/93‑1 (MQ=255)
                                        gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGa                                                 >  2:1664751/1‑93 (MQ=255)
                                        gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGa                                                 >  1:4527680/1‑93 (MQ=255)
                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              >  1:192919/1‑93 (MQ=255)
                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              <  1:361870/93‑1 (MQ=255)
                                            agaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAgg                                             >  1:3249942/1‑93 (MQ=255)
                                            agaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGACATCAGACAgg                                             >  1:833393/1‑93 (MQ=255)
                                              aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                           <  1:421142/93‑1 (MQ=255)
                                                 gCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                        >  2:1166839/1‑93 (MQ=255)
                                                 gCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTg                                        <  1:1322528/93‑1 (MQ=255)
                                                  cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGa                                       >  1:3338396/1‑93 (MQ=255)
                                                   aTTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGAc                                      >  1:109589/1‑93 (MQ=255)
                                                       aGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAg                                  <  2:608835/93‑1 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  1:1136897/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  1:4612185/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 <  1:3014905/93‑1 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  2:1873085/1‑93 (MQ=255)
                                                         gCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGAc                                <  1:2170927/93‑1 (MQ=255)
                                                         gCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGAc                                >  1:801961/1‑93 (MQ=255)
                                                           ccATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgt                              <  1:479342/93‑2 (MQ=255)
                                                             aTAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcg                            <  1:175477/93‑1 (MQ=255)
                                                              tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                           <  1:3789467/93‑1 (MQ=255)
                                                              tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                           <  2:3086045/93‑1 (MQ=255)
                                                              tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                           <  2:4383027/93‑1 (MQ=255)
                                                               aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                          >  2:2984878/1‑93 (MQ=255)
                                                               aGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCt                          >  1:1416518/1‑93 (MQ=255)
                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         <  2:373433/93‑1 (MQ=255)
                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         >  2:3547351/1‑93 (MQ=255)
                                                                  gggTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGt                       >  1:4131892/1‑93 (MQ=255)
                                                                      gCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATcc                   <  1:3764539/93‑1 (MQ=255)
                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  >  2:927780/1‑93 (MQ=255)
                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  <  2:4716002/93‑1 (MQ=255)
                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  <  2:2673526/93‑1 (MQ=255)
                                                                          gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTCTCCATGc               >  2:2450999/1‑93 (MQ=255)
                                                                          gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGc               >  1:596085/1‑93 (MQ=255)
                                                                          gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGc               <  2:4536052/93‑1 (MQ=255)
                                                                           tGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGcc              <  2:1569629/93‑1 (MQ=255)
                                                                            gCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCa             <  2:3682854/93‑1 (MQ=255)
                                                                             cAACTGGCAGTTTGGCAAACAGTCCCATAATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAg            <  2:2492528/93‑1 (MQ=255)
                                                                              aaCTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc           <  2:4585439/93‑1 (MQ=255)
                                                                                cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa         <  2:2899734/93‑1 (MQ=255)
                                                                                cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa         <  1:1708780/93‑1 (MQ=255)
                                                                                 tGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAAc        >  2:2598967/1‑93 (MQ=255)
                                                                                 tGGCAGGTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACg       >  2:2640341/1‑93 (MQ=255)
                                                                                  ggCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGACGCTGGTATCCATGCCAGCAACg       >  2:1289793/1‑93 (MQ=255)
                                                                                    cAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGcc     >  2:1997816/1‑93 (MQ=255)
                                                                                    cAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGcc     <  2:2460951/93‑1 (MQ=255)
                                                                                      ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa   <  1:1810959/93‑1 (MQ=255)
                                                                                      ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa   <  2:3961808/93‑1 (MQ=255)
                                                                                       gTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAg  >  2:3188948/1‑93 (MQ=255)
                                                                                          |                                                                                         
GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAG  >  NC_000913/3893903‑3894082

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