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,4860CG99.8% 1711.5 / NA 436R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (208/227);  total (209/227)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,001,4870GC99.8% 1700.1 / ‑6.3 436R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (209/226);  minor base A (0/1);  total (209/227)
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.78e-01

Only 100 of 451 total aligned reads displayed.

GCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTT  >  NC_000913/4001396‑4001576
                                                                                          ||                                                                                         
 ccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                                                      <  2:9789900/96‑1 (MQ=255)
 ccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                                                      <  1:7737660/96‑1 (MQ=255)
 ccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                                                      <  2:4464951/96‑1 (MQ=255)
   agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCCGc                                                                                    >  1:9738080/1‑96 (MQ=255)
   agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                    <  2:294029/96‑1 (MQ=255)
   agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                                                                                    <  2:3115404/96‑1 (MQ=255)
    gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                                                   >  1:724757/1‑96 (MQ=255)
      caTTGCATCCGGCGGAATGGCTCGCCGAGGAAAATCATCACCAGCACAATGTTCACCAGCGAGTTAATAAAGTAACCAAGGCTCGCTTCCACCata                                                                                 >  1:12643251/1‑96 (MQ=255)
      caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                                 >  1:10875053/1‑96 (MQ=255)
      caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAAACATCCCCCGCACAATGTTCACCAGCGGGCTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                                                 >  2:5996753/1‑96 (MQ=255)
       aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                                >  1:1349339/1‑96 (MQ=255)
        ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg                                                                               <  1:5324363/96‑1 (MQ=255)
        ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg                                                                               <  2:11880525/96‑1 (MQ=255)
          gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTGCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                                             >  2:10296436/1‑96 (MQ=255)
             tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt                                                                          <  1:7917492/96‑1 (MQ=255)
              ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTCCCAGCCTATGGTGAtt                                                                         >  1:1539766/1‑96 (MQ=255)
                ggcggACTCGCTCGCGGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTACCCAAGGCGCGCTTCCAGCATATGGTGATTGt                                                                       <  1:5591291/96‑1 (MQ=255)
                ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt                                                                       <  2:5195566/96‑1 (MQ=255)
                  cggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                                     >  2:10941323/1‑96 (MQ=255)
                  cggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                                     <  1:12525161/96‑1 (MQ=255)
                   ggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCa                                                                    <  2:10528673/96‑1 (MQ=255)
                    gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCCCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                   <  1:8037502/96‑1 (MQ=255)
                    gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                                                   >  2:11048558/1‑96 (MQ=255)
                        cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                               <  2:13344710/93‑1 (MQ=255)
                        cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGcc                                                               <  1:5910181/93‑1 (MQ=255)
                         gctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                                              <  1:9988603/94‑1 (MQ=255)
                         gctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                                              <  2:8645531/94‑1 (MQ=255)
                          ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                                             <  1:14284882/95‑1 (MQ=255)
                          ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                                             >  2:10671081/2‑96 (MQ=255)
                          ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                                             >  2:5815387/2‑96 (MQ=255)
                           gcGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAg                                                            <  1:3628329/95‑1 (MQ=255)
                             gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                                          >  2:1261282/1‑96 (MQ=255)
                             gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                                          >  1:12244063/1‑96 (MQ=255)
                              ccGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                                             <  2:6694336/92‑1 (MQ=255)
                               cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                        >  1:13964211/1‑96 (MQ=255)
                               cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                        >  1:8893230/1‑96 (MQ=255)
                                gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                       >  1:11200099/1‑96 (MQ=255)
                                gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                                       >  1:3877955/1‑96 (MQ=255)
                                  ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACa                                                     >  2:7299116/1‑96 (MQ=255)
                                   gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                                    <  1:3969178/96‑1 (MQ=255)
                                   gAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                                                    >  1:4247985/1‑96 (MQ=255)
                                    aaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                                                   >  2:5545246/1‑96 (MQ=255)
                                    aaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                                                   >  1:8447665/1‑96 (MQ=255)
                                    aaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                                                   <  2:6794443/96‑1 (MQ=255)
                                     aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGTCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                                                  <  1:3495857/96‑1 (MQ=255)
                                     aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                                                  >  1:9158920/1‑96 (MQ=255)
                                      aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                                                 <  1:2809632/96‑1 (MQ=255)
                                        tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                                               >  1:1798957/1‑96 (MQ=255)
                                          ctcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                                             >  1:11434360/2‑96 (MQ=255)
                                          atcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                                             <  2:4812629/96‑1 (MQ=255)
                                          atcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                                             <  2:1660898/96‑1 (MQ=255)
                                            ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                                           >  2:8229912/1‑96 (MQ=255)
                                            ccccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                                           >  1:4903640/1‑96 (MQ=255)
                                              ccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                                         <  2:4205306/96‑1 (MQ=255)
                                              ccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                                         >  2:64226/1‑96 (MQ=255)
                                               cAGCACAATGTTCACCAGCGGGTTAATAAAGTACCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcc                                        <  2:381900/96‑1 (MQ=255)
                                                aGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                                       <  2:1146665/96‑1 (MQ=255)
                                                 gCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                                      >  2:1923775/1‑96 (MQ=255)
                                                 gCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccac                                      <  2:13187300/96‑1 (MQ=255)
                                                  cacaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                     <  2:10804534/96‑1 (MQ=255)
                                                  cacaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                                     <  2:5023445/96‑1 (MQ=255)
                                                  cacaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCACTTGccacc                                     <  2:12956638/96‑1 (MQ=255)
                                                   tcaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca                                    <  1:2195556/95‑1 (MQ=255)
                                                     aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt                                  >  2:1400275/1‑96 (MQ=255)
                                                     aaTCTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt                                  >  1:2464380/1‑96 (MQ=255)
                                                      aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc                                 <  2:9269491/96‑1 (MQ=255)
                                                      aTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc                                 >  1:4404202/1‑96 (MQ=255)
                                                       tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                                >  1:4764825/1‑96 (MQ=255)
                                                       tGTCCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa                                <  1:9397057/96‑1 (MQ=255)
                                                         ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                              >  1:10017308/1‑96 (MQ=255)
                                                         ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                              <  1:13095262/96‑1 (MQ=255)
                                                         ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTGCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                              >  2:3600059/1‑96 (MQ=255)
                                                           cACCAGCGGGTTAATAAAGTACCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc                            <  1:299901/96‑1 (MQ=255)
                                                           cACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc                            >  1:187448/1‑96 (MQ=255)
                                                            ccaaGCGGGTTAATAAAGTAACCAAGGCCTCCTTCCAGCATATGGTGATTTTCCACCCCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg                           <  1:6601348/93‑1 (MQ=255)
                                                            aCCAGCGGGTTAATAAAGTAATCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg                           >  1:12851332/1‑96 (MQ=255)
                                                             ccAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg                          >  1:6841723/1‑96 (MQ=255)
                                                              cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                         >  1:2529544/1‑96 (MQ=255)
                                                              cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGc                         <  1:7393924/96‑1 (MQ=255)
                                                               aGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                        <  2:14303102/96‑1 (MQ=255)
                                                                 cGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaga                      <  1:3754276/96‑1 (MQ=255)
                                                                  gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCaggg                     <  1:10688695/96‑3 (MQ=255)
                                                                   ggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGGGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAgaga                    >  1:623672/1‑96 (MQ=255)
                                                                     ttAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACt                  <  2:709077/96‑1 (MQ=255)
                                                                       aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGc                <  1:8625399/96‑1 (MQ=255)
                                                                       aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGc                >  2:2268371/1‑96 (MQ=255)
                                                                         taaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCa              <  2:934882/96‑1 (MQ=255)
                                                                         taaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCa              >  1:10092753/1‑96 (MQ=255)
                                                                          aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg             <  1:3278438/96‑1 (MQ=255)
                                                                          aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg             >  2:8520476/1‑96 (MQ=255)
                                                                          aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAg             <  1:1195676/96‑1 (MQ=255)
                                                                           aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGc            >  2:8867959/1‑96 (MQ=255)
                                                                           aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGc            <  2:12244063/96‑1 (MQ=255)
                                                                           aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTCCACCGCCCAGATAACCAGTAGCCAGTTCCCGCCAATCAGCACGGCAGAGACTGCCAGc            <  1:7641317/96‑1 (MQ=255)
                                                                              tAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATa         >  2:1123458/1‑96 (MQ=255)
                                                                                aCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaa       <  1:2232887/96‑1 (MQ=255)
                                                                                 ccAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaa      <  2:268203/96‑1 (MQ=255)
                                                                                  cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATaaaaa     >  2:3541023/1‑96 (MQ=255)
                                                                                    aGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAtt   >  2:7565592/1‑96 (MQ=255)
                                                                                     ggCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAttt  >  2:3928292/1‑96 (MQ=255)
                                                                                          ||                                                                                         
GCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTT  >  NC_000913/4001396‑4001576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: