Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,926,036 A→C 5.6% E230A (GAA→GCA)  queF → 7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,0360AC 5.6% 392.6 / 12.6 108E230A (GAA→GCA) queF7‑cyano‑7‑deazaguanine reductase (NADPH‑dependent)
Reads supporting (aligned to +/- strand):  ref base A (51/51);  new base C (3/3);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 110 total aligned reads displayed.

GCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACC  >  NC_000913/2925895‑2926175
                                                                                                                                             |                                                                                                                                           
gCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTg                                                                                                                                       >  4:1116964/1‑148 (MQ=255)
             aTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGa                                                                                                                                 <  2:28309/141‑1 (MQ=255)
             aTCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGa                                                                                                                                 >  1:28309/1‑141 (MQ=255)
              tCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTg                                                                                                                                       <  4:1105923/134‑1 (MQ=255)
              tCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTg                                                                                                                                       <  2:337637/134‑1 (MQ=255)
              tCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTg                                                                                                                                       >  3:1105923/1‑134 (MQ=255)
                 aaCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGaa                                                                                                                                            >  1:966755/1‑126 (MQ=255)
                 aaCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGaa                                                                                                                                            <  2:966755/126‑1 (MQ=255)
                 aaCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTa                                                                                                                      >  2:948215/1‑148 (MQ=255)
                  aCTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTaa                                                                                                                     >  4:173270/1‑148 (MQ=255)
                    tGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATg                                                                                                                   <  1:975176/148‑1 (MQ=255)
                       cTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGaa                                                                                                                                >  1:707907/1‑132 (MQ=255)
                       cTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGaa                                                                                                                                <  2:707907/132‑1 (MQ=255)
                         gATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTg                                                                                                              >  4:683472/1‑148 (MQ=255)
                            cACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTa                                                                                                           <  3:1116964/148‑1 (MQ=255)
                              cccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACg                                                                                                         <  3:683472/148‑1 (MQ=255)
                               ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGt                                                                                                                                    >  3:710672/1‑120 (MQ=255)
                               ccATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGt                                                                                                                                    <  4:710672/120‑1 (MQ=255)
                                  tCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGc                                                                                                                             <  2:168337/124‑1 (MQ=255)
                                  tCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGc                                                                                                                             >  1:168337/1‑124 (MQ=255)
                                      ccAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCa                                                                                                                            <  2:416660/121‑1 (MQ=255)
                                      ccAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCa                                                                                                                            >  1:416660/1‑121 (MQ=255)
                                       cAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGCACGCa                                                                                                                            <  1:539181/120‑1 (MQ=255)
                                       cAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCa                                                                                                                            >  2:539181/1‑120 (MQ=255)
                                       cAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTgcca                                                                                                >  1:932151/1‑148 (MQ=255)
                                        aGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                                                             >  2:50149/1‑134 (MQ=255)
                                        aGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGt                                                                                                             <  1:50149/134‑1 (MQ=255)
                                                gTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAg                                                                                                                                         <  3:746788/98‑1 (MQ=255)
                                                gTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAg                                                                                                                                         >  4:746788/1‑98 (MQ=255)
                                                gTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaaa                                                                                       <  2:932151/148‑1 (MQ=255)
                                                    gCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTg                                                                                   <  3:173270/148‑1 (MQ=255)
                                                     cTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGa                                                                                  <  1:948215/148‑1 (MQ=255)
                                                       ccAAATTCAATATCGTGGACGCCAAATTGACAGAGAAAAACTTCTGCGTTACCTGGTTTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGc                                                                                >  4:21756/1‑148 (MQ=255)
                                                         aaaTTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGt                                                                              <  3:1002941/148‑1 (MQ=255)
                                                                  aTCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCtt                                                                                                      >  1:51152/1‑115 (MQ=255)
                                                                  aTCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCtt                                                                                                      <  2:51152/115‑1 (MQ=255)
                                                                    cGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATg                                                                                                                   >  1:406619/1‑100 (MQ=255)
                                                                    cGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATg                                                                                                                   <  2:406619/100‑1 (MQ=255)
                                                                            ccAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCt                                                                                                                         <  2:858856/86‑1 (MQ=255)
                                                                            ccAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCt                                                                                                                         >  1:858856/1‑86 (MQ=255)
                                                                              aaaTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGa                                                                                                                                 <  2:687421/76‑1 (MQ=38)
                                                                              aaaTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGCACTGTTCGTGGa                                                                                                                                 >  1:687421/1‑76 (MQ=39)
                                                                              aaaTTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCcgtcg                                                         >  2:545507/1‑148 (MQ=255)
                                                                                     cAGAGAAAAACTGCTGCGTTACCTGGTGTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAGCGCATCTTTAATGACCTGTTACGCtt                                                                                                      <  4:666355/96‑1 (MQ=255)
                                                                                     cAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCtt                                                                                                      >  3:666355/1‑96 (MQ=255)
                                                                                      agagAAAGACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTc                                                 <  3:21756/148‑1 (MQ=255)
                                                                                      agagAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTc                                                 <  1:545507/148‑1 (MQ=255)
                                                                                         gaAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGc                                                                                                                             <  2:664904/69‑1 (MQ=255)
                                                                                             aaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTa                                                                                                                      >  1:1057847/1‑72 (MQ=255)
                                                                                             aaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTa                                                                                                           <  4:531680/83‑1 (MQ=255)
                                                                                             aaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCg                                                                               >  2:244342/1‑111 (MQ=255)
                                                                                             aaCTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCg                                                                               <  1:244342/111‑1 (MQ=255)
                                                                                                 gctgcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAg                                                                                                                                         >  3:698045/1‑49 (MQ=255)
                                                                                                 gctgcGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAg                                                                                                                                         <  4:698045/49‑1 (MQ=255)
                                                                                                       ttACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                     >  1:635462/1‑79 (MQ=255)
                                                                                                       ttACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                     <  2:635462/79‑1 (MQ=255)
                                                                                                       ttACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaaaa                                                                                      >  1:874615/1‑94 (MQ=255)
                                                                                                       ttACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaaaa                                                                                      <  2:874615/94‑1 (MQ=255)
                                                                                                       ttACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTgg                                >  1:357803/1‑148 (MQ=255)
                                                                                                        tACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTta                                                                  <  2:785395/113‑1 (MQ=255)
                                                                                                        tACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTta                                                                  >  1:785395/1‑113 (MQ=255)
                                                                                                          ccTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGgcgc                             >  4:1002169/1‑148 (MQ=255)
                                                                                                            tGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaaaa                                                                                      <  2:299876/89‑1 (MQ=255)
                                                                                                            tGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGaaaaa                                                                                      >  1:299876/1‑89 (MQ=255)
                                                                                                            tGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTa                                        >  2:402930/1‑135 (MQ=255)
                                                                                                            tGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTa                                        <  1:402930/135‑1 (MQ=255)
                                                                                                               tctcATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTaa                        >  4:63906/1‑148 (MQ=255)
                                                                                                                  cATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                     >  3:950617/1‑68 (MQ=255)
                                                                                                                  cATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                     <  4:525588/68‑1 (MQ=255)
                                                                                                                  cATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTc                                                                                                     >  3:525588/1‑68 (MQ=255)
                                                                                                                    ttCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTaa                                       <  1:395184/128‑1 (MQ=255)
                                                                                                                    ttCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCg                   <  3:63906/148‑1 (MQ=255)
                                                                                                                         tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACCTTAACCCGTg                                 <  3:653163/129‑1 (MQ=255)
                                                                                                                         tcatcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTg                                 >  4:653163/1‑129 (MQ=255)
                                                                                                                          catcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAg                     >  4:1091117/1‑140 (MQ=255)
                                                                                                                           atcaCAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGt            <  2:357803/148‑1 (MQ=255)
                                                                                                                             cacaACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGAcc                                                                                                                >  4:304248/1‑46 (MQ=255)
                                                                                                                             cacaACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGAcc                                                                                                                <  3:304248/46‑1 (MQ=255)
                                                                                                                             cacaACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTcc          >  1:205570/1‑148 (MQ=255)
                                                                                                                                aaCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCAt       <  2:205570/148‑1 (MQ=255)
                                                                                                                                 aCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCa                                                                       >  4:822871/1‑83 (MQ=255)
                                                                                                                                 aCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCa                                                                       <  3:822871/83‑1 (MQ=255)
                                                                                                                                 aCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTtat                                                                 >  4:1140962/1‑89 (MQ=255)
                                                                                                                                 aCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTtat                                                                 <  3:1140962/89‑1 (MQ=255)
                                                                                                                                 aCGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATc      >  1:259844/1‑148 (MQ=255)
                                                                                                                                   gAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg                                              <  1:264675/106‑1 (MQ=255)
                                                                                                                                   gAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTgg                                              >  2:264675/1‑106 (MQ=255)
                                                                                                                                     gTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGAcc  <  3:1002169/148‑1 (MQ=255)
                                                                                                                                      ttCCACGCACTGTTCGTGGAACTCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCcg                                                            >  4:429084/1‑89 (MQ=39)
                                                                                                                                      ttCCACGCACTGTTCGTGGAACTCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCcg                                                            <  3:429084/89‑1 (MQ=38)
                                                                                                                                      ttCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAAt                       >  4:320550/1‑126 (MQ=255)
                                                                                                                                      ttCCACGCACTGTTCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAAt                       <  3:320550/126‑1 (MQ=255)
                                                                                                                                      ttCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATa                      <  3:452035/127‑1 (MQ=255)
                                                                                                                                          aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTtat                                                                 <  4:1020795/80‑1 (MQ=255)
                                                                                                                                          aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTtat                                                                 >  3:1020795/1‑80 (MQ=255)
                                                                                                                                          aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCa        <  1:257435/137‑1 (MQ=255)
                                                                                                                                          aCGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCa        >  2:257435/1‑137 (MQ=255)
                                                                                                                                           cGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAg                           <  1:177441/117‑1 (MQ=255)
                                                                                                                                           cGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAg                           >  2:177441/1‑117 (MQ=255)
                                                                                                                                           cGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTccc         >  1:699881/1‑135 (MQ=255)
                                                                                                                                             |                                                                                                                                           
GCCACCTGCTGAAATCAAACTGCCTGATCACCCATCAACCAGATTGGGGTTCGCTCCAAATTCAGTATCGTGGACGCCAAATTGACAGAGAAAAACTGCTGCGTTACCTGGTCTCATTCCGTCATCACAACGAGTTCCACGAACAGTGCGTGGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACC  >  NC_000913/2925895‑2926175

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