Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 188,696 A→G 7.8% intergenic (‑46/+16) glnD ← / ← map PII uridylyltransferase/uridylyl removing enzyme/methionine aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913188,6960AG 7.8% 185.3 / 5.6 64intergenic (‑46/+16)glnD/mapPII uridylyltransferase/uridylyl removing enzyme/methionine aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CCGACGGTTAATTCATCACGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCC  >  NC_000913/188550‑188845
                                                                                                                                                  |                                                                                                                                                     
ccGACGGTTAATTCATCACGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTc                                                                                                                                                    >  2:1036295/1‑150 (MQ=255)
   aCGGTTAATTCATCACGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGTGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGTCATTCGcc                                                                                                                                                 <  1:1185187/150‑1 (MQ=255)
       ttAATTCATCACGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCt                                                                                                                                             <  2:197106/150‑1 (MQ=255)
        tAATTCATCACGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCGAATGCCGGCtt                                                                                                                                            <  2:1382220/150‑1 (MQ=255)
            tcatcaCGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATc                                                                                                                                        <  1:355435/150‑1 (MQ=255)
               tcaCGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCGAATGCCGGCttttttaat                                                                                                                                     <  2:97773/150‑8 (MQ=255)
               tcaCGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTa                                                                                                                                     <  1:178100/150‑1 (MQ=255)
                caCGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTAt                                                                                                                                    <  2:1387211/150‑1 (MQ=255)
                caCGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTAt                                                                                                                                    <  2:824822/150‑1 (MQ=255)
                 aCGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTAtt                                                                                                                                   <  2:136321/150‑1 (MQ=255)
                   ggggCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATT‑AAAAAGCCGGCATTCGCCGGCTTCATCTTATtcgt                                                                                                                                >  2:626108/1‑150 (MQ=255)
                    gggCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATtcgt                                                                                                                                >  2:796268/1‑150 (MQ=255)
                      gCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATtcgtcg                                                                                                                              >  2:495507/1‑150 (MQ=255)
                          gACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCg                                                                                                                          >  1:320046/1‑150 (MQ=255)
                             gCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCgaga                                                                                                                       <  2:1113986/150‑1 (MQ=255)
                              cATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCGAATGCCGGCttttt                                                                                                                                         >  1:54371/1‑128 (MQ=255)
                              cATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCGAATGCCGGCttttt                                                                                                                                         <  2:54371/131‑4 (MQ=255)
                              cATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCt                                                                                                                                       >  2:755937/1‑133 (MQ=255)
                              cATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCt                                                                                                                                       <  1:755937/133‑1 (MQ=255)
                                       tGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCgg                                                                                                             <  1:670380/150‑1 (MQ=255)
                                        gCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCggg                                                                                                            <  1:1036295/150‑1 (MQ=255)
                                         cGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCATATTCGTCGTGCGAGATTATAGCCGGGa                                                                                                           >  1:1224837/1‑150 (MQ=255)
                                          ggTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGAt                                                                                                          >  1:907804/1‑150 (MQ=255)
                                           gatGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATg                                                                                                         >  1:193/3‑150 (MQ=255)
                                              gACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCgag                                                                                                                        >  1:876700/1‑132 (MQ=255)
                                              gACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCgag                                                                                                                        <  2:876700/132‑1 (MQ=255)
                                                ccGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTc                                                                                                    <  1:872376/150‑1 (MQ=255)
                                                  gggCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCAt                                                                                                  >  1:93734/1‑150 (MQ=255)
                                                  gggCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCAt                                                                                                  >  2:973925/1‑150 (MQ=255)
                                                         gTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCg                                                                                           <  2:1255123/150‑1 (MQ=255)
                                                           ggggAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTa                                                                                         >  2:41478/1‑150 (MQ=255)
                                                           ggggAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATAATCGCCGGGATGGTGTCATCCTTGCGTa                                                                                         >  1:655325/1‑150 (MQ=255)
                                                              gagaGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCg                                                                                      <  1:549317/150‑1 (MQ=255)
                                                                      tgtTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAAtt                                                                              >  1:1370545/1‑150 (MQ=255)
                                                                         ttGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCg                                                                           <  2:296531/150‑1 (MQ=255)
                                                                              tACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGcc                                                                      <  2:930241/150‑1 (MQ=255)
                                                                                  gTTCTGGAAGGTTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTa                                                                  <  1:1274408/150‑1 (MQ=255)
                                                                                      tGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAg                                                              <  2:31197/150‑1 (MQ=255)
                                                                                       ggAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGt                                                             >  2:994390/1‑150 (MQ=255)
                                                                                          aGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTcac                                                          <  1:302344/150‑1 (MQ=255)
                                                                                             gTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTcaccac                                                       <  2:1230600/150‑1 (MQ=255)
                                                                                             gTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTcaccac                                                       <  2:907804/150‑1 (MQ=255)
                                                                                               aTTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACaa                                                     >  2:467070/1‑150 (MQ=255)
                                                                                                 tCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATa                                                   <  2:1370545/150‑1 (MQ=255)
                                                                                                  cATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAg                                                  >  2:1325295/1‑150 (MQ=255)
                                                                                                   aTTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGt                                                 >  1:340256/1‑150 (MQ=255)
                                                                                                      gtgCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATg                                              >  2:30653/1‑150 (MQ=255)
                                                                                                      gtgCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATg                                              >  2:774185/1‑150 (MQ=255)
                                                                                                       tgCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGc                                             <  2:320046/150‑1 (MQ=255)
                                                                                                         cgcCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTc                                           >  2:1301569/1‑150 (MQ=255)
                                                                                                                cATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATtgt                                    >  2:1417318/1‑150 (MQ=255)
                                                                                                                  tAAGATTAAATTATCGCATTAAAAAAGCCGGCGTTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGc                                  >  1:471967/1‑150 (MQ=255)
                                                                                                                  tAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGc                                  <  2:168637/150‑1 (MQ=255)
                                                                                                                           aTTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCt                         >  1:1037360/1‑150 (MQ=255)
                                                                                                                               tCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGa                     >  1:1392667/1‑150 (MQ=255)
                                                                                                                                 gCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGc                       >  1:530166/1‑146 (MQ=255)
                                                                                                                                 gCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGc                       <  2:530166/146‑1 (MQ=255)
                                                                                                                                   aTTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATTTTGTGCAGACAAGCTGCGATCtt                 >  2:1130824/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTgg              >  1:1238407/1‑150 (MQ=255)
                                                                                                                                         aaaGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGttt           >  2:860005/1‑150 (MQ=255)
                                                                                                                                             ccGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCTTATTGTGCAGACAAGCTGCGATCTTTGGTTTTTAc       >  2:293937/1‑150 (MQ=255)
                                                                                                                                             ccGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTAc       >  1:732747/1‑150 (MQ=255)
                                                                                                                                              cGGCATTCGCAGGCTTCATCTTATTCGTCGTGCGCTTTTATCGGCGGGTTTGTGTCATCGTTGCGTAGCGTAAGAATTTCGCATCCGTTATCATTAACCCCAATAGTATGCTCATATTTTGCATAAAAGCTGCGTTCTTTTGTTTTTAcc      <  2:1224837/150‑1 (MQ=255)
                                                                                                                                               ggCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCg     >  2:963715/1‑150 (MQ=255)
                                                                                                                                                  aTTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTcc  <  1:1296993/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGACGGTTAATTCATCACGGGGCCAGACGCATGGATTTTGCGGTTGACCGGGCAGGGTGGGGAGAGCGGTGTTTGCGTACTGTTCTGGAAGGGTATTCATTGTGCGCCACCCATAAGATTAAATTATCGCATTAAAAAAGCCGGCATTCGCCGGCTTCATCTTATTCGTCGTGCGAGATTATCGCCGGGATGGTGTCATCCTTGCGTAGCGTCAGAATTTCGCAGCCGTTATCAGTCACCACAATAGTATGCTCATATTGTGCAGACAAGCTGCGATCTTTGGTTTTTACCGTCC  >  NC_000913/188550‑188845

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