Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,362,404 G→T 6.6% intergenic (‑193/+27) yfaX ← / ← yfaY putative DNA‑binding transcriptional regulator YfaX/IPR008135 CinA family protein YfaY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,362,4040GT 6.6% 198.9 / 2.5 60intergenic (‑193/+27)yfaX/yfaYputative DNA‑binding transcriptional regulator YfaX/IPR008135 CinA family protein YfaY
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base T (2/2);  total (32/28)
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.42e-01

GGGGAGTAAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGAT  >  NC_000913/2362256‑2362553
                                                                                                                                                    |                                                                                                                                                     
ggggAGTAAAATTGTGAAGCGCATCAAATCCCGCTGCTCGTCCAGCAGACTGGAGAAATTACGCACAATTTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGc                                                                                                                                                      <  2:1385163/150‑1 (MQ=255)
ggggAGTAAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGc                                                                                                                                                      <  1:1317496/150‑1 (MQ=255)
 gggAGTAAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCa                                                                                                                                                     >  2:620362/1‑150 (MQ=255)
   gAGTAAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGt                                                                                                                                                   <  2:244744/150‑1 (MQ=255)
       aaaaTTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAAc                                                                                                                                               <  1:211875/150‑1 (MQ=255)
           ttGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGc                                                                                                                                           <  1:905479/150‑1 (MQ=255)
             gtgAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTTTACTGCACCGGcac                                                                                                                                         <  2:858272/150‑1 (MQ=255)
              tgAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTTTACTGCACCGGcaca                                                                                                                                        <  2:129149/150‑1 (MQ=255)
                 aGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACagag                                                                                                                                     >  1:967304/1‑150 (MQ=255)
                  gcgcATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGt                                                                                                                                    >  1:201219/1‑150 (MQ=255)
                    gcATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTgg                                                                                                                                  >  2:753259/1‑150 (MQ=255)
                     cATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGc                                                                                                                                 >  2:584017/1‑150 (MQ=255)
                      aTCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCATAGAGTGGCg                                                                                                                                >  2:279368/1‑150 (MQ=255)
                      aTCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCg                                                                                                                                >  2:974923/1‑150 (MQ=255)
                         aaaTCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGaa                                                                                                                             >  1:1425452/1‑150 (MQ=255)
                            tCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATca                                                                                                                          <  1:584017/150‑1 (MQ=255)
                              ccGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATca                                                                                                                          >  2:625119/1‑148 (MQ=255)
                              ccGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATca                                                                                                                          <  1:625119/148‑1 (MQ=255)
                               cGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATcagac                                                                                                                       <  2:626308/150‑1 (MQ=255)
                                gCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATcagaca                                                                                                                      >  1:62490/1‑150 (MQ=255)
                                   cctcGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGAt                                                                                                                   >  2:1300233/1‑150 (MQ=255)
                                           agcagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCAt                                                                                                           <  1:1001568/150‑1 (MQ=255)
                                            gcagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATg                                                                                                          >  2:442595/1‑150 (MQ=255)
                                            gcagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATg                                                                                                          <  1:1066989/150‑1 (MQ=255)
                                             cagACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCGGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATgg                                                                                                         <  1:909870/150‑1 (MQ=255)
                                               gACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGAc                                                                                                       >  1:138478/1‑150 (MQ=255)
                                                  tGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTca                                                                                                    >  1:739484/1‑150 (MQ=255)
                                                   ggAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTcaa                                                                                                   <  1:974923/150‑1 (MQ=255)
                                                    gagaAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTcaac                                                                                                  >  1:516329/1‑150 (MQ=255)
                                                        aaTTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAAcc                                                                                              <  2:1425452/150‑1 (MQ=255)
                                                             cGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAAt                                                                                         >  1:794748/1‑150 (MQ=255)
                                                               cacaAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATcc                                                                                       >  1:1159174/1‑150 (MQ=255)
                                                                    aTTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCa                                                                                  <  2:973557/150‑1 (MQ=255)
                                                                    aTTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCa                                                                                  <  2:1255677/150‑1 (MQ=255)
                                                                       tgtgAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGc                                                                               >  2:1107780/1‑150 (MQ=255)
                                                                           aaTCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCAc                                                                           <  1:1308911/150‑1 (MQ=255)
                                                                               aGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGc                                                                       <  1:806932/150‑1 (MQ=255)
                                                                                 cgcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGa                                                                     <  2:739484/150‑1 (MQ=255)
                                                                                 cgcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGa                                                                     <  1:1285761/150‑1 (MQ=255)
                                                                                 cgcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGa                                                                     >  1:1232513/1‑150 (MQ=255)
                                                                                  gcACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGat                                                                    >  2:1352822/1‑150 (MQ=255)
                                                                                           ttttCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGcc                                                           <  2:1318249/150‑1 (MQ=255)
                                                                                                  aTTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAAc                                                    <  2:695594/150‑1 (MQ=255)
                                                                                                    tGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGAACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCa                                                  >  2:726551/1‑150 (MQ=255)
                                                                                                            acaGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgca                                          >  1:1320383/1‑150 (MQ=255)
                                                                                                            acaGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgca                                          <  2:201219/150‑1 (MQ=255)
                                                                                                            acaGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgca                                          <  2:138478/150‑1 (MQ=255)
                                                                                                               ggTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCATAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCgcagca                                       <  1:279368/150‑1 (MQ=255)
                                                                                                                gTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCat                                      >  2:301849/1‑150 (MQ=255)
                                                                                                                     aCAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCa                                 >  1:887009/1‑150 (MQ=255)
                                                                                                                      cAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTATCCAACGGCGCAGCATATTCAg                                >  1:969071/1‑150 (MQ=255)
                                                                                                                         aTTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGc                             <  1:620362/150‑1 (MQ=255)
                                                                                                                          ttttGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGcc                            <  1:1071856/150‑1 (MQ=255)
                                                                                                                               ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                       >  2:1235178/1‑150 (MQ=255)
                                                                                                                               ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                       >  1:626740/1‑150 (MQ=255)
                                                                                                                               ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                       >  2:764900/1‑150 (MQ=255)
                                                                                                                               ccGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCcatcat                       >  2:776300/1‑150 (MQ=255)
                                                                                                                                       actctgtgCCGGTTTACTGCACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGcaca               >  1:1013613/5‑150 (MQ=255)
                                                                                                                                           tgtgCCGGTTTACTGCACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTc           >  1:1065697/1‑150 (MQ=255)
                                                                                                                                             tgCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCt         >  1:1118667/1‑150 (MQ=255)
                                                                                                                                                    gCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGAt  <  2:967304/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGGAGTAAAATTGTGAAGCGCATCAAATCCCGCTCCTCGTCCAGCAGACTGGAGAAATTACGCACAAATTTGTGAATCAGCGCACGGATTTTTTCTTATTGAGTACAACAGGTTGGACAAATTTTGCCGGAGGGGGAATGTGCCGGTGCAGTAAACCGGCACAGAGTGGCGGAATCAGACAGATAAGGTCATGGACTCAACAACCTCAATCCAGCCATGCTCACTGGCGATATCCTGGCCGTTTAACCAACGGCGCAGCATATTCAGTGCCATCATTGCGCACACTTCCTGACGGAT  >  NC_000913/2362256‑2362553

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