Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,486 2 bp→GC coding (316‑317/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4860CG100.0% 1327.0 / NA 410R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (190/220);  total (190/220)
*NC_0009134,001,4870GC99.8% 1321.5 / NA 409R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (188/220);  total (189/220)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

AGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAT  >  NC_000913/4001399‑4001574
                                                                                       ||                                                                                       
agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                                                     >  2:6421223/1‑93 (MQ=255)
   caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                  >  1:117914/1‑93 (MQ=255)
    aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                                                 <  2:5078375/93‑1 (MQ=255)
     ttGCATCCGTGCGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat                                                                                <  1:8220750/93‑1 (MQ=255)
      tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                               >  1:7623182/1‑93 (MQ=255)
      tGCAGCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                               <  1:2243156/93‑1 (MQ=255)
       gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                              <  2:6973508/93‑1 (MQ=255)
       gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTACACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                              <  2:3967429/93‑1 (MQ=255)
         aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                                            <  1:2357832/93‑1 (MQ=255)
          tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                                           <  2:911184/93‑1 (MQ=255)
          tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                                           >  1:5493480/1‑93 (MQ=255)
            cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa                                                                         <  2:5694280/93‑1 (MQ=255)
             ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt                                                                        <  1:8760476/93‑1 (MQ=255)
              gcggAATCGCTCGCCGACGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                                       <  2:2862808/93‑1 (MQ=255)
               cggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTg                                                                      <  1:2722359/93‑1 (MQ=255)
               cggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGAGATTg                                                                      >  1:530916/1‑93 (MQ=255)
                 gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt                                                                    <  2:9292663/93‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                 <  2:4921145/89‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                 <  1:9521308/89‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                 <  2:6684867/89‑1 (MQ=255)
                    tcgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                 >  2:8103680/5‑93 (MQ=255)
                     cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                                                <  1:796314/90‑1 (MQ=255)
                      gctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCg                                                               <  2:3332976/91‑1 (MQ=255)
                           ccGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAg                                                          >  1:639083/1‑93 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                         <  1:1481344/93‑1 (MQ=255)
                            cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                                         >  2:6409943/1‑93 (MQ=255)
                             gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                                        >  1:3697920/1‑93 (MQ=255)
                             gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                                        <  1:7819975/93‑1 (MQ=255)
                              aGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATa                                                       >  2:1340124/1‑93 (MQ=255)
                              aGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATa                                                       <  2:9718147/93‑1 (MQ=255)
                                gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                     >  1:7169991/1‑93 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                   <  2:1876681/93‑1 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                   <  2:9445032/93‑1 (MQ=255)
                                  aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                   >  2:5223626/1‑93 (MQ=255)
                                  aaaCCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                   <  1:5575048/93‑1 (MQ=255)
                                   aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                                  >  1:3608406/1‑93 (MQ=255)
                                    atcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                                                 <  1:1003374/93‑1 (MQ=255)
                                    atcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                                                 <  2:2536211/93‑1 (MQ=255)
                                     tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                                                <  2:4663265/93‑1 (MQ=255)
                                      catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                                               <  2:5247979/93‑1 (MQ=255)
                                         ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCa                                            <  2:7626329/93‑1 (MQ=255)
                                         ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCa                                            <  2:1600301/93‑1 (MQ=255)
                                         ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCa                                            >  1:3271008/1‑93 (MQ=255)
                                         ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCa                                            >  2:6021077/1‑93 (MQ=255)
                                           ccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGt                                          <  1:4180007/93‑1 (MQ=255)
                                           ccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGt                                          <  2:5829345/93‑1 (MQ=255)
                                            cAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                                         >  2:54347/1‑93 (MQ=255)
                                                acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                     >  1:5675383/1‑93 (MQ=255)
                                                acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                     <  2:8778169/93‑1 (MQ=255)
                                                 caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                                    >  2:3030750/1‑93 (MQ=255)
                                                 caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                                    >  1:9233987/1‑93 (MQ=255)
                                                 caATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGCTGccac                                    >  2:1118600/1‑93 (MQ=255)
                                                  aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                   >  1:4127871/1‑93 (MQ=255)
                                                  aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                   >  2:3757437/1‑93 (MQ=255)
                                                   aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                  >  2:8038972/1‑93 (MQ=255)
                                                    tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCaa                                 >  2:1725119/1‑93 (MQ=255)
                                                     gTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt                                >  1:6503905/1‑93 (MQ=255)
                                                      ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc                               >  1:7402471/1‑93 (MQ=255)
                                                       tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                              >  1:6114087/1‑93 (MQ=255)
                                                        cACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAg                             <  1:4995995/93‑1 (MQ=255)
                                                         aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                            >  1:939989/1‑93 (MQ=255)
                                                          ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                           <  2:6243980/93‑1 (MQ=255)
                                                          ccAGCGGGTTAATAAAGGAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa                           <  1:4402291/93‑1 (MQ=255)
                                                            aGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg                         <  1:2943086/93‑1 (MQ=255)
                                                            aGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg                         >  2:7106387/1‑93 (MQ=255)
                                                             gCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAACCAGTAGCCAGTTGCCACCAATCAGCACgg                        <  2:3260621/93‑1 (MQ=255)
                                                             gCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg                        >  1:3788392/1‑93 (MQ=255)
                                                              cGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                       <  1:6166058/93‑1 (MQ=255)
                                                              cGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                       <  1:4690196/93‑1 (MQ=255)
                                                              cGGGTTAATAAAGTAACCAAGGCGCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                       <  2:5533157/93‑1 (MQ=255)
                                                              cGGATTAATAAAGTAACCAAGGGTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                       <  1:5220292/93‑1 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                      >  1:4530003/1‑93 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                      <  2:7934336/93‑1 (MQ=255)
                                                               gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                      <  1:5852734/93‑1 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag                     >  2:2364653/1‑93 (MQ=255)
                                                                ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag                     <  1:9685551/93‑1 (MQ=255)
                                                                 gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                    <  1:6076074/93‑1 (MQ=255)
                                                                 gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                    <  2:4182797/93‑1 (MQ=255)
                                                                 gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                    <  1:1367255/93‑1 (MQ=255)
                                                                 gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                    <  2:8915294/93‑1 (MQ=255)
                                                                 gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                    <  1:484848/93‑1 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                   <  1:8513533/93‑1 (MQ=255)
                                                                  ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCagag                   >  2:8341513/1‑93 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                 <  1:7335416/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                 <  2:9651183/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                 <  2:7395004/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                 <  2:6733161/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                 <  2:1424961/93‑1 (MQ=255)
                                                                    aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGAc                 <  2:4318024/93‑1 (MQ=255)
                                                                        aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGcc             >  2:3422031/1‑93 (MQ=255)
                                                                        aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGcc             >  1:4219952/1‑93 (MQ=255)
                                                                         aGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCa            <  1:8819352/93‑1 (MQ=255)
                                                                          gTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg           >  1:2666717/1‑93 (MQ=255)
                                                                          gTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg           <  2:1821568/93‑1 (MQ=255)
                                                                          gTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg           >  1:1785742/1‑93 (MQ=255)
                                                                            aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCa         >  1:9084836/1‑93 (MQ=255)
                                                                            aaCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCa         <  1:6919537/93‑1 (MQ=255)
                                                                             aCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCaa        >  1:3498013/1‑92 (MQ=255)
                                                                                 aGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaa    <  2:1543686/93‑1 (MQ=255)
                                                                                 aGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaa    <  2:3592712/93‑1 (MQ=255)
                                                                                       ||                                                                                       
AGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAT  >  NC_000913/4001399‑4001574

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