Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,528 T→C Q8R (CAG→CGG)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,5280TC99.5% 677.7 / ‑5.2 212Q8R (CAG→CGG) araDL‑ribulose‑5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (106/105);  total (106/106)
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.89e-01

Only 100 of 222 total aligned reads displayed.

ACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGAT  >  NC_000913/66438‑66617
                                                                                          |                                                                                         
 cgccgcGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGc                                                                                        <  2:5603457/93‑1 (MQ=255)
  gccgcGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCg                                                                                       >  1:4481863/1‑93 (MQ=255)
    cgcgCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGAGCCAGGTTGGCTTCTAATACCCGGCGcc                                                                                     >  2:5036762‑M1/1‑91 (MQ=255)
     gcgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGccc                                                                                    <  1:1159972‑M1/93‑4 (MQ=255)
     gcgcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGccc                                                                                    <  2:422402‑M1/93‑4 (MQ=255)
       gcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccga                                                                                  >  2:5523169‑M1/1‑88 (MQ=255)
        cTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGAAGCGCCAGGTTGGCTTCTAATATCTGGCGcccgag                                                                                 >  1:2810264‑M1/1‑87 (MQ=255)
         tCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgaga                                                                                >  1:560316‑M1/1‑86 (MQ=255)
           gcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacc                                                                              <  1:2887097‑M1/93‑10 (MQ=255)
             gATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccg                                                                            >  1:5270888‑M1/1‑82 (MQ=255)
              aTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgc                                                                           >  1:910244‑M1/1‑81 (MQ=255)
              aTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgc                                                                           >  2:5991316‑M1/1‑81 (MQ=255)
                  aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgc                                                                       >  1:6163792‑M1/1‑77 (MQ=255)
                  aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgc                                                                       <  1:6028280‑M1/93‑17 (MQ=255)
                  aCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgc                                                                       >  2:6176286‑M1/1‑77 (MQ=255)
                    ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                     >  1:757731‑M1/1‑75 (MQ=255)
                    ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                     >  2:293999‑M1/1‑75 (MQ=255)
                    ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                     >  1:5602247‑M1/1‑75 (MQ=255)
                    ggCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgt                                                                     >  2:4430481‑M1/1‑75 (MQ=255)
                       gcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcga                                                                  >  2:428435‑M1/1‑72 (MQ=255)
                       gcTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcga                                                                  <  1:2765620‑M1/93‑22 (MQ=255)
                        cTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgac                                                                 >  2:1917294‑M1/1‑71 (MQ=255)
                         tcacgTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacc                                                                <  2:16507‑M1/91‑24 (MQ=255)
                         tGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgcgacccgcacgcgtcgacc                                                                >  1:2874501‑M1/1‑70 (MQ=255)
                          gACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacct                                                               >  2:3680735‑M1/1‑69 (MQ=255)
                            cGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgc                                                             <  1:2065897‑M1/93‑27 (MQ=255)
                            cGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGNTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgc                                                             <  1:5359083‑M1/93‑27 (MQ=255)
                            cGTCGCGCCATGTGTGCGTGGCTAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgc                                                             <  2:712547‑M1/93‑27 (MQ=255)
                             gTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgca                                                            <  2:3240094‑M1/93‑28 (MQ=255)
                              ttGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcag                                                           >  2:5342245‑M1/1‑65 (MQ=255)
                                 ccccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggta                                                        <  1:5223574‑M1/93‑32 (MQ=255)
                                  cccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtat                                                       >  2:4797847‑M1/1‑61 (MQ=255)
                                   ccATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatg                                                      >  1:5660980‑M1/1‑60 (MQ=255)
                                    cATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatgg                                                     <  2:1565080‑M1/93‑35 (MQ=255)
                                    cATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatgg                                                     <  1:5458760‑M1/93‑35 (MQ=255)
                                    cATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatgg                                                     <  1:3081392‑M1/93‑35 (MQ=255)
                                     aTGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatgga                                                    >  1:4752025‑M1/1‑58 (MQ=255)
                                     aTGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatgga                                                    <  1:2043569‑M1/93‑36 (MQ=255)
                                      tgtgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggag                                                   >  2:1168468‑M1/1‑57 (MQ=255)
                                      tgtgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggag                                                   >  2:2155769‑M1/1‑57 (MQ=255)
                                      tgtgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggag                                                   >  2:5668767‑M1/1‑57 (MQ=255)
                                        tgAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGATGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagga                                                 >  2:4648506‑M1/1‑55 (MQ=255)
                                         gAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaa                                                >  2:3900748‑M1/1‑54 (MQ=255)
                                          aGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaac                                               >  1:3268496‑M1/1‑53 (MQ=255)
                                           gtgtgGCCAGGTTCTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaaca                                              <  2:4054489‑M1/91‑42 (MQ=255)
                                           gCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagtaaca                                              >  1:3941090‑M1/1‑52 (MQ=255)
                                            cGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacag                                             <  1:1096801‑M1/93‑43 (MQ=255)
                                             gTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                            <  1:4242052‑M1/93‑44 (MQ=255)
                                             gTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                            <  1:2459711‑M1/93‑44 (MQ=255)
                                             gTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagt                                            <  1:1499644‑M1/93‑44 (MQ=255)
                                              tGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagta                                           <  2:4930788‑M1/93‑45 (MQ=255)
                                               gACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtag                                          <  1:2928775‑M1/93‑46 (MQ=255)
                                                aCCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtaga                                         >  2:5257755‑M1/1‑47 (MQ=255)
                                                aCCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtaga                                         >  1:629832‑M1/1‑47 (MQ=255)
                                                aCCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtaga                                         <  1:5163410‑M1/93‑47 (MQ=255)
                                                 ccAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagag                                        >  1:1856770‑M1/1‑46 (MQ=255)
                                                 ccAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagag                                        >  1:6139271‑M1/1‑46 (MQ=255)
                                                 ccAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagag                                        <  1:4611874‑M1/93‑48 (MQ=255)
                                                 ccAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagag                                        <  1:4869003‑M1/93‑48 (MQ=255)
                                                  cAGGTTGTGTTGTGGCAGCGCGAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagaga                                       <  1:4975485‑M1/93‑49 (MQ=255)
                                                   aGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagag                                      >  2:416897‑M1/1‑44 (MQ=255)
                                                    ggTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagt                                     <  1:4581845‑M1/93‑51 (MQ=255)
                                                     gTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagtt                                    <  2:737145‑M1/93‑52 (MQ=255)
                                                     gTTGTGTTTTGGCAGCGCCAGGTTGCCTTCTAATACCCGGCGcccgagacccgcacgcgtcgccctgcaggtatggagaaacagtagagagtt                                    <  1:2275044‑M1/93‑52 (MQ=255)
                                                        gtgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcg                                 <  1:2612530‑M1/93‑55 (MQ=255)
                                                        gtgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcg                                 >  2:5190890‑M1/1‑39 (MQ=255)
                                                         tgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcga                                >  1:1976533‑M1/1‑38 (MQ=255)
                                                         tgtTTTGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgacgtcctgtaggtatggagaaacagtagagagttgaga                                <  2:6261014‑M1/93‑56 (MQ=255)
                                                           ttttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgata                              >  2:292822‑M1/1‑36 (MQ=255)
                                                            tttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataa                             <  1:1957382‑M1/93‑59 (MQ=255)
                                                             ttGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaa                            <  1:4998036‑M1/93‑60 (MQ=255)
                                                              tGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaa                           >  2:3338118‑M1/1‑33 (MQ=255)
                                                              tGGCAGCGCCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaa                           <  1:992298‑M1/93‑61 (MQ=255)
                                                                   gcgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgt                      >  2:3559884‑M1/1‑28 (MQ=255)
                                                                    cgcCAGGTTGGCTTCTAATACCCGGCGccggagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtc                     <  1:5253349‑M1/93‑67 (MQ=255)
                                                                    cgcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtc                     <  2:3018440‑M1/93‑67 (MQ=255)
                                                                     gcCAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtca                    <  2:2679261‑M1/93‑68 (MQ=255)
                                                                      ccAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcag                   >  1:3420237‑M1/1‑25 (MQ=255)
                                                                      ccAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcag                   <  1:317698‑M1/93‑69 (MQ=255)
                                                                      ccAGGTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcag                   >  2:6024458‑M1/1‑25 (MQ=255)
                                                                         ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                <  2:198239‑M1/93‑72 (MQ=255)
                                                                         ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                >  1:4214976‑M1/1‑22 (MQ=255)
                                                                         ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                >  2:674400‑M1/1‑22 (MQ=255)
                                                                         ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                >  2:2400686‑M1/1‑22 (MQ=255)
                                                                         ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                <  1:1301996‑M1/93‑72 (MQ=255)
                                                                         ggTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggta                >  1:2197756‑M1/1‑22 (MQ=255)
                                                                          gTTGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtag               >  2:249836‑M1/1‑21 (MQ=255)
                                                                           ttGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagg              >  2:4295400‑M1/1‑20 (MQ=255)
                                                                           ttGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagg              <  2:2606729‑M1/93‑74 (MQ=255)
                                                                            tGGCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtagga             <  2:5784674‑M1/93‑75 (MQ=255)
                                                                             ggCTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggat            <  2:5152207‑M1/93‑76 (MQ=255)
                                                                               cTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatcc          >  1:5781562‑M1/1‑16 (MQ=255)
                                                                               cTTCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatcc          >  1:5178503‑M1/1‑16 (MQ=255)
                                                                                ttCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagggagttgcgataaaaagcgtcagggaggatccg         >  2:272997‑M1/1‑15 (MQ=255)
                                                                                 tCTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgc        <  1:4165755‑M1/93‑80 (MQ=255)
                                                                                  cTAATACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgct       >  1:2535108‑M1/1‑13 (MQ=255)
                                                                                    aaTACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaa     <  2:3460444‑M1/93‑83 (MQ=255)
                                                                                     aTACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaat    >  2:3524339‑M1/1‑10 (MQ=255)
                                                                                      tACCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatc   >  1:6215100‑M1/1‑9 (MQ=255)
                                                                                       aCCCGGCGcccgagacccgcacgcgtcgacctgcaggtatggagaaacagtagagagttgcgataaaaagcgtcaggtaggatccgctaatct  >  1:3494510‑M1/1‑8 (MQ=255)
                                                                                          |                                                                                         
ACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGAT  >  NC_000913/66438‑66617

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