Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,279,415 T→A 5.4% intergenic (+67/‑449) narK → / → narG nitrate/nitrite transporter/nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,279,4150TA 5.4% 264.9 / 2.8 84intergenic (+67/‑449)narK/narGnitrate/nitrite transporter/nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (37/42);  new base A (3/2);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTAT  >  NC_000913/1279267‑1279562
                                                                                                                                                    |                                                                                                                                                   
cgcAATGAAGTTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCacttactt                                                                                                                                                     <  6:321364/149‑1 (MQ=255)
cgcAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCacttactt                                                                                                                                                     <  7:225582/149‑1 (MQ=255)
  cAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGTTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGTTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTTACTTACTTTa                                                                                                                                                   <  4:86621/149‑1 (MQ=255)
  cAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCacttactt                                                                                                                                                     >  1:44940/1‑147 (MQ=255)
  cAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTa                                                                                                                                                   <  4:328157/149‑1 (MQ=255)
    aTGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGt                                                                                                                                                 <  6:188065/149‑1 (MQ=255)
    aTGAACGTATTTTTGAATTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATGCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGt                                                                                                                                                 <  1:283014/149‑1 (MQ=255)
        aGGTATTTTTGATTTTCTATATCGCCTGCGAAGTGATAACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAg                                                                                                                                              >  4:276174/1‑148 (MQ=255)
          gTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTc                                                                                                                                           >  2:292290/1‑149 (MQ=255)
             tttttGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTg                                                                                                                                        <  8:234030/149‑1 (MQ=255)
               tttGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTggg                                                                                                                                      >  4:226471/1‑149 (MQ=255)
               tttGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTACGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTgg                                                                                                                                       >  1:327497/1‑148 (MQ=255)
                    ttttCTATATCGCCTGCGTAGTGATTAACTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACAGTTACTCCTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGattca                                                                                                                                 >  2:138971/1‑149 (MQ=255)
                     tttCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGattcat                                                                                                                                >  3:108173/1‑149 (MQ=255)
                      ttCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGattcatt                                                                                                                               >  6:96455/1‑149 (MQ=255)
                      ttCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGattcat                                                                                                                                <  3:54334/148‑1 (MQ=255)
                          atcccGCCTGCCTCGTGATTACCTGGGCGGTATATGGTAGGCCGTCTACAAAATAACCCGTTACTCGTCATACTTCGGGTTACATGGGCTGCGGCTGCGGTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACt                                                                                                                           <  8:239345/145‑1 (MQ=255)
                           tatCTCCTGCGTAGTGATTACCTGTGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACtt                                                                                                                          >  6:135101/1‑149 (MQ=255)
                           tatCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGAttcattca                                                                                                                             <  7:253439/146‑1 (MQ=255)
                           tatCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACtt                                                                                                                          >  7:66038/1‑149 (MQ=255)
                             tCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTg                                                                                                                         <  7:74375/148‑1 (MQ=255)
                             tCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTg                                                                                                                         <  3:226471/148‑1 (MQ=255)
                              cGCCTTCGTAGTGATTACCTGGGCGGTATATGGTCGGCTTTCTAAAAATTTAACCGTTACTCGTCATTCTTCGGGTTTCATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTgc                                                                                                                        <  1:255320/148‑1 (MQ=255)
                               gCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTgccg                                                                                                                      >  7:113081/1‑149 (MQ=255)
                                ccTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCGCTTACTAAAGTAAgtgactg                                                                                                                                        <  7:149021/130‑7 (MQ=255)
                                ccTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCGCTTACTAAAGTAAgtgactg                                                                                                                                        >  8:149021/1‑124 (MQ=255)
                                  tGCGTAGTGATTACCTGGGCGGTATATGGTCTGCATTCAAAAAAATAAACCGTTACTCGTCATACTTCGGGTTTCATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCCCTTGCCGACt                                                                                                                   >  5:84580/1‑149 (MQ=37)
                                    cGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCGCTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTc                                                                                                                 <  5:36389/149‑1 (MQ=255)
                                      tAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCt                                                                                                               <  8:221730/149‑1 (MQ=255)
                                      gaGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAACTAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCt                                                                                                               <  2:247281/148‑1 (MQ=255)
                                                cTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTATAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGaa                                                                                                     >  3:67071/1‑149 (MQ=255)
                                                   ggCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATta                                                                                                  >  1:275009/1‑149 (MQ=255)
                                                    gCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATta                                                                                                  >  5:131649/1‑148 (MQ=255)
                                                       gTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATtatata                                                                                              <  4:108173/149‑1 (MQ=255)
                                                       gTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATtatata                                                                                              <  5:96455/149‑1 (MQ=255)
                                                       gTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGACTTCCTGTAAACCGAATtatata                                                                                              <  2:207977/149‑1 (MQ=255)
                                                          tatGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCTTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATaga                                                                                            <  5:86662/148‑1 (MQ=255)
                                                            tGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGt                                                                                          >  5:14514/1‑148 (MQ=255)
                                                              gTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTAATTTAGTAAGCTCCTGGGGTTCATTCACTTGCCGCCTTCCTGTAAACAGAAGTATATAGAGTaa                                                                                        >  5:228439/1‑148 (MQ=255)
                                                               tCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTaaaa                                                                                      <  6:14514/149‑1 (MQ=255)
                                                                    aTTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTAACTTGCCGCCTTAATGTAAACCGAATTATATATAGTAAAATAttt                                                                                 >  2:351933/1‑149 (MQ=255)
                                                                     ttCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATAttt                                                                                 >  7:305857/1‑148 (MQ=255)
                                                                        tAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGAt                                                                              >  1:143945/1‑148 (MQ=255)
                                                                        tAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGAt                                                                              >  2:54512/1‑148 (MQ=255)
                                                                          aaaaaTACACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACAATAGAAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAAATATATAGAGAAAAAAATTAGATTAt                                                                           >  6:188289/1‑149 (MQ=255)
                                                                          aaaaaTAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTa                                                                            <  7:285087/148‑1 (MQ=255)
                                                                          aaaaaTAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTa                                                                            <  3:276174/148‑1 (MQ=255)
                                                                          aaaaaTAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTAt                                                                           >  5:325316/1‑149 (MQ=255)
                                                                          aaaaaTAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTAt                                                                           >  2:75865/1‑149 (MQ=255)
                                                                          aaaaaTAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTAt                                                                           >  6:318307/1‑149 (MQ=255)
                                                                               tAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCttt                                                                      <  8:282586/149‑1 (MQ=255)
                                                                                aaaCCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTg                                                                     >  6:34448/1‑149 (MQ=255)
                                                                                 aaCCTTTACTCGTCATACTTCGGGTTACATGTGCTGAGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGAAAGCTCCTGGGATTAATTCAATTGCCGCCTTCCAGTAAACCGAATAATAAAGAGTAAAAAATTTGATTATCCttt                                                                      >  5:160376/1‑147 (MQ=25)
                                                                                    cGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAGATATTTGATTATCCTTTGcgcg                                                                 >  7:308470/1‑149 (MQ=255)
                                                                                            gTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCatgatg                                                         >  8:196836/1‑149 (MQ=255)
                                                                                             tCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATAttt                                                                                 <  4:56045/124‑1 (MQ=255)
                                                                                             tCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATAttt                                                                                 >  3:56045/1‑124 (MQ=255)
                                                                                               aTACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTc                                                       <  1:75865/148‑1 (MQ=255)
                                                                                                tACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTcgc                                                     <  8:107555/149‑1 (MQ=255)
                                                                                                 aCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTcgcg                                                    >  7:295144/1‑149 (MQ=255)
                                                                                                  cTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTAAGTTTGTAAGCGACGGGGATTCATTCACTTGCCGCCGTCCGGTAAACCGCATTATATAGAGTAAAATATTTTATTATACTTTGCGCGGCATGATTTagcgg                                                   >  8:57058/1‑148 (MQ=255)
                                                                                                      gggTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCtttt                                               <  2:275009/149‑1 (MQ=255)
                                                                                                       ggTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCt                                                                                                               >  7:105976/1‑84 (MQ=255)
                                                                                                       ggTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCt                                                                                                               <  8:105976/84‑1 (MQ=255)
                                                                                                        gTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCtttttt                                             <  8:305857/149‑1 (MQ=255)
                                                                                                            cATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTCACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTAt                                          <  7:12764/148‑1 (MQ=255)
                                                                                                             aTGTGCTGCGGCTGCGTTGATTGACCCCAGTCACTTACTTTAGTTTGCTCGTGGGTTTGTTTCAGTTGCGGCCTTCCTGTAATCGGTATTATATAGAGTATAATATTTGATTATGCTTTGGGCGGCATGATGTCGCGCTTTTTTTATgg                                        <  8:113081/149‑2 (MQ=255)
                                                                                                             aTGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCatga                                                           >  8:303623/1‑130 (MQ=255)
                                                                                                             aTGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCatga                                                           <  7:303623/130‑1 (MQ=255)
                                                                                                             aTGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGc                                        >  5:137045/1‑149 (MQ=255)
                                                                                                              tgtgCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCg                                       <  1:90756/149‑1 (MQ=255)
                                                                                                               gtgCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGt                                      >  5:15740/1‑149 (MQ=255)
                                                                                                                tgctgcGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTTGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTc                                     <  1:196239/149‑1 (MQ=255)
                                                                                                                  ctgcggctgcgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCAt                                   <  7:196836/149‑1 (MQ=255)
                                                                                                                     cggctgcgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCAtt                                  <  6:137045/147‑1 (MQ=255)
                                                                                                                       gctgcgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATGCATTCACTGGCCGCCTTCCAGTAAACCGCATTATATAGAGTACAATAGTTGCTTATCCTTTGCTCGGCATGATGTCGCGCTTCTTTTCTGCGGAATTTAGt                              >  3:262669/1‑149 (MQ=255)
                                                                                                                        ctgcgTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCAGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGtt                             <  1:138971/149‑1 (MQ=255)
                                                                                                                                attcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATa                     >  1:54413/1‑149 (MQ=255)
                                                                                                                                attcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATa                     <  8:295144/149‑1 (MQ=255)
                                                                                                                                attcaCCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATa                     <  5:34448/149‑1 (MQ=255)
                                                                                                                                     ccccAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATTGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAAt                >  3:19833/1‑149 (MQ=255)
                                                                                                                                     ccccAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAAt                <  2:327497/149‑1 (MQ=255)
                                                                                                                                      ccccGTCACATACCTTCCCAAGCTCCAGGGCTGCCTACACAAGCCGCCCTCCTCTCAACCCACTTATATCGCGGAAAAAACTTCACTATCCTTTGCCCGGCATGAGGTCGCGCTTTTTTTATGCGCCACTTAGTTACACCATACTAATg               <  8:319240/149‑1 (MQ=255)
                                                                                                                                        cAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGtt             <  7:42352/149‑1 (MQ=255)
                                                                                                                                           tCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAAt                >  6:247161/1‑143 (MQ=255)
                                                                                                                                           tCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAAt                <  5:247161/143‑1 (MQ=255)
                                                                                                                                           tCACTTACTATAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTtata          <  4:67071/149‑1 (MQ=255)
                                                                                                                                                   tttAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTAt  >  7:48605/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAAAATATTTGATTATCCTTTGCGCGGCATGATGTCGCGCTTTTTTTATGCGTCATTTAGTTACAACATACTAATGTTATATGGTTTAT  >  NC_000913/1279267‑1279562

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