Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,541,077 G→C 7.1% A505G (GCC→GGC)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,541,0770GC 7.1% 267.7 / 10.2 84A505G (GCC→GGC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (33/45);  new base C (3/3);  total (36/48)
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.92e-01

TAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGG  >  NC_000913/1540930‑1541224
                                                                                                                                                   |                                                                                                                                                    
tAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCt                                                                                                                                                    >  2:343379/1‑150 (MQ=255)
     cATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgt                                                                                                                                               <  1:851226/150‑1 (MQ=255)
        gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgag                                                                                                                                            <  1:1234098/150‑3 (MQ=255)
         tCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTAc                                                                                                                                           >  1:1312186/1‑150 (MQ=255)
         tCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtgagc                                                                                                                                           <  2:92627/150‑4 (MQ=255)
          ccATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACg                                                                                                                                          <  2:1414426/150‑1 (MQ=255)
             tgtgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGtt                                                                                                                                       >  2:1088194/1‑150 (MQ=255)
               tgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                     <  2:105460/150‑1 (MQ=255)
               tgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTg                                                                                                                                     >  2:1510436/1‑150 (MQ=255)
                gATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGa                                                                                                                                    <  2:1383992/150‑1 (MQ=255)
                    ccAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTc                                                                                                                                >  1:662868/1‑150 (MQ=255)
                      aGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                                                                                                              >  1:608657/1‑150 (MQ=255)
                       gTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGc                                                                                                                             <  1:589095/150‑1 (MQ=255)
                        tGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCg                                                                                                                            <  2:757279/150‑1 (MQ=255)
                          gTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAt                                                                                                                                   >  2:1657490/1‑141 (MQ=255)
                          gTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAt                                                                                                                                   <  1:1657490/141‑1 (MQ=255)
                                  cTGCGCCGAGGATAATCATCGAGCGCCCATGGGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCg                                                                                                                  <  2:314912/150‑1 (MQ=255)
                                     cgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCAc                                                                                                               <  2:407235/150‑1 (MQ=255)
                                       ccGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTg                                                                                                             >  1:1108062/1‑150 (MQ=255)
                                        cGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑t                                                                                                           <  2:1383872/150‑1 (MQ=255)
                                          aGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑ttt                                                                                                         <  2:1530894/150‑1 (MQ=255)
                                           ggATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCTATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑tttt                                                                                                        <  1:1091889/150‑1 (MQ=255)
                                           ggATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑tttt                                                                                                        <  2:1628014/150‑1 (MQ=255)
                                            gATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTc                                                                                                       <  2:1312186/150‑1 (MQ=255)
                                             aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCa                                                                                                      >  1:1495143/1‑150 (MQ=255)
                                             aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCCATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCa                                                                                                      >  1:230007/1‑150 (MQ=255)
                                                 tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCtt                                                                                                  <  1:347588/150‑1 (MQ=255)
                                                   atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTcc                                                                                                <  2:1262047/150‑1 (MQ=255)
                                                   atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTcc                                                                                                >  2:747195/1‑150 (MQ=255)
                                                       aGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGcc                                                                                            >  2:18871/1‑150 (MQ=255)
                                                       aGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGcc                                                                                            <  1:343379/150‑1 (MQ=255)
                                                       aGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGcc                                                                                            >  2:49175/1‑150 (MQ=255)
                                                        gcgcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGccc                                                                                           >  2:490541/1‑150 (MQ=255)
                                                                 cggtttATGGGCGGTAGCGGCAAATTCACGGTCGATGGTTTCAATATACTGGAGCGGCACGGCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAg                                                                                  <  1:1061273/147‑1 (MQ=255)
                                                                 cGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAg                                                                                  <  2:662868/150‑1 (MQ=255)
                                                                  gTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAgg                                                                                 >  1:1229687/1‑150 (MQ=255)
                                                                      tacgggCGGAATCGGCAAATTCACGGGCGAGGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCGGGGGGGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGt                                                                             <  2:1324250/147‑1 (MQ=255)
                                                                       aTGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                            >  2:1490143/1‑150 (MQ=255)
                                                                        tGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCATTATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTATGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTaa                                                                           >  2:1523459/1‑150 (MQ=255)
                                                                            cGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCg                                                                       <  1:1088194/150‑1 (MQ=255)
                                                                              gTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCgtc                                                                     <  1:1474130/150‑3 (MQ=255)
                                                                              gTAGCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGcc                                                                     <  1:133704/150‑1 (MQ=255)
                                                                                 tCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGt                                                                  <  2:1190023/150‑1 (MQ=255)
                                                                                    gCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTAcc                                                               <  2:563577/150‑1 (MQ=255)
                                                                                     cAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCa                                                              >  1:1139814/1‑150 (MQ=255)
                                                                                         ttCACGGGCGATGGTTTCAATATACAGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATc                                                          >  1:1695113/1‑150 (MQ=255)
                                                                                          tCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAATTCa                                                         >  2:952543/1‑150 (MQ=255)
                                                                                           cACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt                                                        >  2:304194/1‑150 (MQ=255)
                                                                                           cACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt                                                        <  1:1088490/150‑1 (MQ=255)
                                                                                                gCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACg                                                   >  2:1115643/1‑150 (MQ=255)
                                                                                                  gATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGc                                                  >  2:584526/1‑150 (MQ=255)
                                                                                                   aTGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTg                                                <  1:490541/150‑1 (MQ=255)
                                                                                                   aTGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTg                                                <  2:1650707/150‑1 (MQ=255)
                                                                                                     ggTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGac                                              <  2:541012/150‑1 (MQ=255)
                                                                                                      gTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGacc                                             <  2:1229687/150‑1 (MQ=255)
                                                                                                        ttCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGaccac                                           <  2:608657/150‑1 (MQ=255)
                                                                                                             tataCTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGTGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGAc                                      <  1:49175/150‑1 (MQ=255)
                                                                                                                 cTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgt                                  >  2:1463136/1‑150 (MQ=255)
                                                                                                                  tGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtg                                 <  1:1510436/150‑1 (MQ=255)
                                                                                                                   ggCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                >  2:1707197/1‑150 (MQ=255)
                                                                                                                   ggCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                >  1:246329/1‑150 (MQ=255)
                                                                                                                    gcgcGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGtt                               >  2:1480256/1‑150 (MQ=255)
                                                                                                                        ggCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGGGTTGCCa                           <  1:1356771/150‑1 (MQ=255)
                                                                                                                         gCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCAt                          <  1:1690427/150‑1 (MQ=255)
                                                                                                                             gCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGAc                      <  1:1463136/150‑1 (MQ=255)
                                                                                                                              ccGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACt                     >  2:557401/1‑150 (MQ=255)
                                                                                                                               cGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTa                    >  1:414331/1‑150 (MQ=255)
                                                                                                                                 gTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAAc                  <  1:1115643/150‑1 (MQ=255)
                                                                                                                                  tAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACg                 <  2:1139814/150‑1 (MQ=255)
                                                                                                                                  tAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACg                 <  1:533153/150‑1 (MQ=255)
                                                                                                                                   aaTTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTTATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGt                <  1:18871/150‑1 (MQ=255)
                                                                                                                                    aTTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTc               >  2:311567/1‑150 (MQ=255)
                                                                                                                                      ttGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCATATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACGAACGTCAg             >  2:1699799/1‑150 (MQ=255)
                                                                                                                                      ttGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAg             >  1:975072/1‑150 (MQ=255)
                                                                                                                                       tGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAgg            <  2:585476/150‑1 (MQ=255)
                                                                                                                                       tGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAgg            >  1:690420/1‑150 (MQ=255)
                                                                                                                                         tacACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtt          >  1:940771/3‑150 (MQ=255)
                                                                                                                                             ccccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGcc                            >  1:877817/1‑128 (MQ=255)
                                                                                                                                             ccccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGcc                            <  2:877817/128‑1 (MQ=255)
                                                                                                                                             ccccAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTa      >  2:547530/1‑150 (MQ=255)
                                                                                                                                               ccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAAc    <  1:1210861/150‑1 (MQ=255)
                                                                                                                                                cAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACg   <  2:1695113/150‑1 (MQ=255)
                                                                                                                                                cAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACg   <  1:1441498/150‑1 (MQ=255)
                                                                                                                                                 aGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACgg  <  1:1174030/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTG‑TTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGG  >  NC_000913/1540930‑1541224

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: ‑