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)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCGCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCtttt                                                                                    >  1:1427662/1‑100 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCtttt                                                                                    <  2:481514/100‑1 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCtttt                                                                                    <  2:1493034/100‑1 (MQ=255)
            cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCtttt                                                                                    >  1:318801/1‑100 (MQ=255)
             aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCttttt                                                                                   <  2:2136640/100‑1 (MQ=255)
             aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCttttt                                                                                   <  2:1053728/100‑1 (MQ=255)
               cGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCa                                                                                 <  2:2478240/100‑1 (MQ=255)
                gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCaa                                                                                <  1:578368/100‑1 (MQ=255)
                gCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCaa                                                                                >  1:267902/1‑100 (MQ=255)
                 ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAAc                                                                               <  2:2555961/100‑1 (MQ=255)
                 acgctcttCCGATCT‑CCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACg                                                                              <  1:2295258/94‑1 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGAt                                                                            <  2:655434/100‑1 (MQ=255)
                    gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGAt                                                                            <  2:205029/100‑1 (MQ=255)
                      tttCTGATCTACCGCAAATGCAAATCCCGCTTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTa                                                                          <  1:659152/100‑1 (MQ=255)
                       ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTa                                                                          <  1:1278681/100‑1 (MQ=255)
                       ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTaa                                                                         >  2:396583/1‑100 (MQ=255)
                         cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgc                                                                       >  1:852554/1‑100 (MQ=255)
                         cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgc                                                                       >  2:335289/1‑100 (MQ=255)
                           gATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                     >  2:2192576/1‑100 (MQ=255)
                           gATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCGCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAgcgc                                                                     >  1:1431889/1‑100 (MQ=255)
                            aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcc                                                                    >  1:1033409/1‑100 (MQ=255)
                            aTCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcc                                                                    >  2:968516/1‑100 (MQ=255)
                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                 >  2:1175626/1‑100 (MQ=255)
                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                 >  1:130957/1‑100 (MQ=255)
                               tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGcgcccgc                                                                 >  2:1541038/1‑100 (MQ=255)
                                 ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGcgcccgcc                                                                <  2:2521343/100‑1 (MQ=255)
                                  cGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                              >  2:1741300/1‑100 (MQ=255)
                                   gcaaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCttt                                                             >  1:2334643/1‑100 (MQ=255)
                                    caaatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGAC‑GCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTa                                                            <  1:2398208/100‑1 (MQ=255)
                                      aatgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCGCCCATCTTTTTCAACGATTAACCGCACGCCTTTAAt                                                          >  1:929989/1‑100 (MQ=255)
                                        tgcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                        >  2:107869/1‑100 (MQ=255)
                                         gcaaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGt                                                       <  2:235029/100‑1 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                      <  2:2349548/100‑1 (MQ=255)
                                          caaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTg                                                      <  2:2368784/100‑1 (MQ=255)
                                           aaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                     <  2:30661/100‑1 (MQ=255)
                                           aaatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                     <  1:1217629/100‑1 (MQ=255)
                                            aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTgg                                                     <  1:232950/100‑1 (MQ=255)
                                               ccccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTAt                                                  <  1:991912/100‑1 (MQ=255)
                                               ccccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATc                                                 <  1:2419322/100‑1 (MQ=255)
                                                cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCa                                                >  1:1636525/1‑100 (MQ=255)
                                                cccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCa                                                >  2:1143812/1‑100 (MQ=255)
                                                 ccccTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCAt                                               >  2:1819375/1‑100 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCGTTAATCGTGGTATCATTCGc                                           <  1:773535/100‑1 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                          >  2:1850334/1‑100 (MQ=255)
                                                      ctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCt                                          <  2:1105277/100‑1 (MQ=255)
                                                         tAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcg                                       <  2:1076777/100‑1 (MQ=255)
                                                          aaCCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcg                                       >  2:1574453/1‑100 (MQ=255)
                                                             cTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcgcaa                                    <  1:1940913/100‑3 (MQ=255)
                                                              tCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                  >  2:1128017/1‑100 (MQ=255)
                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                  >  2:1759834/1‑100 (MQ=255)
                                                               ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                 >  1:209517/1‑100 (MQ=255)
                                                                cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaaa                                >  2:1696295/1‑100 (MQ=255)
                                                                 agtagtCCTCCCTTTTCGATAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaaa                                >  2:1389146/1‑100 (MQ=255)
                                                                 agtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAc                               >  2:1205409/1‑100 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                              <  2:1903911/100‑1 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                              <  2:169009/100‑1 (MQ=255)
                                                                  gtagtCCTCCCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAAcc                              <  1:2237818/100‑1 (MQ=255)
                                                                      tcctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                          <  2:130958/100‑1 (MQ=255)
                                                                       cctccCTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGGGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                         <  1:1438502/100‑1 (MQ=255)
                                                                          cccTTTTCTATAGCGACAAAGAC‑GCCCCTCCATCTGTTTCAACGATTAAGCGCCTGCCTTTATTCGTGTTATTATTCTCTGCGCCATAAAcc                              <  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)
                                                                            cTTTTCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCacac                    >  1:1112348/1‑100 (MQ=255)
                                                                                tCGATAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCat                 >  1:896171/1‑99 (MQ=255)
                                                                                   aTAGCGACAAAGACGGCCCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCattat              >  2:1796927/1‑100 (MQ=255)
                                                                                    tAGCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTg            <  2:1636534/100‑1 (MQ=255)
                                                                                      gCGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAg          <  2:1675966/100‑1 (MQ=255)
                                                                                       cGACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGc         <  2:2235822/100‑1 (MQ=255)
                                                                                        gACAAAGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCg        <  2:1332938/100‑1 (MQ=255)
                                                                                             aGACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTc   >  2:2280555/1‑100 (MQ=255)
                                                                                              gACGG‑CCCCCCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCt  >  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: ‑