Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,026,650 G→A D42N (GAT→AAT)  fre → NAD(P)H‑flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,026,6500GA98.6% 241.4 / ‑6.7 72D42N (GAT→AAT) freNAD(P)H‑flavin reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (38/33);  minor base T (1/0);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGC  >  NC_000913/4026410‑4026879
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
ccAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGcgcg                                                                                                                                                                                                                              <  1:410181/250‑1 (MQ=255)
ccAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGa                                                                                                                                                                                                                             <  1:385749/251‑1 (MQ=255)
      gttgtCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACg                                                                                                                                                                                                                       >  2:86058/1‑251 (MQ=255)
        tgtCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGt                                                                                                                                                                                                                      <  1:266954/250‑1 (MQ=255)
         gtCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTcc                                                                                                                                                                                                                    >  1:250900/1‑251 (MQ=255)
              gcATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAACCGTCCGTtct                                                                                                                                                                                                               >  1:198018/1‑251 (MQ=255)
              gcATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTtct                                                                                                                                                                                                               >  1:504810/1‑251 (MQ=255)
              gcATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTtct                                                                                                                                                                                                               <  1:468916/251‑1 (MQ=255)
                atatTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCATGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGACGCTATCACGGATACCGTATATCGTGTCCGCATAGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGttatc                                                                                                                                                                                                              >  1:261337/1‑247 (MQ=255)
                  atTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGCCAAACGTCCGTTCTCAAt                                                                                                                                                                                                           >  2:174067/1‑251 (MQ=255)
                          cATCTGTGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTGTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGa                                                                                                                                                                                                   <  1:351730/251‑1 (MQ=255)
                               gggATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTGCTTTTCGTGCTGGTCAGTATTTGAGGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGcc                                                                                                                                                                                               <  2:534951/250‑1 (MQ=255)
                                ggATGAACTGGCTATTTTTAACATCGGTAAAAGCGCCTGATGCGCGTTTGTTGTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACACCCTTAAGCTGTAAAGTGACCGCGGCAGAAGCTATCAGGGATACCGGATATCGTGTCCGCATCGTGCCAGACGCGGTCCTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAATCGTCCGTTCTCAATGGCTTCTACGCCgg                                                                                                                                                                                             <  1:354199/251‑1 (MQ=255)
                                ggATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGAAGAGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCgg                                                                                                                                                                                             >  2:71630/1‑251 (MQ=255)
                                     aaCTGGATATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACGACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGAATTTGATGGTAGTGATGACTGAGCGCGACAAACGTCCGGTCTCAATGGCGTCGACGCCGGATgc                                                                                                                                                                                         >  2:369458/1‑249 (MQ=255)
                                        tGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCACTGGCTTCGACGCCGGATGaaac                                                                                                                                                                                      >  1:196505/1‑249 (MQ=255)
                                             aTTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAggg                                                                                                                                                                                  <  1:404071/249‑1 (MQ=255)
                                                tttAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCGCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTa                                                                                                                                                                              <  2:293958/250‑1 (MQ=255)
                                                 ttAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTAt                                                                                                                                                                             <  2:99387/250‑1 (MQ=255)
                                                    acaacGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGGGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCg                                                                                                                                                                           >  1:159083/1‑249 (MQ=255)
                                                     caacGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATAGATCCGCCAGAGAAAGAGCATGACAAACTTAGGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCAGCATCGTGACAGACGCGGCCTTTTTTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGTATGGACGCGACACACGTACGTTCTCAATGGCTTCGACGCCGGCTGAAAATGGGTTTATCGAGc                                                                                                                                                                        >  1:549782/1‑251 (MQ=255)
                                                            aaaaGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACACACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCatat                                                                                                                                                                 >  2:447898/1‑251 (MQ=255)
                                                            aaaaGAGCCTGATGCGCGTTTGTTTTGCCCTATTTAACGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGTCCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGTGGCCTTTTCTTTGCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCTCGACAAACGTCCGTTCTCAATGGCTTCTACGCCGGATGAAAAAGGGTTTATCGAGCTGCatat                                                                                                                                                                 <  2:175939/251‑1 (MQ=255)
                                                             aaaGCGCCTGATGCGCGTGTGTTTTGCCCTATTTATCGATCCGACAGATAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGGTCTCTATGGCTGCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATAtt                                                                                                                                                                <  2:261337/251‑1 (MQ=255)
                                                               aGCGCCTGATGAGCGTTTGTTTTGCCCCAATTATCGATCCGACAGCGACAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGCTACCGTATATCGTGTCCGCATCGTGCCAGACGCGGACTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGCTGAAAACGGGTTTATCGAGCTGCATATTgg                                                                                                                                                              >  2:553118/1‑251 (MQ=255)
                                                                      gATGCGCGGTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTc                                                                                                                                                        <  1:185217/250‑1 (MQ=255)
                                                                                tgtttTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAAc                                                                                                                                              <  1:126066/250‑1 (MQ=255)
                                                                                tgtttTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAAc                                                                                                                                              <  1:158523/250‑1 (MQ=255)
                                                                                     tGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTAc                                                                                                                                        >  2:137469/1‑251 (MQ=255)
                                                                                       cccTATTTATCGATCCGACAGAGACAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTAcgc                                                                                                                                      >  2:273141/1‑251 (MQ=255)
                                                                                          tATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAACGTGCCCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCAGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGTGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAACGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGa                                                                                                                                    >  1:484189/1‑250 (MQ=255)
                                                                                          tATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGa                                                                                                                                    >  1:351534/1‑250 (MQ=255)
                                                                                                     ccGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATg                                                                                                                         <  1:155051/250‑1 (MQ=255)
                                                                                                               aaGCGCATGACAACCTTAAGCGGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTAGATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTGTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGGCATGGACCGCATCCt                                                                                                              <  2:196505/251‑1 (MQ=255)
                                                                                                                            ccTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCACGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGAtcatca                                                                                                  >  1:78269/1‑250 (MQ=255)
                                                                                                                             cTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAACAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCaa                                                                                                 >  2:201334/1‑250 (MQ=255)
                                                                                                                                      tAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCCTCGTGCCCGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTTTTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAACCGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGCACGCAGTCATGGACAGCATCCTCAAAGATCATCAAATCGTGGTCg                                                                                       >  2:280051/1‑251 (MQ=255)
                                                                                                                                           tGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTAGAcctt                                                                                  >  2:449974/1‑248 (MQ=255)
                                                                                                                                                 cGGTAGAAGCTATCACGGATTCCGTATATCGTGTCCGCATCGTGCCAGACTCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGTATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCAc                                                                            <  1:369458/251‑1 (MQ=255)
                                                                                                                                                  ggTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACg                                                                           <  1:108636/251‑1 (MQ=255)
                                                                                                                                                     aGAAGCTATCACGGATACCGTATATCGTGTCCGCATAGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTCTTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTAGACGCCGGATGAAAAAGGGTTTATAGCGCTGCTTATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGAATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACgg                                                                          >  1:493537/1‑249 (MQ=255)
                                                                                                                                                       aaGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGaa                                                                      <  1:221811/251‑1 (MQ=255)
                                                                                                                                                       aaGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTAATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAACAGGGTTTATCGAGCTGCATATTTGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCAAGGACCGCATCCTCAAAGATCATCACATCGTGGTCGACATTCCCCACGgaga                                                                       >  2:136991/1‑250 (MQ=255)
                                                                                                                                                       aaGCTATAACGGATAACGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTATGAAATCAACCTTTACGCGAAAGCAGTCCTGGACCGCATCCTCAACGATCATCAAATCGTGGTCGACATTCACCAAGgaga                                                                       >  2:544622/1‑250 (MQ=255)
                                                                                                                                                          cTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatc                                                                                                      <  1:46883/216‑1 (MQ=255)
                                                                                                                                                          cTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTATTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGCAATCAACCTTTACGCGAAAGCAGTCTTGGACCGCATCCTCAAAgctc                                                                                                      >  2:46883/1‑213 (MQ=255)
                                                                                                                                                           tATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAACAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGt                                                                  >  1:108538/1‑251 (MQ=255)
                                                                                                                                                           tATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGCCATTCCCCACGGAGAAGCGt                                                                  >  2:320151/1‑251 (MQ=255)
                                                                                                                                                              cACGGATACCGGATATCGTGTCCGCAGCGTGCCAGACGTGGCCTTTTCTTTGCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGGGCGCGACAAACGTCCGGGCTCAATGGCTTCGTGGCCGTATGAAAAAGGGTTTATCGGGCTGCATATTGGCGCCTCTGAAATCAACCTTTTCGCGAAAGCAGTCATGGACCGCACCCTTAAATATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTgg                                                                <  1:497290/250‑1 (MQ=255)
                                                                                                                                                               aCGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATTAGCGCGACAAACGTCCGTTCCCAATGGCGTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCGTCTTAAATCAACCTTTACGCTAAATCATTCATGGACCGCATCCTCAAAGATCAGCAAATCGTGTTCTACATTCCCCACGTTGAAGCGTGGCt                                                              <  2:314344/251‑1 (MQ=255)
                                                                                                                                                                 ggATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTgc                                                            <  1:273141/251‑1 (MQ=255)
                                                                                                                                                                    tACCGTATATCGTGTCCGCATCGTGCCAGCCGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGcgcg                                                         >  2:378282/1‑251 (MQ=255)
                                                                                                                                                                      ccGTATTTCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTGCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGTGTGGCTGCGCgat                                                       <  2:183713/251‑1 (MQ=255)
                                                                                                                                                                       cGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGCTGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGAATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatg                                                      >  2:505685/1‑251 (MQ=255)
                                                                                                                                                                            atCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGAtgatga                                                  <  1:137469/250‑1 (MQ=255)
                                                                                                                                                                            atCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGaa                                                 <  2:472762/251‑1 (MQ=255)
                                                                                                                                                                                    cGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTcc                                         <  1:463529/251‑1 (MQ=255)
                                                                                                                                                                                       aTCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCgat                                      >  2:116873/1‑251 (MQ=255)
                                                                                                                                                                                             ccAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGAtttt                                >  1:522367/1‑251 (MQ=255)
                                                                                                                                                                                              cAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTg                               >  1:11377/1‑251 (MQ=255)
                                                                                                                                                                                               aGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                                                                       >  1:248541/1‑210 (MQ=255)
                                                                                                                                                                                               aGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                                                                       <  2:248541/210‑1 (MQ=255)
                                                                                                                                                                                                   gcgGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTgc                          <  1:227772/251‑1 (MQ=255)
                                                                                                                                                                                                         ttttcttttcGTGCTGGAAAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGCTGAAAAAGGGTTTATCGCGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGCAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTgcgggcgg                    >  1:101261/1‑251 (MQ=255)
                                                                                                                                                                                                              ttttcGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAACCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGCGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAAt                                                                                               >  1:40537/1‑171 (MQ=255)
                                                                                                                                                                                                              ttttcGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCAc                                                                            <  2:454408/190‑1 (MQ=255)
                                                                                                                                                                                                              ttttcGTGCTGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCAGGGACCGCATCCTCAAAGATCATCAAAt                                                                                               <  2:40537/171‑1 (MQ=255)
                                                                                                                                                                                                              ttttcGTGCTGGTAAGTCTTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCAc                                                                            >  1:454408/1‑190 (MQ=255)
                                                                                                                                                                                                                  cGTGCTGGCCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACATCCGTTCTCAATGGCTTCGACGCCGGATGAAAACGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGCTGATTTTGATTGCGGGCGGCACCGGGtt           >  1:339983/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGGTCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCACTGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGATGCCTATTGGCGCTTATGAACTCAACCTTTACGCGAAAGAAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCtt      >  1:515531/1‑251 (MQ=255)
                                                                                                                                                                                                                          tCAGTATTTGATGGTAGTGATGAATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATg   >  1:298162/1‑251 (MQ=255)
                                                                                                                                                                                                                           cAGTATTTGATGGTAGGGAGGAATGAGCGCGACAAACGGCCGTTCTCAATGGCGTCGACGCCGGATGAAAAAGGGTGTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTGACGCGAAAGCAGTCATGGACCGCATCCTCAAATATCATCAAATCGTGGTCGACATGCCCCACGGAGAAGCGTGGCGGCGCGATGATGAAGATCGGCCGATGATTTTGATTGCCTGGGGCACCGGTTTCTCTTATGc  <  2:339983/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
CCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAGCGTCCGATGATTTTGATTGCGGGCGGCACCGGGTTCTCTTATGC  >  NC_000913/4026410‑4026879

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