Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,819,816 A→C 5.2% L452R (CTT→CGT)  ligB ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,819,8160AC 5.2% 262.4 / 3.4 86L452R (CTT→CGT) ligBDNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base C (2/3);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGT  >  NC_000913/3819669‑3819963
                                                                                                                                                   |                                                                                                                                                   
aGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATaaaa                                                                                                                                                     <  6:190914/148‑1 (MQ=255)
aGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATaaaa                                                                                                                                                     <  4:100399/148‑1 (MQ=255)
 gcgcCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGc                                                                                                                                                   <  5:41683/149‑1 (MQ=255)
    ccgccCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAg                                                                                                                                                <  5:218112/149‑1 (MQ=255)
    ccgccCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAg                                                                                                                                                <  6:52478/149‑1 (MQ=255)
         cGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgagcaa                                                                                                                                           <  7:59447/149‑6 (MQ=255)
            gTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGat                                                                                                                                        <  7:146854/149‑1 (MQ=255)
            gTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGa                                                                                                                                         >  7:293476/1‑148 (MQ=255)
                gCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGt                                                                                                                                    <  8:293476/149‑1 (MQ=255)
                      ttCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAg                                                                                                                              <  8:32001/149‑1 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATaaaa                                                                                                                                                     >  1:53835/1‑123 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATaaaa                                                                                                                                                     <  2:53835/123‑1 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCa                                                                                                                                                 <  8:294360/127‑1 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCa                                                                                                                                                 >  7:294360/1‑127 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTc                                                                                                                                  <  6:189325/142‑1 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTc                                                                                                                                  >  5:189325/1‑142 (MQ=255)
                         ccATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGa                                                                                                                           >  4:370093/1‑149 (MQ=255)
                           aTTGCCATCACCCAGCGAGTAACAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATg                                                                                                                         <  6:81007/149‑1 (MQ=255)
                                cATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCt                                                                                                                    >  2:124970/1‑149 (MQ=255)
                                  tCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGa                                                                                                                  <  3:206675/149‑1 (MQ=255)
                                    aCCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATg                                                                                                                >  5:365853/1‑149 (MQ=255)
                                      ccAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCa                                                                                                              <  5:368299/149‑1 (MQ=255)
                                       catcGAGTAAAAGTCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCTATCCACGGCGTTTTCTGTAATTGCTCTGGCGTTAATAAAATCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAg                                                                                                             <  5:339212/146‑1 (MQ=255)
                                       cAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAg                                                                                                             >  7:248178/1‑149 (MQ=255)
                                                aaaGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCa                                                                                                              >  4:360284/1‑139 (MQ=255)
                                                aaaGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCa                                                                                                              <  3:360284/139‑1 (MQ=255)
                                                 aaGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgagcaa                                                                                                                                           >  8:218864/1‑104 (MQ=255)
                                                 aaGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCTTTTATTAACGCCAgagcaa                                                                                                                                           <  7:218864/109‑6 (MQ=255)
                                                      agctTACGAGCAAGATTAAACTGATTACATAGCTCAGCACTTTTACATATCGCGATACACCTACTGATCTGTAATAGAACTGGCGTTAATAAAAGACAGGAAAAGATATGTTCAAAGCGATCAGTCTGAAGCAGCTCGCGACACCCGGc                                                                                              >  6:184754/2‑149 (MQ=255)
                                                        ctTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGc                                                                                              <  6:366457/147‑1 (MQ=255)
                                                        ctTACGAGCCAGATTAAACTATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGc                                                                                              <  7:267142/146‑1 (MQ=255)
                                                        ctTACGAGCCAGATTAAACTATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGc                                                                                              >  8:267142/1‑146 (MQ=255)
                                                         ttACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCt                                                                                            >  4:81347/1‑148 (MQ=255)
                                                               gCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAAt                                                                                     >  8:175988/1‑149 (MQ=255)
                                                                 cAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGc                                                                                   <  8:248178/149‑1 (MQ=255)
                                                                   gATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCa                                                                                 >  5:269566/1‑149 (MQ=255)
                                                                     ttAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAATCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATc                                                                               >  7:190383/1‑149 (MQ=255)
                                                                      tAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATc                                                                               >  8:235411/1‑148 (MQ=255)
                                                                       aaaCTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATc                                                                               >  8:374237/1‑147 (MQ=255)
                                                                             aTGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGccc                                                                         >  5:389358/1‑147 (MQ=255)
                                                                             aTGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCa                                                                        <  8:413103/148‑1 (MQ=255)
                                                                              tGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGccc                                                                         >  8:326760/1‑146 (MQ=255)
                                                                               gCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAg                                                                      >  1:311308/1‑148 (MQ=255)
                                                                                 cATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAAc                                                                   >  2:193726/1‑149 (MQ=255)
                                                                                  aTAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTCCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAAcc                                                                  >  4:321109/1‑149 (MQ=255)
                                                                                   tAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGACATCCAGCCCAAGAACCt                                                                 >  3:19937/1‑149 (MQ=255)
                                                                                     gCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGt                                                               >  1:389863/1‑149 (MQ=255)
                                                                                     gCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTCCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGt                                                               <  3:321109/149‑1 (MQ=255)
                                                                                       tGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAATCCAGGAAAAGATATGTTCAAATCGATGAGTCTGATGCAGCGCGCGCCAACCTGCCTCACAAATGCCATCCAGCCCAAGAAACTGttt                                                             >  6:175551/1‑149 (MQ=255)
                                                                                         cgcACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGttt                                                             <  6:269566/147‑1 (MQ=255)
                                                                                                ttacttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGccc                                                        >  7:68792/1‑145 (MQ=255)
                                                                                                ttacttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGccc                                                        <  8:68792/145‑1 (MQ=255)
                                                                                                    ttttCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCccagc                                                   >  3:288519/1‑146 (MQ=255)
                                                                                                        agcgATCCCAGGAGTGTTCTGTAATGTCTCTGGCGTTAATAAAAGCCAGGAAAAGAAATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCACAGCCAGACt                                             <  5:175551/147‑1 (MQ=255)
                                                                                                         gcgATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAACAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCACCCGGCCTCACCAATGCCAACCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACt                                             >  3:56214/1‑147 (MQ=255)
                                                                                                         gcgATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACt                                             >  1:218263/1‑147 (MQ=255)
                                                                                                         gcgATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTaa                                           >  7:217826/1‑149 (MQ=255)
                                                                                                               cccGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTg                                      <  1:193726/148‑1 (MQ=255)
                                                                                                               cccGGCCTGTTCCGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCGAGCCAGACTAAGCGTg                                      <  6:336874/148‑1 (MQ=255)
                                                                                                                      tgtTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAAc                              <  2:311308/149‑1 (MQ=255)
                                                                                                                          cTGTAATTGCTCTGGCGTTAATAAACGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACAAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGtt                          >  5:205184/1‑149 (MQ=255)
                                                                                                                           tGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCAc                                                                                         >  6:173779/1‑85 (MQ=255)
                                                                                                                           tGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCAc                                                                                         >  6:179683/1‑85 (MQ=255)
                                                                                                                           tGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCAc                                                                                         <  5:179683/85‑1 (MQ=255)
                                                                                                                           tGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCAc                                                                                         <  5:173779/85‑1 (MQ=255)
                                                                                                                            gTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCCACCGGCCTCCCCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTcc                        >  6:163597/1‑149 (MQ=255)
                                                                                                                             tAATTTCTCTGGCGTTAATAAACGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACt                                             <  7:290404/127‑1 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCCCCAATGCCATCCAGCCCAAGAACCTCATATGCAACCAGCCAGACt                                             >  8:290404/1‑127 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTc                         <  7:235411/147‑1 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTc                         <  5:142944/147‑1 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTc                         <  7:326760/147‑1 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTc                         <  7:374237/147‑1 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCt                       <  8:190383/149‑1 (MQ=255)
                                                                                                                             tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCt                       <  2:218263/149‑1 (MQ=255)
                                                                                                                               aTTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGa                     >  1:318880/1‑149 (MQ=255)
                                                                                                                                ttGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCGGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCccagcca                                                 <  7:281913/120‑1 (MQ=255)
                                                                                                                                ttGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCCGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCccagcca                                                 >  8:281913/1‑120 (MQ=255)
                                                                                                                                ttGATCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGAc                    <  5:184754/149‑1 (MQ=255)
                                                                                                                                        gCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATc             <  3:81347/148‑1 (MQ=255)
                                                                                                                                        gCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTACTTACAAACATCa            <  8:217826/149‑1 (MQ=255)
                                                                                                                                        gCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCCATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATCCCATCCAGCCCAAGAACCTGTTATGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATc             <  1:124970/148‑1 (MQ=255)
                                                                                                                                            taataaAAGCCAGTAAAAGATCTGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTTTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGaa         <  6:205184/148‑1 (MQ=255)
                                                                                                                                             aataaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCccagc                                                   >  5:123152/1‑105 (MQ=255)
                                                                                                                                             aataaAAGACAGGAAAAGATATGTTCAAAGAGATGAGTCTGATGCAGCTCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCccagc                                                   <  6:123152/105‑1 (MQ=255)
                                                                                                                                              ataaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAAt                                                                                     <  7:367296/70‑1 (MQ=255)
                                                                                                                                              ataaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAAt                                                                                     >  8:367296/1‑70 (MQ=255)
                                                                                                                                              ataaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTc                         >  2:413713/1‑130 (MQ=255)
                                                                                                                                              ataaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTc                         <  1:413713/130‑1 (MQ=255)
                                                                                                                                               taaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGc                                                          >  5:224728/1‑96 (MQ=255)
                                                                                                                                               taaAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGc                                                          <  6:224728/96‑1 (MQ=255)
                                                                                                                                                 aaaGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCTCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAg   <  6:389358/149‑1 (MQ=255)
                                                                                                                                                  aaGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGt  <  2:318880/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCGCCGCCCGGGTTAGCGGTATTCCCATTGCCATCACCCAGCGAGTAAAAGGCTGCTTACGAGCCAGATTAAACTGATGCCATAGCTGCGCACTTTTACTTTTCGCGATCCCCGGCGTGTTCTGTAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGT  >  NC_000913/3819669‑3819963

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