Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,992 A→G E117G (GAG→GGG)  wzzE → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,9920AG98.6% 456.6 / ‑2.4 146E117G (GAG→GGG) wzzEputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (68/75);  total (68/78)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 153 total aligned reads displayed.

CAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACA  >  NC_000913/3966902‑3967081
                                                                                          |                                                                                         
  aaCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTc                                                                                       >  1:1345357/1‑93 (MQ=255)
  aaCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTc                                                                                       <  1:3474373/93‑1 (MQ=255)
   aCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCt                                                                                      >  1:466698/1‑93 (MQ=255)
   aCATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCt                                                                                      <  1:2800390/93‑1 (MQ=255)
    cATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTg                                                                                     <  2:2403217/93‑1 (MQ=255)
     aTGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTgg                                                                                    <  1:2089374/93‑1 (MQ=255)
      tGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGc                                                                                   <  1:3134576/93‑1 (MQ=255)
      tGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGc                                                                                   <  2:3941621/93‑1 (MQ=255)
      tGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGc                                                                                   <  2:3180144/93‑1 (MQ=255)
       ggCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCt                                                                                  >  2:2266120/1‑93 (MQ=255)
       ggCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCt                                                                                  >  1:1119050/1‑93 (MQ=255)
             tGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAAc                                                                            <  2:3157500/93‑1 (MQ=255)
             tGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAAc                                                                            >  1:3751588/1‑93 (MQ=255)
              gCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAAcc                                                                           >  1:578214/1‑93 (MQ=255)
                 gACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAAATGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGat                                                                        >  1:4164658/1‑92 (MQ=255)
                  accaaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACt                                                                       >  2:2019024/1‑93 (MQ=255)
                    caaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTAt                                                                     >  2:3688097/1‑93 (MQ=255)
                     aaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTAtt                                                                    >  1:80561/1‑93 (MQ=255)
                     aaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTAtt                                                                    >  2:3037722/1‑93 (MQ=255)
                     aaccaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTAtt                                                                    >  2:2728485/1‑93 (MQ=255)
                      accaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTa                                                                   >  2:4153489/1‑93 (MQ=255)
                      accaTCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTa                                                                   <  2:296643/93‑1 (MQ=255)
                      accaGCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTa                                                                   <  2:296644/93‑1 (MQ=255)
                          tCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACaaa                                                               >  2:1659819/1‑93 (MQ=255)
                          tCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACaaa                                                               >  1:3347883/1‑93 (MQ=255)
                              tCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATCCCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGc                                                           <  2:2225203/93‑1 (MQ=255)
                              tCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTCAAAACAGc                                                           >  2:3435607/1‑93 (MQ=255)
                              tCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGc                                                           >  2:4085473/1‑93 (MQ=255)
                              tCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCAGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGc                                                           <  1:4578854/93‑1 (MQ=255)
                              gcATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGc                                                           <  1:3525584/92‑1 (MQ=255)
                                aTGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCgg                                                         <  1:1275276/93‑1 (MQ=255)
                                aTGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCgg                                                         >  2:3584933/1‑93 (MQ=255)
                                 tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGa                                                        >  1:3441153/1‑93 (MQ=255)
                                 tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGa                                                        >  2:4382370/1‑93 (MQ=255)
                                 tGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGa                                                        <  1:4481409/93‑1 (MQ=255)
                                  ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAt                                                       <  2:1324865/93‑1 (MQ=255)
                                  ggACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAt                                                       <  2:4370891/93‑1 (MQ=255)
                                    aCGAAGTCTATAAAGAGTTTGTTATGGAAGTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtgg                                                     <  1:2754532/93‑1 (MQ=255)
                                    aCGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtgg                                                     <  1:2754534/93‑1 (MQ=255)
                                     cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtggt                                                    >  2:731219/1‑93 (MQ=255)
                                     cGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtggt                                                    <  1:4562071/93‑1 (MQ=255)
                                      gAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtggtg                                                   >  2:2861804/1‑93 (MQ=255)
                                      gAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtggtg                                                   <  2:4290357/93‑1 (MQ=255)
                                      gAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtggtg                                                   <  1:2855250/93‑1 (MQ=255)
                                       aaGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGAtggtgg                                                  >  2:4313033/1‑93 (MQ=255)
                                        aGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTggg                                                 <  1:2264414/93‑1 (MQ=255)
                                        aGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTggg                                                 >  1:3969198/1‑93 (MQ=255)
                                         gCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTggga                                                >  1:2649787/1‑92 (MQ=255)
                                         gCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGc                                                <  1:2363220/93‑1 (MQ=255)
                                          ctgataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCa                                               <  1:4487988/90‑1 (MQ=255)
                                           cTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCaa                                              >  1:3536007/1‑93 (MQ=255)
                                            tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCACACCGACTATTACAAACAGCGGATGGTGGGCAAc                                             >  2:448939/1‑93 (MQ=255)
                                            tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAAc                                             >  2:349965/1‑93 (MQ=255)
                                            tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAAc                                             <  1:3509625/93‑1 (MQ=255)
                                            tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAAc                                             <  2:1900865/93‑1 (MQ=255)
                                            tataAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAAc                                             >  1:3421727/1‑93 (MQ=255)
                                               aaaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGc                                          <  1:1128922/93‑1 (MQ=255)
                                                aaGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCa                                         <  2:4451405/93‑1 (MQ=255)
                                                   agTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAg                                      <  2:1583160/93‑1 (MQ=255)
                                                   agTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAg                                      >  2:2677296/1‑93 (MQ=255)
                                                   agTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAg                                      <  1:694740/93‑1 (MQ=255)
                                                   agTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAg                                      <  2:2239606/93‑1 (MQ=255)
                                                    gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGc                                     >  1:1579751/1‑93 (MQ=255)
                                                    gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGc                                     <  2:4233552/93‑1 (MQ=255)
                                                    gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGc                                     <  2:2805840/93‑1 (MQ=255)
                                                    gTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGCCAAAAGAAAAGc                                     <  2:2347779/93‑1 (MQ=255)
                                                     tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAgac                                    >  1:3324174/1‑91 (MQ=255)
                                                     tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGcc                                    <  1:29243/93‑1 (MQ=255)
                                                     tttGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGcc                                    <  2:3220930/93‑1 (MQ=255)
                                                       tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGa                                  <  1:1512439/93‑1 (MQ=255)
                                                       tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGa                                  <  1:1025053/93‑1 (MQ=255)
                                                       tGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGa                                  <  1:3266141/93‑1 (MQ=255)
                                                         gtATGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATg                                <  2:2062736/92‑1 (MQ=255)
                                                           aTGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcg                              >  1:2953456/1‑93 (MQ=255)
                                                            tGCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcgg                             >  2:682747/1‑93 (MQ=255)
                                                             gCAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcggc                            <  1:724815/93‑1 (MQ=255)
                                                              cAACTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATgcggcg                           <  2:1184778/93‑1 (MQ=255)
                                                               aaCTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCACCAGCAACGCCGATGCGGCGt                          >  2:1328274/1‑93 (MQ=255)
                                                               aaCTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGt                          <  1:2106702/93‑1 (MQ=255)
                                                               aaCTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGt                          >  1:484711/1‑93 (MQ=255)
                                                                aCTGGCCTCGTGTTATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGtt                         <  2:2293964/93‑1 (MQ=255)
                                                                aCTGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGtt                         <  1:1491608/93‑1 (MQ=255)
                                                                aCTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACAGACTATTACAAACATCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGtt                         <  1:967438/93‑1 (MQ=255)
                                                                  tGGCCTCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGc                       <  2:2979817/93‑1 (MQ=255)
                                                                   ggCCTCGTGGGATACCCGCAGAGCGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCt                      <  1:1606369/93‑1 (MQ=255)
                                                                       tCGTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGAt                  >  2:2118639/1‑93 (MQ=255)
                                                                         gTGGGATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGa                >  2:151304/1‑93 (MQ=255)
                                                                           gggATACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGaaa              >  1:4136327/1‑93 (MQ=255)
                                                                            ggATACCCGCAGAGGGTTCTGTCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAAt             >  1:850846/1‑93 (MQ=255)
                                                                             gATACCCGCAGAGGTTGCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATg            <  2:1287698/93‑1 (MQ=255)
                                                                              aTACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGa           <  2:1437605/93‑1 (MQ=255)
                                                                              aTACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGa           <  1:2270627/93‑1 (MQ=255)
                                                                               tACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt          <  2:1292105/93‑1 (MQ=255)
                                                                               tACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt          >  1:2483219/1‑93 (MQ=255)
                                                                               tACCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAt          >  2:4547667/1‑93 (MQ=255)
                                                                                aCCCGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGAtt         >  2:4718954/1‑93 (MQ=255)
                                                                                 cccGCAGAGGGTTCTGGCTGGAAACCGACTATTACAAACAGCGGATGGTGGGCAATAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTa        <  1:4068270/93‑1 (MQ=255)
                                                                                 cccGCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGTAATGATTa        >  2:84137/1‑93 (MQ=255)
                                                                                    gCAGAGGGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaaca     <  1:1025682/93‑1 (MQ=255)
                                                                                      agagggTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTaacaac   >  2:1773641/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAAACATGGCTTCTGCCGACCAACCATCGGTCATGGACGAAGCCTATAAAGAGTTTGTTATGCAACTGGCCTCGTGGGATACCCGCAGAGAGTTCTGGCTGCAAACCGACTATTACAAACAGCGGATGGTGGGCAACAGCAAAGCCGATGCGGCGTTGCTGGATGAAATGATTAACAACA  >  NC_000913/3966902‑3967081

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