Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 526,302 (G)6→7 coding (3818/4281 nt) rhsD → RhsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913526,2961.G89.4% 463.9 / 50.5 160coding (3812/4281 nt)rhsDRhsD
Reads supporting (aligned to +/- strand):  ref base . (6/11);  new base G (79/64);  total (85/75)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAAC  >  NC_000913/526204‑526386
                                                                                              |                                                                                           
gAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAc                                                                                              >  1:978666/1‑93 (MQ=255)
gAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAc                                                                                              <  1:4655229/93‑1 (MQ=255)
 aaTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑g                                                                                            <  2:206517/93‑1 (MQ=255)
 aaTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑g                                                                                            >  1:4517420/1‑93 (MQ=255)
  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gg                                                                                           <  1:1868790/93‑1 (MQ=255)
  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gg                                                                                           <  1:3381714/93‑1 (MQ=255)
  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gg                                                                                           <  1:2384502/93‑1 (MQ=255)
   tttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTCTTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑ggg                                                                                          <  2:6114485/93‑1 (MQ=255)
   tttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑ggg                                                                                          <  1:3515911/93‑1 (MQ=255)
   tttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑ggg                                                                                          <  2:1653261/93‑1 (MQ=255)
       tatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑ggggggg                                                                                      <  2:4936107/93‑2 (MQ=255)
        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gggggggg                                                                                     <  2:2077551/93‑3 (MQ=255)
         tCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCAAGATCTGGAGCA‑TGTAC‑ggggggggg                                                                                    <  2:5025511/93‑4 (MQ=255)
          cAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gggggggggt                                                                                   >  2:6269773/1‑89 (MQ=255)
          cAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gggggggggt                                                                                   <  2:568345/93‑5 (MQ=255)
          cAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑gggggggggt                                                                                   <  2:830804/93‑5 (MQ=255)
          cAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGG‑Gtt                                                                                  <  2:4802918/93‑1 (MQ=255)
               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTtgtg                                                                              >  2:5813869/1‑93 (MQ=255)
               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTtgtg                                                                              >  2:5418741/1‑93 (MQ=255)
               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGG‑GTTTgtgg                                                                             >  1:307468/1‑93 (MQ=255)
                ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTgtgg                                                                             >  2:4248101/1‑93 (MQ=255)
                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTgtggtg                                                                           >  2:3759852/1‑93 (MQ=255)
                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTgtgggg                                                                           >  2:3297235/1‑91 (MQ=255)
                   ttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGgtgt                                                                          >  1:1764080/1‑93 (MQ=255)
                   ttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGgtgt                                                                          >  1:5102395/1‑93 (MQ=255)
                   ttAAATCCACTAAAAAAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGgtgt                                                                          <  2:3549259/93‑1 (MQ=255)
                      aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCt                                                                       <  1:2634312/93‑1 (MQ=255)
                       aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGGTCtt                                                                      >  2:5883803/1‑93 (MQ=255)
                        tCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCttt                                                                     >  2:3927865/1‑93 (MQ=255)
                           aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCAc                                                                  <  2:326099/93‑1 (MQ=255)
                           aCTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GGTTGTGGTGTTCTTTCAc                                                                  >  2:1450725/1‑93 (MQ=255)
                            cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACg                                                                 >  1:3148582/1‑93 (MQ=255)
                            cTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACg                                                                 >  1:6113989/1‑93 (MQ=255)
                             tACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGt                                                                >  2:2740209/1‑93 (MQ=255)
                             tACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGt                                                                >  2:5219910/1‑93 (MQ=255)
                               caacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGtat                                                              >  1:3703236/1‑93 (MQ=255)
                                aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGtata                                                             >  2:6199105/1‑93 (MQ=255)
                                aacaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGGGGTTTGTGGTGTTCTTTCACGtata                                                             <  2:5922601/93‑1 (MQ=255)
                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTataa                                                            >  2:102034/1‑93 (MQ=255)
                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTataa                                                            >  1:1173173/1‑93 (MQ=255)
                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTataa                                                            >  1:5340545/1‑93 (MQ=255)
                                  caaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTataat                                                           <  1:866925/93‑1 (MQ=255)
                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAg                                                         <  1:5030252/93‑1 (MQ=255)
                                      ttGCCCTTATGGGATTATGGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGa                                                       <  1:2541220/93‑1 (MQ=255)
                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGGGTTCTTTCACGTATAATAGGa                                                       >  1:1299379/1‑93 (MQ=255)
                                       tGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGAc                                                      <  2:1027022/93‑1 (MQ=255)
                                       tGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGAc                                                      >  1:3890187/1‑93 (MQ=255)
                                        gACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATATGAcc                                                     >  2:4349823/1‑93 (MQ=255)
                                         gctgtATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCa                                                    >  2:3163456/5‑93 (MQ=255)
                                          cccTATGGGATTATTGCAGACTTGGGAAGATGCCAGGTCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                                   <  2:2022366/93‑1 (MQ=255)
                                          cccTATGGGATTATAGCAGATTTGGGATGATGCCAGATCTGGAGCA‑TGTCCGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                                   <  2:2375927/93‑1 (MQ=255)
                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                >  1:1917284/1‑93 (MQ=255)
                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                >  1:2164481/1‑93 (MQ=255)
                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GCTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                                >  2:927185/1‑93 (MQ=255)
                                              aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGG‑GCTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                              >  2:1755225/1‑93 (MQ=255)
                                               tGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                              >  1:3823042/1‑93 (MQ=255)
                                                gggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAt                                             <  2:5589483/93‑1 (MQ=255)
                                                 ggATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAtt                                            >  1:5851404/1‑93 (MQ=255)
                                                   attattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                          <  1:4190405/93‑1 (MQ=255)
                                                    ttattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATACTAGGACCAAGTAAATTTGa                                         >  1:593476/1‑93 (MQ=255)
                                                    ttattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGa                                         <  1:1229032/93‑1 (MQ=255)
                                                     tattGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGAt                                        <  1:3773577/93‑1 (MQ=255)
                                                      attGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATa                                       >  2:1486245/1‑93 (MQ=255)
                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                      <  2:766363/93‑1 (MQ=255)
                                                         gCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                    >  2:2257563/1‑93 (MQ=255)
                                                         gCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                    <  1:909482/93‑1 (MQ=255)
                                                         gCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                    <  1:4675748/93‑1 (MQ=255)
                                                         gCAGACTTGGGATGATGCCAGAGCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                    <  1:2126251/93‑1 (MQ=255)
                                                           agtgttGGGATGAGGCCAGATGTGGACCAGTGCGG‑GGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                                  <  1:1169877/89‑1 (MQ=255)
                                                           aGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                                 <  2:50422/93‑1 (MQ=255)
                                                             aCTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                               >  1:5826917/1‑93 (MQ=255)
                                                              cTTGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                               >  1:1863056/1‑93 (MQ=255)
                                                               ttGGGATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                              >  2:4180404/1‑93 (MQ=255)
                                                                 gggATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAt                            >  2:1844925/1‑93 (MQ=255)
                                                                 gggATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAt                            >  1:5495235/1‑93 (MQ=255)
                                                                  ggATGATGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtg                           <  2:2517080/93‑1 (MQ=255)
                                                                   gatgatGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                          >  2:3066013/1‑93 (MQ=255)
                                                                   gatgatGCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                          <  1:2962357/93‑1 (MQ=255)
                                                                    atgatgCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                         >  1:4004988/1‑93 (MQ=255)
                                                                       atgCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                      >  1:5241005/1‑93 (MQ=255)
                                                                       atgCCAGATCGGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                      <  1:6248862/93‑1 (MQ=255)
                                                                        tgCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGt                     >  1:1451079/1‑93 (MQ=255)
                                                                         gCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGtt                    >  1:3952365/1‑93 (MQ=255)
                                                                         gCCAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGtt                    >  1:792698/1‑93 (MQ=255)
                                                                          ccAGATCTGGAGCA‑TGTACGGGGGGGT‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                   <  1:3056042/93‑1 (MQ=255)
                                                                          ccAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                   >  2:4050355/1‑93 (MQ=255)
                                                                          ccAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                   >  2:89575/1‑93 (MQ=255)
                                                                           cAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                  <  1:4988427/93‑1 (MQ=255)
                                                                           cAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                  >  2:2803716/1‑93 (MQ=255)
                                                                           cAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                  >  1:5684925/1‑93 (MQ=255)
                                                                           cAGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAg                  >  2:2349493/1‑93 (MQ=255)
                                                                            aGATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGa                 <  2:2860378/93‑1 (MQ=255)
                                                                             gATCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt                >  2:456380/1‑93 (MQ=255)
                                                                               tCTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGc              <  1:4833250/93‑1 (MQ=255)
                                                                                cTGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGTCCAAGTAAATTTGATAGTACTGCAGTTGCTGCGTTAGATGCt             <  2:5227769/93‑1 (MQ=255)
                                                                                 tGGAGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCtt            >  1:2559559/1‑93 (MQ=255)
                                                                                  ggAGCA‑TGTGCGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCttt           <  2:126007/93‑1 (MQ=255)
                                                                                    aGCA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTga         <  2:3624930/93‑1 (MQ=255)
                                                                                      cA‑TGTACGGGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaa       <  1:3722435/93‑1 (MQ=255)
                                                                                           tAC‑GGGGGGG‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAAc  <  2:4721770/93‑1 (MQ=255)
                                                                                              |                                                                                           
GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCA‑TGTAC‑GGGGGGA‑GTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGAAAC  >  NC_000913/526204‑526386

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