Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC43.6% 130.0 / 32.6 78G219G (GGT→GGGnanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (40/4);  new base C (0/34);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 5.12e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537964‑4538249
                                                                                                                                                   |                                                                                                                                               
ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                            >  2:134522/1‑151 (MQ=255)
  gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                         >  2:916868/1‑151 (MQ=255)
   cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                        >  2:622402/1‑151 (MQ=255)
   cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                         <  2:139413/150‑1 (MQ=255)
    aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                        <  2:28160/150‑1 (MQ=255)
    aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                       >  1:12801/1‑151 (MQ=255)
        ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTcaca                                                                                                                                   >  1:862930/1‑150 (MQ=255)
         cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                                <  2:905627/151‑1 (MQ=255)
              gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgccg                                                                                                                           >  1:1144075/1‑151 (MQ=255)
                  ctccATCTGGTGATGAGCCGGGAGCAAGCTTATTCCATTTTAATTTATTACCGGTGAAATCAAGTTAAAATACGGCTTCCGTTCGCAATCCTGTTTTGGCTTCGCCTTTAATAAGCCGGGTTTTATCCCCTTTATTCACAACcgccgc                                                                                                                          <  2:1133127/146‑1 (MQ=255)
                      aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGcc                                                                                                                   >  1:5939/1‑151 (MQ=255)
                              atgaGACGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTAAAATCAAGTTAAAATCCGGCATCCTTTCGCAATCCGGTTTTGGTTTCGCATTTAATAAGCCAGGTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTa                                                                                                           <  1:1144530/151‑1 (MQ=255)
                                gagaCGGGAGCAAGTTTTTTCCATTTTAAATTATTCCCGGTAAAATAAATTTAAAATACGGCACCCGTTCGCAACCCTGTTTTGGCTCCCCCATAAATAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                          <  2:834885/150‑1 (MQ=255)
                                         gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                >  2:409035/1‑151 (MQ=255)
                                          cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                               >  2:862070/1‑151 (MQ=255)
                                           aaGCTTATTCTATTTAAAATTATTCCCGGTAAAATAAAGTTAAAATACGGCACCCTTCCGCAATCCTGTTTTGGCTTCGCCATAAAAAGGCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                               <  2:1164386/150‑1 (MQ=255)
                                           aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                              >  2:1029462/1‑151 (MQ=255)
                                             gtttATTCCATTTTAAATTATTCCGGGTGAAATCAAGTTCAAACACGGCATCCGTCCGCAATCCTGTTTGGGCTTCGCCATTAATAAGCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                             <  1:11047/148‑1 (MQ=255)
                                              cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  1:333356/1‑150 (MQ=255)
                                               ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  2:595635/1‑149 (MQ=255)
                                                tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                         >  1:780473/1‑151 (MQ=255)
                                                    ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                       >  1:745628/1‑149 (MQ=255)
                                                     cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                    >  2:921329/1‑151 (MQ=255)
                                                       ttttAAATTATTCCCGGTAAATTCAAGTTAAAACACGGCATCCGTTCGCATTTCTGGTTTGGTTTCCCAATTAATAGCCCAGGTTTTACCCCTTTTATTCAAAACCGCCGCGCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                   <  2:89139/150‑1 (MQ=255)
                                                       tttaAAATTATTACCGGTAAAATAAAGTTAAAATACGGCATCCGTTCGCAACCCGGTTTGGGTTTCGCCATAAATAAGCCAGGTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                      <  2:40384/147‑1 (MQ=255)
                                                         ttAAATTATTACGGGTGAAACAAATTTAAAATCCGCCATCCGTCCGCAATCCTGTTTTGGCTCCGCCATAAAAAAGCCAGTTTTTACCCCCTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                 <  2:436830/150‑1 (MQ=255)
                                                         ttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTCCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                 >  1:755577/1‑150 (MQ=255)
                                                           aaaTTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                              >  2:641321/1‑151 (MQ=255)
                                                                    cgggtaaaacaaaGTTAAAATAGGGCACCCTTTCGCAACCCGGTTTGGGCTTCGCCATAAAAAAGCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                      <  1:231934/139‑1 (MQ=255)
                                                                    ccGGTGAAATCAATTTAAAATACGGCATCCTTTCGCAACCCTGTTTTGGCTTCGCCATAAAAAAGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                       <  1:827551/149‑1 (MQ=255)
                                                                    ccGGTAAAATCAATTTCAATTCCGGCACCCTTTCCCAATCCTGGTTTGGTTCCCCCATTAATAACCCGGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                       <  1:231333/149‑1 (MQ=255)
                                                                        taaaaTCAAGTCAAATTACGGCATCCGTCCTCAATCCTTTTTTTGTTTCGCCTTTAAAAAGCCAGGTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                  <  1:464644/148‑1 (MQ=255)
                                                                          aaaaCAAGTTAAATTACGGCATCCGTTCGCAATCCGGTTTTGGCTTCGCCATTAAAAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                <  1:947794/150‑1 (MQ=255)
                                                                          aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                               >  1:870448/1‑151 (MQ=255)
                                                                           aaTCAAGTTCAAATACGGCATCCGTTCACAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                              >  1:377848/1‑151 (MQ=255)
                                                                              cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                            >  2:1096687/1‑150 (MQ=255)
                                                                                 gTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                           >  1:155079/1‑148 (MQ=255)
                                                                                  ttCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                       >  1:187387/1‑151 (MQ=255)
                                                                                    cAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                      >  2:897122/1‑150 (MQ=255)
                                                                                     aaatcgggCACCGGTTCGCAATCCTGTTTGGGCTTCGCAATTAAAAGGCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                       <  1:1064024/142‑1 (MQ=255)
                                                                                        tACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                    >  1:98620/1‑148 (MQ=255)
                                                                                          cGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                >  2:1018061/1‑150 (MQ=255)
                                                                                           ggCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                              >  1:795050/1‑151 (MQ=255)
                                                                                           ggCACCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                <  1:910769/149‑1 (MQ=255)
                                                                                                ccGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                         >  1:422868/1‑151 (MQ=255)
                                                                                                ccGTTCGAAACCCTGTTTTGGTTTCGCCTTAAAAAAGCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                               <  1:223890/145‑1 (MQ=255)
                                                                                                    tCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       >  2:243179/1‑149 (MQ=255)
                                                                                                      gCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                    >  1:522473/1‑150 (MQ=255)
                                                                                                      cgaaccccggTTTGGGTTCCCCAAAAAAAAACCCAGTTTTTACCCCGTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                <  2:860425/131‑1 (MQ=255)
                                                                                                      cccaTCCCGGTTTGGCCTCCCCTATTAAAAAGCCTGTTTTTATCCCCTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       <  2:264305/145‑1 (MQ=255)
                                                                                                        atccttggtttggttTCCCCATAAAAAACCCGGGTTTTACCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                    <  2:337366/143‑1 (MQ=255)
                                                                                                         acccTGGTTGGGCAACACCTTTAATCGGCCATGTTTTATCCCTTTTATCCATAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 <  1:191228/148‑1 (MQ=255)
                                                                                                          tcttgttttGCCTCCCCCATTAAAAAGCCGTGTTTAGCCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       <  2:351065/137‑1 (MQ=255)
                                                                                                           ccTGTTTTGTTTTCGCCATAAATAACCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                <  1:853167/149‑1 (MQ=255)
                                                                                                           ccTGTTTTGGCTTCCCCATAAAAAAGCCAGTTTTTACCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                              <  2:165160/151‑1 (MQ=255)
                                                                                                           ccTGTTTGGGCTCCCCCTTTAAAAGGCCGGTTTTTACCCCCTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                <  1:354917/149‑1 (MQ=255)
                                                                                                             tggtttggttCCGCAATAAAAAAGCAGGTTTTTACCCCCTTAATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                              <  1:914227/149‑1 (MQ=255)
                                                                                                                tttGGTTTCGCCATTAATAACCCGGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                     <  1:895137/139‑1 (MQ=255)
                                                                                                                tttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                         >  1:84438/1‑151 (MQ=255)
                                                                                                                   ggCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      >  1:440263/1‑151 (MQ=255)
                                                                                                                   ggCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      >  2:466942/1‑151 (MQ=255)
                                                                                                                   ggCTCCGCCATTAAAAACCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                               <  2:969522/142‑1 (MQ=255)
                                                                                                                    gtatcctcattaaataaACCAGGTTTTATCACTTTATTCAAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                        <  1:133335/137‑1 (MQ=255)
                                                                                                                     gttccccctttAATAAGCCAGGTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                      <  2:397389/141‑1 (MQ=255)
                                                                                                                      tcccccATTAAAAACCCAGTTTTTACCCCTTTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                     <  2:1144046/145‑1 (MQ=255)
                                                                                                                       tCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                  >  2:358053/1‑151 (MQ=255)
                                                                                                                        ccccATTAATAAGCCAGGTTTTATCCCCTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                   <  1:1160103/147‑1 (MQ=255)
                                                                                                                        cGCCATAAAAAAGCCGGTTTTTATCCCTTTTTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  2:413100/151‑1 (MQ=255)
                                                                                                                         gCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCAATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                >  2:241912/1‑151 (MQ=255)
                                                                                                                         cgcATTAAAAAAGCCGGTTTTTTCCCCTTTAATTCAAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                 <  1:86164/149‑1 (MQ=255)
                                                                                                                           cATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc              >  1:999616/1‑151 (MQ=255)
                                                                                                                                aataaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa           >  2:700436/1‑150 (MQ=255)
                                                                                                                                 aagaGCCGGGTTTTTTCCCCCTTCTTCCCAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa            <  2:143424/146‑1 (MQ=255)
                                                                                                                                    aGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      >  2:50093/1‑151 (MQ=255)
                                                                                                                                    aGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      >  1:334117/1‑151 (MQ=255)
                                                                                                                                       ccAGGTTTTTCCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg     <  2:78667/150‑1 (MQ=255)
                                                                                                                                        cAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg     >  1:1043660/1‑149 (MQ=255)
                                                                                                                                         aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc  >  2:602612/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537964‑4538249

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: