Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,292 C→T intergenic (+96/‑16) queF → / → ppnN 7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,2920CT93.0% 158.2 / 1.5 57intergenic (+96/‑16)queF/ppnN7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (26/27);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAG  >  NC_000913/2926047‑2926526
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
gAACGCATCTTTAATGACCTTTTACGCTTCTGCCAGCAAGAAAAATTGAGCGTTTACGCACGGTATACCCGTCGTGGCGGTCTTGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATaa                                                                                                                                                                                                                                       <  2:661033/251‑1 (MQ=255)
 aaCGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAg                                                                                                                                                                                                                                      <  1:492684/251‑1 (MQ=255)
      aTCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAACGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGtt                                                                                                                                                                                                                                 >  1:17202/1‑251 (MQ=255)
          ttAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTGTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAAGAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTAGAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGtt                                                                                                                                                                                                                                 <  2:226721/247‑1 (MQ=255)
                  cTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTa                                                                                                                                                                                                                      >  1:121753/1‑250 (MQ=255)
                              cgccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCTACTTTGTCCCATCGACCACAAGACTTGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGccc                                                                                                                                                                                                         <  1:235809/250‑1 (MQ=255)
                               gccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGccc                                                                                                                                                                                                         <  1:428075/250‑1 (MQ=255)
                                  agccagAAAAAATGATCGTTTACGCCCGGTATTCCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCtt                                                                                                                                                                                                     <  2:508468/251‑1 (MQ=255)
                                   gccagAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTg                                                                                                                                                                                                    <  1:546344/251‑1 (MQ=255)
                                                  cGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgt                                                                                                                                                                                     <  1:420214/251‑1 (MQ=255)
                                                   gTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                    >  1:617994/1‑251 (MQ=255)
                                                       aCGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCg                                                                                                                                                                                <  1:65369/251‑1 (MQ=255)
                                                          cACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGc                                                                                                                                                                               >  1:367443/1‑249 (MQ=255)
                                                                ataCCCGTCGTGGCGGGCTGGACATTATCCCGTGGCGCAGTAATCGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGa                                                                                                                                                                        <  2:515462/250‑1 (MQ=255)
                                                                    ccGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTgg                                                                                                                                                                   <  1:127189/251‑1 (MQ=255)
                                                                      gtcgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGAt                                                                                                                                                                                                                        <  1:629344/196‑1 (MQ=255)
                                                                        cgtGGCGGTCCGTACATTAACCCGTGGCGGAAGAATAGCGGTTTTGTCCCATCGACCACAAAACTGGGTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGCTTCACGCAGACGGA‑TGTGAAAGATCATCAGGCCGGGCTATTGTAATCCAAGGGAAAGACGATATTCTTCTCAAAAGGAGTTTTTCTTTATTCCACATATTAGACCGCTTGGCTCCAGAGAGATGTTGTCGCAGCTGGAAGTGGaaat                                                                                                                                                               >  2:153615/1‑248 (MQ=255)
                                                                              ggTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTGTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaa                                                                                                                                                         <  1:697542/251‑1 (MQ=255)
                                                                                              cGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACGAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCagcagc                                                                                                                                         >  1:270816/1‑251 (MQ=255)
                                                                                                      aGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTa                                                                                                                                 <  1:627940/251‑1 (MQ=255)
                                                                                                       gTAATAGCGATTTTGTCCCATCGACCACAAGACGGGTTCGGCAATAAATTTTTTCTAAATTTTGCGTGTTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTAt                                                                                                                                <  1:372575/251‑1 (MQ=255)
                                                                                                              cGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGt                                                                                                                         >  1:4307/1‑251 (MQ=255)
                                                                                                                  tttGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTACTCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCg                                                                                                                     >  2:34254/1‑251 (MQ=255)
                                                                                                                      tCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGCAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCCTCAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAAc                                                                                                                 >  2:408779/1‑251 (MQ=255)
                                                                                                                             gACCACAAGACTGGTTCGGCAATAACTTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGACGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACTTCTATCAACTGTTTCGCAACTGTTCAc                                                                                                          >  2:487769/1‑251 (MQ=255)
                                                                                                                              aCCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGTAAGGGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGAAACTGTTCACt                                                                                                         >  2:145054/1‑251 (MQ=255)
                                                                                                                                  caAGACTGGATCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGATTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGCTATGCTTAAACGCACCGCCCGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGcc                                                                                                     >  2:532063/1‑251 (MQ=255)
                                                                                                                                        ggtttCGGAAATAAATTTTTTCTAAATTTTGCTGGCGGGATTCCGGCAGA‑GGTGTGTGAAAGGGCATAAGGCAGGGCTATTGTACACAAGGGGAATGTCTTTTGTCGTCTCATAAGGAGTTTTTCTTGATGACACATATTCGCCCGGTTGGGTCCATGGATAGGTTGTCGCGGCTGAATGTGGATATGCTTAACCGACCCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACGTGCCGTACTg                                                                                               <  2:275055/248‑1 (MQ=255)
                                                                                                                                           ttCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCACCTGTTTCGCACCTGTTCACTTGCCGTACTGAAc                                                                                            >  2:327423/1‑251 (MQ=255)
                                                                                                                                            tCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGGCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAAc                                                                                            <  1:153615/250‑1 (MQ=255)
                                                                                                                                              ggCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTcc                                                                                         >  1:164702/1‑251 (MQ=255)
                                                                                                                                                  aTAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAAc                                                                                                                            >  2:360283/1‑212 (MQ=255)
                                                                                                                                                  aTAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAAc                                                                                                                            <  1:360283/212‑1 (MQ=255)
                                                                                                                                                      atttgttcgtaaTTTTGCGTGCTGGATTCACGCAGAATGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTGGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGTCCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGttt                                                                                 <  1:640714/241‑1 (MQ=255)
                                                                                                                                                                     tgtgggtgggatttaggCAGAAGGT‑TGTGAAAGGCCTCCAGGCAGGGGTATTGTAATCAAAGGGAATGACGATATTGTTCTCATAAGGAGTTTTTCTTGATGACAAATATTGGCCCGTTTGGCTCCATGGATATGGTGTCGCAGCTGGAAGGGGATGTGCTTAAACCCCCCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAgtt                                                                   <  2:135095/234‑3 (MQ=255)
                                                                                                                                                                     tGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGACATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCa                                                                   >  2:691809/1‑250 (MQ=255)
                                                                                                                                                                        gTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAAGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCTGCCTGAACTACGGTAGTTTGACCGATacaagcaa                                                                  >  2:471089/1‑241 (MQ=255)
                                                                                                                                                                               aTTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTTTAATCAAAGGGAATGACGATATGCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAGTTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGaa                                                                                             <  2:512683/214‑1 (MQ=255)
                                                                                                                                                                               aTTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGaa                                                                                             >  1:512683/1‑214 (MQ=255)
                                                                                                                                                                                ttCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGt                                                       >  2:183571/1‑251 (MQ=255)
                                                                                                                                                                                     gCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGt                                                  >  1:374547/1‑251 (MQ=255)
                                                                                                                                                                                       aGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCACGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGATGTCTCGttt                                                >  2:411523/1‑251 (MQ=255)
                                                                                                                                                                                          aGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGaa                                                                                             >  1:632459/1‑203 (MQ=255)
                                                                                                                                                                                               tttgAAAGTTCATCCGGCAGGGCTATTTTAATAAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGTCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAACCGCACCGCCAGCAGCGACCTCTATCAAc                                                                                                                            <  2:155798/166‑1 (MQ=255)
                                                                                                                                                                                               tgtgAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAAc                                                                                                                            >  1:155798/1‑168 (MQ=255)
                                                                                                                                                                                                   aaaGGTCATCAGGAAGGGCTATTTTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGCCGCAGCTGGAAGTGGATATGCTTAAACGAACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCg                                     <  1:549309/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTCATAAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCTCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGata                                  <  2:270816/251‑1 (MQ=255)
                                                                                                                                                                                                        tcatcaGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACCGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATAtt                                >  2:505080/1‑251 (MQ=255)
                                                                                                                                                                                                          atcaGGCAGGGCTACTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGCTTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAACATTTCGATATTa                               >  2:292970/1‑250 (MQ=255)
                                                                                                                                                                                                             aggcaggGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAGGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGt                           <  1:34254/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTgcg                     <  1:183571/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aatcaaaTGACGATATTCGTCTCAAAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACgcggcg       <  2:48075/246‑1 (MQ=255)
                                                                                                                                                                                                                                     gAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATCTGTTGTCGCAGATGGAAGTGGATATGCTTAAAc                                                                                                                                                       >  2:551787/1‑103 (MQ=255)
                                                                                                                                                                                                                                     gAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTAGGCTGCATGGATATGGTGGCGCAGCTGGAAGTGGATATGCTTAAAc                                                                                                                                                       <  1:551787/103‑1 (MQ=255)
                                                                                                                                                                                                                                      aaggACGATATTCGTCTCATAAGGAGTTTTTCTTGATAACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGGGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAg  <  1:292970/248‑1 (MQ=255)
                                                                                                                                                                                                                                      aaTGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACgcggcg       >  1:48075/1‑246 (MQ=255)
                                                                                                                                                                                                                                      aaTGACGATATCCGTACTATAAGGAGTTTTTCTGGATTACACATATTAGCCCTCTTGGCTCCAGGGATAGGTTGTCGCAGCGGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGTCGGGAACGCGGCGTAAAg  <  1:201217/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
GAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGT‑TGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAG  >  NC_000913/2926047‑2926526

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑