Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,203,421 A→G V30A (GTT→GCT)  purD ← phosphoribosylglycinamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,203,4210AG99.6% 828.3 / ‑5.0 261V30A (GTT→GCT) purDphosphoribosylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (143/117);  total (144/117)
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.15e-01

Only 100 of 281 total aligned reads displayed.

CGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATC  >  NC_000913/4203332‑4203510
                                                                                         |                                                                                         
cGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGcagcaa                                                                                        <  2:7358960/93‑1 (MQ=255)
  aaaTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCaaaa                                                                                      >  1:1980007/1‑93 (MQ=255)
   aaTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCaaaaa                                                                                     >  1:81788/1‑93 (MQ=255)
   aaTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCaaaaa                                                                                     <  1:454687/93‑1 (MQ=255)
   aaTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCaaaaa                                                                                     <  1:2056157/93‑1 (MQ=255)
     tCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACa                                                                                   >  1:3404681/1‑93 (MQ=255)
     tCCAACAGCGCCGGGATATCGGTCACGCCAATACCAACGTTTTACAGCGCGGGGTCCAAGGCAGGGCCTGCATTACCCGGCGCCGCAAAAACa                                                                                   >  2:5985908/1‑93 (MQ=255)
      ccAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAg                                                                                  >  1:508887/1‑93 (MQ=255)
       cAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGt                                                                                 >  2:6549500/1‑93 (MQ=255)
       cAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGt                                                                                 <  2:6040928/93‑1 (MQ=255)
        aaCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtc                                                                                >  2:3145918/1‑93 (MQ=255)
         aCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtct                                                                               >  2:106199/1‑93 (MQ=255)
         aCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtct                                                                               >  2:5137821/1‑93 (MQ=255)
           aGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCa                                                                             <  1:3722742/93‑1 (MQ=255)
           aGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCa                                                                             <  1:7363943/93‑1 (MQ=255)
           aGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCa                                                                             >  1:2197258/1‑93 (MQ=255)
           aGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCa                                                                             <  1:1717418/93‑1 (MQ=255)
            gcgcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCaa                                                                            >  1:6948352/1‑93 (MQ=255)
               ccGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCa                                                                         <  1:6775841/93‑1 (MQ=255)
               ccGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCa                                                                         >  1:722854/1‑93 (MQ=255)
                 gggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgc                                                                       >  2:4194499/1‑93 (MQ=255)
                  ggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcg                                                                      <  1:3300480/93‑1 (MQ=255)
                   gATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcgg                                                                     <  2:721008/93‑1 (MQ=255)
                    atatCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcggc                                                                    >  1:364151/1‑93 (MQ=255)
                     tatCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcggcg                                                                   >  1:2351205/1‑93 (MQ=255)
                     tatCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAgcggcg                                                                   >  1:6851039/1‑93 (MQ=255)
                      atCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGa                                                                  <  2:263430/93‑1 (MQ=255)
                      atCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGa                                                                  <  1:1417909/93‑1 (MQ=255)
                      atCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGa                                                                  <  2:7576080/93‑1 (MQ=255)
                       tCGGTCACGCCAATAGCAACGTTTTTCAGCGCGTGTTGCAGTGCACTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGAc                                                                 >  2:6368617/1‑93 (MQ=255)
                       tCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGAc                                                                 <  1:5264327/93‑1 (MQ=255)
                       tCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGAc                                                                 <  1:5709944/93‑1 (MQ=255)
                             aCGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGcc                                                           >  2:334863/1‑93 (MQ=255)
                             aCGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGcc                                                           <  1:417372/93‑1 (MQ=255)
                                 cAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCtt                                                       >  2:1852134/1‑93 (MQ=255)
                                   aTAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTc                                                     >  2:6478835/1‑93 (MQ=255)
                                    tAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTcc                                                    >  1:2639958/1‑93 (MQ=255)
                                    tAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTcc                                                    <  1:749289/93‑1 (MQ=255)
                                     aGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCa                                                   <  2:3628188/93‑1 (MQ=255)
                                     aGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCa                                                   <  1:1789847/93‑1 (MQ=255)
                                      gCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAg                                                  >  2:1487878/1‑93 (MQ=255)
                                      gCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAg                                                  >  1:7661765/1‑93 (MQ=255)
                                        aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGc                                                >  1:4098560/1‑93 (MQ=255)
                                        aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGc                                                >  2:4016464/1‑93 (MQ=255)
                                        aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGc                                                <  2:2313630/93‑1 (MQ=255)
                                        aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGc                                                >  2:4767911/1‑93 (MQ=255)
                                         aCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGcc                                               >  1:5670802/1‑93 (MQ=255)
                                         aCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGcc                                               <  1:5947504/93‑1 (MQ=255)
                                          cGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCa                                              <  1:3965960/93‑1 (MQ=255)
                                          cGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCa                                              <  2:810854/93‑1 (MQ=255)
                                          cGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCa                                              <  1:6718166/93‑1 (MQ=255)
                                           gTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAg                                             >  2:5074650/1‑93 (MQ=255)
                                            ttttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCCGGCCAgc                                            >  2:6323792/1‑93 (MQ=255)
                                            ttttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAgc                                            <  1:707815/93‑1 (MQ=255)
                                            ttttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAgc                                            <  1:6215432/93‑1 (MQ=255)
                                              ttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAgcgc                                          >  2:6181062/1‑93 (MQ=255)
                                               tGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGcgcg                                         <  2:2840426/93‑1 (MQ=255)
                                                 cAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTg                                       <  1:4124865/93‑1 (MQ=255)
                                                     gcgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTcg                                   >  2:1289554/1‑93 (MQ=255)
                                                     gcgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTcc                                   >  1:5393269/1‑92 (MQ=255)
                                                     gcgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTAGGCCGCTTTCCAGGCCAGCGCGTGCTcg                                   >  2:6127218/1‑93 (MQ=255)
                                                         gTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcc                               >  1:5832875/1‑93 (MQ=255)
                                                          ttCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgccc                              >  2:5260207/1‑93 (MQ=255)
                                                             cagtgcagtgCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccgcc                           >  1:939345/1‑93 (MQ=255)
                                                             cagtgcagtgCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccgcc                           >  2:5647813/1‑93 (MQ=255)
                                                               gtgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGt                         <  2:203746/93‑1 (MQ=255)
                                                                tgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCACCCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGtt                        <  2:7775083/93‑1 (MQ=255)
                                                                tgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGtt                        <  2:306519/93‑1 (MQ=255)
                                                                tgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGtt                        >  2:93490/1‑93 (MQ=255)
                                                                tgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGtt                        <  1:2595314/93‑1 (MQ=255)
                                                                 gcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTa                       >  2:1227249/1‑93 (MQ=255)
                                                                  cagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTAc                      >  1:7891913/1‑93 (MQ=255)
                                                                   agtgcCTGCATTCCCCGGAGCAGCAAAAGCAGTCTCAACCAGCGGCGACTGGGCCGTTTTCCAGCCCAGCGCGTGCTCGCGCCCGCCGTTAcc                     <  2:2145752/93‑1 (MQ=255)
                                                                   agtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTAcc                     >  1:1564662/1‑93 (MQ=255)
                                                                   agtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTAcc                     >  1:106148/1‑93 (MQ=255)
                                                                   agtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTAcc                     >  2:1707365/1‑93 (MQ=255)
                                                                     tgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCaa                   >  1:7334293/1‑93 (MQ=255)
                                                                     tgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCaa                   <  2:4844922/93‑1 (MQ=255)
                                                                      gcCTGCATTCCCCGGAGCAGAAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAAt                  <  2:4492500/93‑1 (MQ=255)
                                                                       ccTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATc                 <  1:120163/93‑1 (MQ=255)
                                                                       ccTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATc                 <  1:3053682/93‑1 (MQ=255)
                                                                       ccTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATc                 <  1:3812538/93‑1 (MQ=255)
                                                                       ccTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATc                 >  2:4054317/1‑93 (MQ=255)
                                                                        cTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCa                >  2:297335/1‑93 (MQ=255)
                                                                         tGCATTACCCGGAGCGGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCAc               >  2:3539979/1‑93 (MQ=255)
                                                                         tGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCAc               <  1:1531034/93‑1 (MQ=255)
                                                                         tGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGCCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCAc               <  2:69491/93‑1 (MQ=255)
                                                                            aTTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTaa            >  1:4888455/1‑93 (MQ=255)
                                                                            aTTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACGGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTaa            <  1:1314489/93‑1 (MQ=255)
                                                                             ttACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAAt           <  2:2743266/93‑1 (MQ=255)
                                                                             ttACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAAt           <  1:7119643/93‑1 (MQ=255)
                                                                             ttACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAAt           >  1:7589625/1‑93 (MQ=255)
                                                                                 ccGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACtt       <  2:6673145/93‑1 (MQ=255)
                                                                                  cgaagcagcAAAACCAGCCTGACCCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACttt      <  2:3126681/90‑1 (MQ=255)
                                                                                  cGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACttt      >  1:4838331/1‑93 (MQ=255)
                                                                                   ggAGCAGCAAAAACAGTCTCAACTAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTc     <  1:4741756/93‑1 (MQ=255)
                                                                                   ggAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGCCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTc     <  1:630202/93‑1 (MQ=255)
                                                                                    gAGCAGCAAAAGCAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCa    <  2:2442910/93‑1 (MQ=255)
                                                                                    gAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGCCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCa    <  1:7763889/93‑1 (MQ=255)
                                                                                      gcagcaAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATc  <  2:4565953/93‑1 (MQ=255)
                                                                                         |                                                                                         
CGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATC  >  NC_000913/4203332‑4203510

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