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,9930AG99.8% 1274.6 / ‑5.5 398A82A (GCT→GCCyieGputative membrane/transport protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (185/212);  total (186/212)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 409 total aligned reads displayed.

GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAG  >  NC_000913/3893903‑3894082
                                                                                          |                                                                                         
gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         >  2:6742783/1‑93 (MQ=255)
gACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGc                                                                                         >  2:5235113/1‑93 (MQ=255)
 aCCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCa                                                                                        <  2:7524113/93‑1 (MQ=255)
  ccTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCaa                                                                                       >  2:1029789/1‑93 (MQ=255)
    tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGGAAAc                                                                                     >  2:5762389/1‑93 (MQ=255)
    tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCaagc                                                                                     <  1:1154829/93‑3 (MQ=255)
    tGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAAc                                                                                     >  1:5355309/1‑93 (MQ=255)
     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                    <  2:3821960/93‑1 (MQ=255)
     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                    <  1:6343925/93‑1 (MQ=255)
     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                    <  1:7317343/93‑1 (MQ=255)
     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                    <  2:4502917/93‑1 (MQ=255)
     gCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACa                                                                                    <  1:1269087/93‑1 (MQ=255)
         cGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTcc                                                                                <  2:2828781/93‑1 (MQ=255)
           gCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCca                                                                              >  2:250194/1‑93 (MQ=255)
              aGCCCATCGCCTGTACAACGACAAAAGCGAAGAAAGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatca                                                                           >  1:3820218/1‑93 (MQ=255)
               gCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCcatcat                                                                          >  2:421633/1‑93 (MQ=255)
                 ccATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataa                                                                        >  1:7269565/1‑93 (MQ=255)
                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                       >  2:848036/1‑93 (MQ=255)
                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                       >  2:2216639/1‑93 (MQ=255)
                  cATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataat                                                                       >  2:5965813/1‑93 (MQ=255)
                   aTCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCataata                                                                      <  2:3956735/93‑1 (MQ=255)
                    tCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATAc                                                                     >  1:610691/1‑93 (MQ=255)
                     cGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACt                                                                    >  2:5834472/1‑93 (MQ=255)
                      gCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTg                                                                   >  1:7935325/1‑93 (MQ=255)
                       ccTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGc                                                                  <  1:649977/93‑1 (MQ=255)
                        cTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGcc                                                                 >  2:5132641/1‑93 (MQ=255)
                         tGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCg                                                                >  1:976686/1‑93 (MQ=255)
                          gTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGa                                                               >  2:3572949/1‑93 (MQ=255)
                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGaa                                                              >  2:7200591/1‑93 (MQ=255)
                           tACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGaa                                                              <  1:3583769/93‑1 (MQ=255)
                            aCAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAAt                                                             >  1:5216265/1‑93 (MQ=255)
                             cAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATg                                                            >  1:5960017/1‑93 (MQ=255)
                              caCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                           <  1:2564435/92‑1 (MQ=255)
                              aaCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                           <  1:4153001/93‑1 (MQ=255)
                              aaCGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGc                                                           >  1:4309784/1‑93 (MQ=255)
                                 gACAAAAGCGAAGAACGCATTCAAGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGc                                                        <  2:5394805/93‑1 (MQ=255)
                                  aCAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCg                                                       >  2:5665727/1‑93 (MQ=255)
                                   cAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGa                                                      >  2:20202/1‑93 (MQ=255)
                                    aaaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGAt                                                     <  2:730804/93‑1 (MQ=255)
                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                    >  1:6372450/1‑93 (MQ=255)
                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                    <  2:7390844/93‑1 (MQ=255)
                                     aaaGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATc                                                    <  2:3332035/93‑1 (MQ=255)
                                       aGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAg                                                  <  1:2920954/93‑1 (MQ=255)
                                        gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGa                                                 >  2:6934050/1‑93 (MQ=255)
                                        gCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGa                                                 <  1:13569/93‑1 (MQ=255)
                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              <  2:5011041/93‑1 (MQ=255)
                                           aagaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAg                                              <  1:7169343/93‑1 (MQ=255)
                                            agaaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAgg                                             <  1:1215531/93‑1 (MQ=255)
                                              aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                           <  1:908415/93‑1 (MQ=255)
                                              aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                           >  1:2360903/1‑93 (MQ=255)
                                              aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                           >  2:6592099/1‑93 (MQ=255)
                                              aaCGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTa                                           <  2:1665100/93‑1 (MQ=255)
                                               aCGCATTCAGGCCAATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATCCTGCCGAATGCAGCGATCAGACAGGTAg                                          <  2:7150700/93‑1 (MQ=255)
                                                cGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAgg                                         >  2:6819959/1‑92 (MQ=255)
                                                cGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGt                                         <  1:6293893/93‑1 (MQ=255)
                                                  cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGa                                       >  1:2097192/1‑93 (MQ=255)
                                                  cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGa                                       >  2:4873011/1‑93 (MQ=255)
                                                  cATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGa                                       >  1:2836613/1‑93 (MQ=255)
                                                   aTTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGAc                                      <  1:4434178/93‑1 (MQ=255)
                                                    ttCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACg                                     >  2:2848523/1‑93 (MQ=255)
                                                     tCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGa                                    <  2:4632550/93‑1 (MQ=255)
                                                      cAGGGCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGaa                                   <  2:5511080/93‑1 (MQ=255)
                                                       aGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAg                                  >  1:3360130/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 <  2:2288228/93‑1 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  1:7793938/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  2:7713483/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 <  1:4643982/93‑1 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  1:6688013/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 >  1:3249366/1‑93 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 <  2:1501259/93‑1 (MQ=255)
                                                        ggCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGa                                 <  2:1894205/93‑1 (MQ=255)
                                                          cccATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACg                               <  1:2209933/93‑1 (MQ=255)
                                                           ccATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgt                              >  2:3631550/1‑92 (MQ=255)
                                                           ccATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgt                              <  2:7858543/93‑2 (MQ=255)
                                                           ccATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgt                              <  1:3202514/93‑2 (MQ=255)
                                                            cATAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACgtc                             <  1:1821420/93‑3 (MQ=255)
                                                              tAGCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGtcgc                           >  1:1438742/1‑93 (MQ=255)
                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         <  1:3613102/93‑1 (MQ=255)
                                                                gCGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTg                         <  1:3086496/93‑1 (MQ=255)
                                                                 cGGGTGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTgg                        <  1:1536750/93‑1 (MQ=255)
                                                                     tGCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATc                    <  2:1333105/93‑1 (MQ=255)
                                                                      gCCAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATcc                   <  2:468857/93‑1 (MQ=255)
                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  >  1:5636701/1‑93 (MQ=255)
                                                                       ccAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCa                  <  1:5487632/93‑1 (MQ=255)
                                                                        cAGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCAt                 <  1:4805579/93‑1 (MQ=255)
                                                                         aGTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATg                >  2:3133429/1‑93 (MQ=255)
                                                                          gTGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGc               >  2:3453612/1‑93 (MQ=255)
                                                                           tGCAACTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGcc              <  1:4938204/93‑1 (MQ=255)
                                                                              aaCTGGCAGGTTGGCAAACCGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc           <  1:7695420/93‑1 (MQ=255)
                                                                              aaCTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGc           >  1:2442568/1‑93 (MQ=255)
                                                                                cTGGCCGGTTGGCAAGAAGTCCCATCATAATCCTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa         <  1:6624630/93‑1 (MQ=255)
                                                                                cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa         <  2:6091285/93‑1 (MQ=255)
                                                                                cTGGCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCaa         >  2:6416771/1‑93 (MQ=255)
                                                                                cTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCaa         <  1:7458649/93‑1 (MQ=255)
                                                                                 tGGCAGGTAGGCAAACAGTCCCATCATAAAACTGCCGAATGCAGAGATCAGACAGGAAGAGACGAAGACGTCGCTGGTATCCATGACAGCAAc        <  1:6994578/93‑1 (MQ=255)
                                                                                  gacAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGAAACg       <  2:4095081/91‑1 (MQ=255)
                                                                                   gCAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGc      <  2:1254303/93‑1 (MQ=255)
                                                                                    cAGGTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGcc     >  2:6241089/1‑93 (MQ=255)
                                                                                      ggTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCaa   >  1:6869893/1‑93 (MQ=255)
                                                                                       gTTGGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGTCGCTGGTATCCATGCCAGCAACGCCAAg  >  2:1422554/1‑93 (MQ=255)
                                                                                          |                                                                                         
GACCTGCCACGGCAAGCCCATCGCCTGTACAACGACAAAAGCGAAGAACGCATTCAGGCCCATAGCGGGTGCCAGTGCAACTGGCAGGTTAGCAAACAGTCCCATCATAATACTGCCGAATGCAGCGATCAGACAGGTAGTGACGAAGACGGCGCTGGTATCCATGCCAGCAACGCCAAG  >  NC_000913/3893903‑3894082

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