Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,288,561 (C)7→6 100% 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.78.1% 181.1 / 66.5 73pseudogene (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 (8/8);  new base . (29/28);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCT  >  NC_000913/2288456‑2288650
                                                                                                   |                                                                                               
ttgttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGc                                                                                                 <  1:1962910/100‑1 (MQ=255)
 tgttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcc                                                                                                >  2:1855205/1‑100 (MQ=255)
  gttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                               >  2:2214727/1‑100 (MQ=255)
  gttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                               >  2:768403/1‑100 (MQ=255)
  gttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccc                                                                                               >  2:1940923/1‑100 (MQ=255)
   ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccc                                                                                              >  2:41448/1‑100 (MQ=255)
   ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccc                                                                                              >  2:415504/1‑100 (MQ=255)
   ttttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCCTTAACCTCACGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccc                                                                                              <  1:1600681/100‑1 (MQ=255)
    tttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccccc                                                                                             <  2:1923979/100‑1 (MQ=255)
    tttAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGccccc                                                                                             <  2:955671/100‑1 (MQ=255)
      tAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGcccccca                                                                                           >  2:231723/1‑99 (MQ=255)
          aGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCt                                                                                       >  2:1278686/1‑100 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCGCTTTTCGATAGCGACAAAGACGGCCCCCCATCtttt                                                                                    >  1:1427662/1‑100 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCtttt                                                                                    <  2:481514/100‑1 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCtttt                                                                                    <  2:1493034/100‑1 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCtttt                                                                                    >  1:318801/1‑100 (MQ=255)
             aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCttttt                                                                                   <  2:2136640/100‑1 (MQ=255)
             aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCttttt                                                                                   <  2:1053728/100‑1 (MQ=255)
               cGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCa                                                                                 <  2:2478240/100‑1 (MQ=255)
                gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCaa                                                                                <  1:578368/100‑1 (MQ=255)
                gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCaa                                                                                >  1:267902/1‑100 (MQ=255)
                 ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAAc                                                                               <  2:2555961/100‑1 (MQ=255)
                 acgctcttCCGATCTCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACg                                                                              <  1:2295258/94‑1 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAt                                                                            <  2:655434/100‑1 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGAt                                                                            <  2:205029/100‑1 (MQ=255)
                      tttCTGATCTACCGCAAATGCAAATCCCGCTTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTa                                                                          <  1:659152/100‑1 (MQ=255)
                       ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTa                                                                          <  1:1278681/100‑1 (MQ=255)
                       ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTaa                                                                         >  2:396583/1‑100 (MQ=255)
                         cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgc                                                                       >  1:852554/1‑100 (MQ=255)
                         cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgc                                                                       >  2:335289/1‑100 (MQ=255)
                           gATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgcgc                                                                     >  2:2192576/1‑100 (MQ=255)
                           gATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCGCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAgcgc                                                                     >  1:1431889/1‑100 (MQ=255)
                            aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                    >  1:1033409/1‑100 (MQ=255)
                            aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcc                                                                    >  2:968516/1‑100 (MQ=255)
                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                 >  2:1175626/1‑100 (MQ=255)
                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                 >  1:130957/1‑100 (MQ=255)
                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                 >  2:1541038/1‑100 (MQ=255)
                                 ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGcgcccgcc                                                                <  2:2521343/100‑1 (MQ=255)
                                  cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                              >  2:1741300/1‑100 (MQ=255)
                                   gcaaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                             >  1:2334643/1‑100 (MQ=255)
                                    caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTa                                                            <  1:2398208/100‑1 (MQ=255)
                                      aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCGCCCATCTTTTTCAACGATTAACCGCACGCCTTTAAt                                                          >  1:929989/1‑100 (MQ=255)
                                        tgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                        >  2:107869/1‑100 (MQ=255)
                                         gcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGt                                                       <  2:235029/100‑1 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                      <  2:2349548/100‑1 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                      <  2:2368784/100‑1 (MQ=255)
                                           aaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                     <  2:30661/100‑1 (MQ=255)
                                           aaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                     <  1:1217629/100‑1 (MQ=255)
                                            aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                     <  1:232950/100‑1 (MQ=255)
                                               ccccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTAt                                                  <  1:991912/100‑1 (MQ=255)
                                               ccccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATc                                                 <  1:2419322/100‑1 (MQ=255)
                                                cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCa                                                >  1:1636525/1‑100 (MQ=255)
                                                cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCa                                                >  2:1143812/1‑100 (MQ=255)
                                                 ccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCAt                                               >  2:1819375/1‑100 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCGTTAATCGTGGTATCATTCGc                                           <  1:773535/100‑1 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                          >  2:1850334/1‑100 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                          <  2:1105277/100‑1 (MQ=255)
                                                         tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcg                                       <  2:1076777/100‑1 (MQ=255)
                                                          aaCCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcg                                       >  2:1574453/1‑100 (MQ=255)
                                                             cTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcgcaa                                    <  1:1940913/100‑3 (MQ=255)
                                                              tCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                  >  2:1128017/1‑100 (MQ=255)
                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                  >  2:1759834/1‑100 (MQ=255)
                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                 >  1:209517/1‑100 (MQ=255)
                                                                cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaaa                                >  2:1696295/1‑100 (MQ=255)
                                                                 agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaaa                                >  2:1389146/1‑100 (MQ=255)
                                                                 agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                               >  2:1205409/1‑100 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                              <  2:1903911/100‑1 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                              <  2:169009/100‑1 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                              <  1:2237818/100‑1 (MQ=255)
                                                                      tcctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                          <  2:130958/100‑1 (MQ=255)
                                                                       cctccCTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGGGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                         <  1:1438502/100‑1 (MQ=255)
                                                                          cccTTTTCTATAGCGACAAAGACGCCCCTCCATCTGTTTCAACGATTAAGCGCCTGCCTTTATTCGTGTTATTATTCTCTGCGCCATAAAcc                              <  1:818096/92‑1 (MQ=255)
                                                                           ccTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcac                      >  2:1192175/1‑100 (MQ=255)
                                                                           ccTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcac                      >  1:491707/1‑100 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                     >  2:295392/1‑100 (MQ=255)
                                                                            cTTTTCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                    >  1:1112348/1‑100 (MQ=255)
                                                                                tCGATAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCat                 >  1:896171/1‑99 (MQ=255)
                                                                                   aTAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattat              >  2:1796927/1‑100 (MQ=255)
                                                                                    tAGCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTg            <  2:1636534/100‑1 (MQ=255)
                                                                                      gCGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAg          <  2:1675966/100‑1 (MQ=255)
                                                                                       cGACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGc         <  2:2235822/100‑1 (MQ=255)
                                                                                        gACAAAGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCg        <  2:1332938/100‑1 (MQ=255)
                                                                                             aGACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTc   >  2:2280555/1‑100 (MQ=255)
                                                                                              gACGGCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCt  >  2:298392/1‑100 (MQ=255)
                                                                                                   |                                                                                               
TTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCT  >  NC_000913/2288456‑2288650

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