Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,995,873 C→A pseudogene (148/707 nt) ygeN ← putative type III secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,995,8730CA96.9% 217.5 / 0.5 66pseudogene (148/707 nt)ygeNputative type III secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (48/15);  total (51/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTATCATGTTTATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTC  >  NC_000913/2995739‑2996000
                                                                                                                                      |                                                                                                                               
ttATCATGTTTATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCaa                                                                                                                            >  1:537247/1‑140 (MQ=255)
    cATGTTTATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATAct                                                                                                                        <  2:463691/140‑1 (MQ=255)
    cATGTTTATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATAct                                                                                                                        <  2:463692/140‑1 (MQ=255)
     aTGTTTATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATActc                                                                                                                       >  1:533836/1‑140 (MQ=255)
        tttATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGa                                                                                                                    >  2:298348/1‑140 (MQ=255)
                 ttttttCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTGCTCAATACTCTGAATTTTTGCg                                                                                                           >  2:851000/1‑140 (MQ=255)
                 ttttttCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCg                                                                                                           >  1:264867/1‑140 (MQ=255)
                      tCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCggtggt                                                                                                      <  2:557036/140‑1 (MQ=255)
                             cTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCaaga                                                                                               >  2:60971/1‑140 (MQ=255)
                                     aGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGAtttt                                                                                       >  1:521745/1‑140 (MQ=255)
                                       cAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGAtttttt                                                                                     >  1:207044/1‑140 (MQ=255)
                                               tCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAaga                                                                             >  2:611866/1‑140 (MQ=255)
                                                    ggTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAAc                                                                        >  1:592983/1‑140 (MQ=255)
                                                      tataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTa                                                                      >  2:93835/1‑140 (MQ=255)
                                                      tataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTa                                                                      >  2:576041/1‑140 (MQ=255)
                                                      tataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTa                                                                      >  2:687799/1‑140 (MQ=255)
                                                      tataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTa                                                                      >  1:483750/1‑140 (MQ=255)
                                                      tataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTa                                                                      >  1:296721/1‑140 (MQ=255)
                                                       ataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAt                                                                     <  1:901638/140‑1 (MQ=255)
                                                       ataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAt                                                                     <  2:470359/140‑1 (MQ=255)
                                                       ataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAt                                                                     <  1:298348/140‑1 (MQ=255)
                                                              tATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAtttctttc                                                              >  1:545100/1‑140 (MQ=255)
                                                               aTCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCa                                                             >  1:287696/1‑140 (MQ=255)
                                                               aTCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCa                                                             >  1:750036/1‑140 (MQ=255)
                                                                     ctcATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCAt                                                       >  2:733679/1‑140 (MQ=255)
                                                                         tcaATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGa                                                   <  2:783372/140‑1 (MQ=255)
                                                                         tcaATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGa                                                   <  2:592983/140‑1 (MQ=255)
                                                                               atGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAAtttt                                             >  2:212632/1‑140 (MQ=255)
                                                                                tGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAAttttt                                            >  2:833554/1‑140 (MQ=255)
                                                                                  ccTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCt                                          >  1:336968/1‑140 (MQ=255)
                                                                                  ccTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCt                                          >  2:804054/1‑140 (MQ=255)
                                                                                  ccTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCt                                          >  2:712385/1‑140 (MQ=255)
                                                                                  ccTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCt                                          >  2:866358/1‑140 (MQ=255)
                                                                                  ccTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCt                                          >  2:263525/1‑140 (MQ=255)
                                                                                    ttttGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCttt                                        >  1:383959/1‑140 (MQ=255)
                                                                                       tGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTaa                                     >  2:260472/1‑140 (MQ=255)
                                                                                       tGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTaa                                     >  2:279774/1‑140 (MQ=255)
                                                                                         tataGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATc                                   >  1:872296/1‑140 (MQ=255)
                                                                                           taGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAATTCTTTCAATCCATCAGAAATTTTTCTTTTAATcac                                 >  1:183304/1‑140 (MQ=255)
                                                                                            aGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATcaca                                >  2:41473/1‑140 (MQ=255)
                                                                                              ccATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCacaac                              <  1:212562/140‑1 (MQ=255)
                                                                                               catcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCacaaca                             <  1:708378/140‑1 (MQ=255)
                                                                                                atcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCTATCCATCAGAAATTTTTCTTTTAATCACAacac                            >  2:888983/1‑140 (MQ=255)
                                                                                                 tcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCTATCCATCAGAAATTTTTCTTTTAATCACAACAcc                           <  1:888983/140‑1 (MQ=255)
                                                                                                 tcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACAcc                           >  2:530805/1‑140 (MQ=255)
                                                                                                   ataataACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAt                         >  2:753763/1‑140 (MQ=255)
                                                                                                   ataataACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAt                         >  2:492080/1‑140 (MQ=255)
                                                                                                   ataataACCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAt                         >  2:672010/1‑140 (MQ=255)
                                                                                                   ataataACATTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAt                         >  2:735607/1‑140 (MQ=255)
                                                                                                        aaCCTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCa                    <  1:804054/140‑1 (MQ=255)
                                                                                                          ccTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAgg                  >  2:307184/1‑140 (MQ=255)
                                                                                                          ccTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAgg                  >  1:723261/1‑140 (MQ=255)
                                                                                                           cTTGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGga                 >  1:452007/1‑140 (MQ=255)
                                                                                                             tGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGgaga               >  2:368149/1‑140 (MQ=255)
                                                                                                             tGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGgaga               >  2:133379/1‑140 (MQ=255)
                                                                                                             tGAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGgaga               >  1:356941/1‑140 (MQ=255)
                                                                                                              gAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAt              >  1:806697/1‑140 (MQ=255)
                                                                                                              gAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAt              >  2:57695/1‑140 (MQ=255)
                                                                                                              gAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAt              <  1:576041/140‑1 (MQ=255)
                                                                                                              gAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAt              >  1:293471/1‑140 (MQ=255)
                                                                                                              gAATAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAt              >  2:40414/1‑140 (MQ=255)
                                                                                                               aaTAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAtt             >  1:85278/1‑140 (MQ=255)
                                                                                                               aaTAAATTTTTCTTCTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAtt             >  1:128767/1‑140 (MQ=255)
                                                                                                                        ttcttcTCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGc    <  2:522929/140‑1 (MQ=255)
                                                                                                                         tcttctCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCt   >  1:841227/1‑140 (MQ=255)
                                                                                                                          cttctCTCTTTTACTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTc  <  1:263525/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                               
TTATCATGTTTATTCTTTTTTTTCAAGTTCTGAGCATAGCAAGCTTATCAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTC  >  NC_000913/2995739‑2996000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: