Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,288,561 (C)7→6 pseudogene (354/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,288,5550C.99.2% 543.0 / ‑2.7 122pseudogene (360/2525 nt)yejOpseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (76/45);  total (76/46)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 134 total aligned reads displayed.

TTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTA  >  NC_000913/2288456‑2288646
                                                                                                   |                                                                                           
 tgttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGc                                                                                             <  2:390485/99‑1 (MQ=255)
  gttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                           <  2:334736/100‑1 (MQ=255)
  gttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcag                                                                                           >  1:390486/1‑98 (MQ=255)
   ttttAATAGCACCTCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                           >  2:287897/1‑99 (MQ=255)
   tcttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAACTCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccc                                                                                          >  1:1956142/3‑100 (MQ=255)
    tttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCTGTAGTCCTCCCTTTTCGATAGCGATAAAGACGGccccc                                                                                         <  2:1784795/100‑1 (MQ=255)
    tttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccccc                                                                                         >  1:2586731/1‑100 (MQ=255)
      tAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccccca                                                                                       >  1:1650939/1‑99 (MQ=255)
        aTAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccccccatc                                                                                     <  1:273288/100‑4 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCtttt                                                                                >  2:2325642/1‑100 (MQ=255)
              ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAGCCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTc                                                                              <  2:652016/100‑1 (MQ=255)
              ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTc                                                                              >  2:1202203/1‑100 (MQ=255)
               cGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCa                                                                             <  2:735381/100‑1 (MQ=255)
                 ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAAc                                                                           <  1:388146/100‑1 (MQ=255)
                  cTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACg                                                                          <  2:812050/100‑1 (MQ=255)
                  cTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACg                                                                          <  1:2753265/100‑1 (MQ=255)
                  cTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACg                                                                          <  2:227430/100‑1 (MQ=255)
                   tGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCtttttt                                                                              >  2:2378132/1‑94 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAt                                                                        <  2:1776231/100‑1 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAt                                                                        >  2:976539/1‑100 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAt                                                                        >  1:1270958/1‑100 (MQ=255)
                     cTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAtt                                                                       <  1:2066144/100‑1 (MQ=255)
                     cTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAtt                                                                       <  2:2539124/100‑1 (MQ=255)
                      tttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTa                                                                      <  1:2345264/100‑1 (MQ=255)
                         cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgc                                                                   <  1:1243362/100‑1 (MQ=255)
                         cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgc                                                                   >  1:322171/1‑100 (MQ=255)
                          tGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgcg                                                                  >  2:905368/1‑100 (MQ=255)
                            aTCTACCGCAAGTGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                >  2:2220163/1‑100 (MQ=255)
                             tCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgccc                                                               <  2:2128575/100‑1 (MQ=255)
                              cTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGCTAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAtt                                                                       >  2:2351727/1‑91 (MQ=255)
                                aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcccgcc                                                            <  2:2386139/100‑1 (MQ=255)
                                aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcacccgcc                                                            <  2:554132/100‑1 (MQ=255)
                                 ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCt                                                           >  2:624301/1‑100 (MQ=255)
                                   gcaaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                         >  2:1398156/1‑100 (MQ=255)
                                   gcaaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                         >  1:1412886/1‑100 (MQ=255)
                                    caaatgcaagtCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                         <  1:2297980/99‑1 (MQ=255)
                                     aaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTaa                                                       >  2:1875765/1‑100 (MQ=255)
                                      aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTa                                                        >  2:2877778/1‑98 (MQ=255)
                                      aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                      >  2:2212274/1‑100 (MQ=255)
                                      aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                      >  1:746439/1‑100 (MQ=255)
                                      aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                      >  2:1273697/1‑100 (MQ=255)
                                       atgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                      <  1:1837218/100‑1 (MQ=255)
                                        tgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                    >  1:1030462/1‑100 (MQ=255)
                                        tgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                    >  1:2124484/1‑100 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                  >  2:2242594/1‑100 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                  >  1:2406171/1‑100 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                  >  2:396963/1‑100 (MQ=255)
                                           aaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                 <  2:2759581/100‑1 (MQ=255)
                                           aaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                 <  2:1867929/100‑1 (MQ=255)
                                            aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGt                                                >  1:659904/1‑100 (MQ=255)
                                            aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGt                                                >  1:958575/1‑100 (MQ=255)
                                            aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGt                                                >  1:659808/1‑100 (MQ=255)
                                             atCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTa                                               >  2:1539300/1‑100 (MQ=255)
                                                cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCa                                            >  1:877495/1‑100 (MQ=255)
                                                 ccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCAt                                           >  1:1555286/1‑100 (MQ=255)
                                                    ctctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCg                                        >  2:919417/1‑100 (MQ=255)
                                                    ctctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCg                                        >  2:122825/1‑100 (MQ=255)
                                                    ctctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCg                                        >  1:380394/1‑100 (MQ=255)
                                                     tctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGc                                       <  1:2320965/100‑1 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                      >  1:839573/1‑100 (MQ=255)
                                                       tttAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                     <  1:2139367/100‑1 (MQ=255)
                                                           aCCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGcc                                 >  2:2654712/1‑100 (MQ=255)
                                                            ccTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATGTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                <  2:958573/100‑1 (MQ=255)
                                                             cTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                               >  1:789709/1‑100 (MQ=255)
                                                              tCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                              <  2:2832132/100‑1 (MQ=255)
                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                             <  2:2327775/100‑1 (MQ=255)
                                                                 agtagtCCTCCCTTTTCGATAGCGACAAAGATGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                           >  1:2785266/1‑100 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGTTATCATTCGCTGCTCCATAAAcc                          >  2:875018/1‑100 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                          >  1:1989880/1‑100 (MQ=255)
                                                                     gtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGt                       >  1:1499526/1‑100 (MQ=255)
                                                                     gtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGt                       >  1:1506083/1‑100 (MQ=255)
                                                                     gtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGt                       >  1:2525546/1‑100 (MQ=255)
                                                                       cctccCTTTTTGATATCTACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCAGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                      >  2:2032125/1‑99 (MQ=255)
                                                                        ctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGc                    >  1:1020961/1‑100 (MQ=255)
                                                                        ctaccTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGc                    <  2:2635254/97‑1 (MQ=255)
                                                                         tccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGca                   <  2:2124471/100‑1 (MQ=255)
                                                                           ccTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                 <  1:2418334/100‑1 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                <  2:2809809/100‑1 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                <  1:764298/100‑1 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                >  1:1720836/1‑100 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                >  2:2423973/1‑100 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                >  1:110488/1‑100 (MQ=255)
                                                                             ttttCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACAcc               <  2:1701253/100‑1 (MQ=255)
                                                                               ttCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCat             >  1:2201020/1‑100 (MQ=255)
                                                                                 cGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCatta           >  1:2411508/1‑100 (MQ=255)
                                                                                  gATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattat          >  2:2733716/1‑100 (MQ=255)
                                                                                  gATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattat          <  1:1781300/100‑1 (MQ=255)
                                                                                  gATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattat          >  1:1640852/1‑100 (MQ=255)
                                                                                  gATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattat          <  2:789129/100‑1 (MQ=255)
                                                                                   aTAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattatt         >  2:733480/1‑100 (MQ=255)
                                                                                   aTAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattatt         >  2:907741/1‑100 (MQ=255)
                                                                                   aTAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattatt         <  1:777710/100‑1 (MQ=255)
                                                                                   aTAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCATGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattatt         >  2:913832/1‑100 (MQ=255)
                                                                                     aGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAGGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTg        <  2:1574848/99‑1 (MQ=255)
                                                                                      gCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAg      >  1:2004282/1‑100 (MQ=255)
                                                                                      gCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAg      >  2:1708463/1‑100 (MQ=255)
                                                                                      gCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAg      >  1:965452/1‑100 (MQ=255)
                                                                                        gACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCg    <  2:938499/100‑1 (MQ=255)
                                                                                         aCAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGt   <  1:52141/100‑1 (MQ=255)
                                                                                          cAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTa  >  2:2775906/1‑100 (MQ=255)
                                                                                                   |                                                                                           
TTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTA  >  NC_000913/2288456‑2288646

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: