Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 394,041 T→G pseudogene (1072/1449 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,0410TG95.9% 474.1 / 7.7 149pseudogene (1072/1449 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (71/70);  total (74/75)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 150 total aligned reads displayed.

TACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAAC  >  NC_000913/393944‑394140
                                                                                                 |                                                                                                   
tACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGAc                                                                                                   <  1:848604/100‑1 (MQ=255)
  cTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACgg                                                                                                 >  1:1046210/1‑100 (MQ=255)
  cTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACgg                                                                                                 >  1:1042676/1‑100 (MQ=255)
  cTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACgg                                                                                                 >  1:3225452/1‑100 (MQ=255)
  cTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACgg                                                                                                 >  2:1581151/1‑100 (MQ=255)
   ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGc                                                                                                <  1:3125632/100‑1 (MQ=255)
   ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGc                                                                                                <  2:2316360/100‑1 (MQ=255)
   ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGc                                                                                                >  1:1639318/1‑100 (MQ=255)
   ttACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGc                                                                                                <  1:1108942/100‑1 (MQ=255)
    tACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCa                                                                                               >  1:2038420/1‑100 (MQ=255)
     aCGTGGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAg                                                                                              <  2:3796235/100‑1 (MQ=255)
      cGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGc                                                                                             >  2:2834219/1‑100 (MQ=255)
      cGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGc                                                                                             <  2:2727438/100‑1 (MQ=255)
      cGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGc                                                                                             <  2:2727420/100‑1 (MQ=255)
      cGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGc                                                                                             >  2:633818/1‑100 (MQ=255)
      cGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGc                                                                                             <  1:4043603/100‑1 (MQ=255)
       gTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCg                                                                                            >  2:2790905/1‑100 (MQ=255)
             ggTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTg                                                                                      <  2:3151526/100‑1 (MQ=255)
             ggTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTg                                                                                      >  1:1663563/1‑100 (MQ=255)
                aGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTc                                                                                   <  2:1611872/100‑1 (MQ=255)
                aGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTc                                                                                   <  1:885751/100‑1 (MQ=255)
                   aCCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGt                                                                                >  2:2549615/1‑100 (MQ=255)
                   aCCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGt                                                                                >  2:3460766/1‑100 (MQ=255)
                    ccAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGtt                                                                               <  2:943657/100‑1 (MQ=255)
                     cAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTc                                                                              <  1:238293/100‑1 (MQ=255)
                     cAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTc                                                                              <  2:3821777/100‑1 (MQ=255)
                      aaCTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTcc                                                                             <  1:3040485/100‑1 (MQ=255)
                        cTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGGGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAg                                                                           >  1:2369153/1‑100 (MQ=255)
                        cTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTCCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAg                                                                           >  2:1483545/1‑100 (MQ=255)
                          ccGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTc                                                                         >  2:292593/1‑100 (MQ=255)
                          ccGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTc                                                                         >  1:3372361/1‑100 (MQ=255)
                            gACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTg                                                                       <  1:1549718/100‑1 (MQ=255)
                             aCGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTCCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTgg                                                                      <  1:1509590/100‑1 (MQ=255)
                              cGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGt                                                                     <  2:3175688/100‑1 (MQ=255)
                               gCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtg                                                                    <  2:3038386/100‑1 (MQ=255)
                               gCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtg                                                                    <  2:1262955/100‑1 (MQ=255)
                                  tGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGtgatg                                                                 >  1:293637/1‑100 (MQ=255)
                                    gggCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTATCTGGTCTGTTCCAGTCTGGTGATGAc                                                               >  2:1812330/1‑100 (MQ=255)
                                     ggCTTCGGTTTGAAAGCCGGTTACGACTTCAGACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACt                                                              >  2:1240707/1‑100 (MQ=255)
                                     ggCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACt                                                              >  2:661416/1‑100 (MQ=255)
                                      gCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTa                                                             >  1:63300/1‑100 (MQ=255)
                                      gCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTa                                                             >  1:3144158/1‑100 (MQ=255)
                                        ttCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTAcc                                                           <  1:437490/100‑1 (MQ=255)
                                        ttCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTAcc                                                           >  1:2003035/1‑100 (MQ=255)
                                           ggTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGc                                                        >  1:3317930/1‑100 (MQ=255)
                                           ggTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGc                                                        >  1:1612516/1‑100 (MQ=255)
                                           ggTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGc                                                        >  1:3609789/1‑100 (MQ=255)
                                           ggTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGc                                                        <  2:1250219/100‑1 (MQ=255)
                                           ggTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGc                                                        <  1:2304073/100‑1 (MQ=255)
                                            gTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGcg                                                       >  2:2109952/1‑99 (MQ=255)
                                             tttGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTg                                                      <  1:2613766/100‑1 (MQ=255)
                                               tGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAg                                                    <  1:723392/100‑1 (MQ=255)
                                               tGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAg                                                    <  1:724053/100‑1 (MQ=255)
                                                gAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGc                                                   >  2:2348191/1‑100 (MQ=255)
                                                 aaaGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCa                                                  >  2:1059406/1‑100 (MQ=255)
                                                    gCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACg                                               <  2:2482304/100‑1 (MQ=255)
                                                    gCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACg                                               <  2:92351/100‑1 (MQ=255)
                                                     ccGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGa                                              >  1:4085536/1‑100 (MQ=255)
                                                      cGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGAc                                             <  1:3982887/100‑1 (MQ=255)
                                                      cGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGAc                                             <  2:1756861/100‑1 (MQ=255)
                                                      cGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGAc                                             <  2:2291763/100‑1 (MQ=255)
                                                            ngACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGaaa                                       >  1:3253399/2‑100 (MQ=255)
                                                            cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGCCGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGCCTACCAGCTGAGCAACGACATGaaa                                       >  2:788887/1‑100 (MQ=255)
                                                            cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGaaa                                       >  1:2126951/1‑100 (MQ=255)
                                                            cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGaaa                                       >  1:3084513/1‑100 (MQ=255)
                                                             gACTGCAACGTGGGTGATGCTGGTTACGTGACTCCTCACGGCAGCGTTGCTGGTCTGTTCCAGTCTCGTGATGACTACCAGCTTGGCAACGACATGAAAg                                      <  2:360185/100‑1 (MQ=255)
                                                              aCTTCAAACTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGt                                     >  2:2089283/1‑100 (MQ=255)
                                                                   aaaCTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACg                                <  1:1187650/100‑1 (MQ=255)
                                                                    aaCTGGGTGATGCTGGTTACGTGACTTCTGACGGCAGCGCTTCTGCTCTGTTCCAGTCTGGTGATGACTATCAGCTGAGCAACGACATGAAAGTTGACgg                               <  1:3811224/100‑1 (MQ=255)
                                                                    aaCTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACgg                               >  2:312437/1‑100 (MQ=255)
                                                                    aaCTGGGTGATGCTGGATACGTGACTCCTGACGGCAGCGTTTCTGTTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACgg                               <  2:305249/100‑1 (MQ=255)
                                                                      cTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACggcc                             <  2:793051/100‑3 (MQ=255)
                                                                      cTGGGTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTc                             <  1:2109282/100‑1 (MQ=255)
                                                                      cTGGGGGATGCTGGTGACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTc                             <  2:3053630/100‑1 (MQ=255)
                                                                        gggTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGACAGTTGACGGTCAg                           >  1:3676022/1‑100 (MQ=255)
                                                                        gggTGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAg                           >  2:1539259/1‑100 (MQ=255)
                                                                           tGATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCt                        >  2:556965/1‑100 (MQ=255)
                                                                            gATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCtt                       >  1:2026508/1‑100 (MQ=255)
                                                                            gATGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCtt                       >  2:1173065/1‑100 (MQ=255)
                                                                             aTGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTa                      <  1:537942/100‑1 (MQ=255)
                                                                             aTGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTa                      <  1:982069/100‑1 (MQ=255)
                                                                              tGCTGGTTACGTGACTCCTGACGTCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTAc                     <  2:1074096/100‑1 (MQ=255)
                                                                              tGCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTAc                     >  2:2758078/1‑100 (MQ=255)
                                                                               gCTGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACg                    <  1:4009004/100‑1 (MQ=255)
                                                                                cTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGa                   <  2:3263214/100‑1 (MQ=255)
                                                                                 tGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGAc                  <  2:1021034/100‑1 (MQ=255)
                                                                                 tGGTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGAc                  <  1:2350700/100‑1 (MQ=255)
                                                                                  ggTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACa                 >  2:3663376/1‑100 (MQ=255)
                                                                                   gTTACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAg                <  1:1310440/100‑1 (MQ=255)
                                                                                    ttACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGc               <  1:1081370/100‑1 (MQ=255)
                                                                                    ttACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGc               <  2:2269775/100‑1 (MQ=255)
                                                                                     tACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCa              >  2:1731238/1‑100 (MQ=255)
                                                                                     tACGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCa              >  2:1730055/1‑100 (MQ=255)
                                                                                       cGTGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATg            >  2:505315/1‑100 (MQ=255)
                                                                                         tGACTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCg          >  1:194141/1‑100 (MQ=255)
                                                                                            cTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAGGTTGACGGTCAGTCTTACGACAGCATGCGTTa       <  2:1501075/100‑1 (MQ=255)
                                                                                            cTCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTa       >  1:2629218/1‑100 (MQ=255)
                                                                                             tCCTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTAt      <  1:2561777/100‑1 (MQ=255)
                                                                                              ccTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATg     >  1:1830327/1‑100 (MQ=255)
                                                                                              ccTGACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATg     >  2:1658903/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
TACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAAC  >  NC_000913/393944‑394140

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