Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,369,807 C→A 5.4% G182C (GGC→TGC)  nanK ← N‑acetylmannosamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,369,8070CA 5.4% 257.5 / 2.3 74G182C (GGC→TGC) nanKN‑acetylmannosamine kinase
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base A (2/2);  total (39/35)
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

TTTAATATCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTACACCGCCGCCAA  >  NC_000913/3369664‑3369953
                                                                                                                                               |                                                                                                                                                  
tttAATATCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGTAATGCCGCGAc                                                                                                                                               >  2:250069/1‑149 (MQ=255)
tttAATATCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGAc                                                                                                                                               >  4:513152/1‑149 (MQ=255)
  tAATATCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCa                                                                                                                                             <  4:553667/149‑1 (MQ=255)
     tatCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGaa                                                                                                                                          <  1:261528/149‑1 (MQ=255)
     tatCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGaa                                                                                                                                          <  3:135223/149‑1 (MQ=255)
     tatCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGaa                                                                                                                                          <  2:395050/149‑1 (MQ=255)
     tatCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCAGCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGaa                                                                                                                                          <  2:75711/149‑1 (MQ=255)
     tagcAGCGATCATCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGaa                                                                                                                                          >  2:90190/4‑149 (MQ=255)
       tCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGc                                                                                                                                        >  4:326674/1‑149 (MQ=255)
       tCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGc                                                                                                                                        >  2:477304/1‑149 (MQ=255)
         aGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCaa                                                                                                                                      >  2:241026/1‑149 (MQ=255)
         aGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCaa                                                                                                                                      >  3:48129/1‑149 (MQ=255)
                aGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTc                                                                                                                               <  3:513152/149‑1 (MQ=255)
                            gTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCtgt                                                                                                                   >  4:520745/1‑149 (MQ=255)
                            gTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGACAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCtgt                                                                                                                   <  4:48129/149‑1 (MQ=255)
                              aCGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGc                                                                                                                 <  3:326674/149‑1 (MQ=255)
                              aCGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGc                                                                                                                 <  4:446584/149‑1 (MQ=255)
                               cGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCg                                                                                                                <  1:241026/149‑1 (MQ=255)
                                 tGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTc                                                                                                              >  1:394121/1‑149 (MQ=255)
                                              aTCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGAc                                                                                                 >  2:1942/1‑149 (MQ=255)
                                               tCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACt                                                                                                >  3:434771/1‑149 (MQ=255)
                                               tCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACt                                                                                                >  2:307223/1‑149 (MQ=255)
                                                  gctgctGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTggg                                                                                             >  2:286923/1‑149 (MQ=255)
                                                    tgctgcGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGcc                                                                                           <  4:522908/149‑1 (MQ=255)
                                                      ctgcGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGCCCACAGCCGCAGACTGGGCCgt                                                                                         <  1:90190/149‑1 (MQ=255)
                                                            cTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGAt                                                                                   >  1:488304/1‑149 (MQ=255)
                                                            cTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGAt                                                                                   >  3:171030/1‑149 (MQ=255)
                                                            cTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGAt                                                                                   >  2:538752/1‑149 (MQ=255)
                                                             tGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCACAGACTGGGCCGTGTGGATc                                                                                  <  2:182038/149‑1 (MQ=255)
                                                                tCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGc                                                                               <  2:339455/149‑1 (MQ=255)
                                                                  gTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCaa                                                                             <  2:44959/149‑1 (MQ=255)
                                                                  gTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCaa                                                                             <  1:477304/149‑1 (MQ=255)
                                                                  gTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCaa                                                                             >  4:504266/1‑149 (MQ=255)
                                                                    cACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGc                                                                           <  4:270062/149‑1 (MQ=255)
                                                                       ccTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTa                                                                        >  4:580627/1‑149 (MQ=255)
                                                                       ccTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTa                                                                        >  3:329764/1‑149 (MQ=255)
                                                                        cTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGGGTCCAAAGCGGCAGACTGGTCCGTGTGGATCGGGAAGCGTAt                                                                       >  3:15676/1‑149 (MQ=255)
                                                                         tGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATg                                                                      <  4:202749/149‑1 (MQ=255)
                                                                           cccGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGcc                                                                    >  1:555819/1‑149 (MQ=255)
                                                                              ggCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGa                                                                 <  4:109524/149‑1 (MQ=255)
                                                                                   gcgTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGc                                                            <  1:538752/149‑1 (MQ=255)
                                                                                           aTAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAg                                                    >  3:283825/1‑149 (MQ=255)
                                                                                            tAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGa                                                   >  1:685168/1‑149 (MQ=255)
                                                                                               ttttCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCgtgt                                                                                       <  1:421984/110‑1 (MQ=255)
                                                                                               ttttCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCgtgt                                                                                       >  2:421984/1‑110 (MQ=255)
                                                                                               gtttCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCAGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGCGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCAAGACCg                                                >  4:76679/2‑149 (MQ=255)
                                                                                                 ttCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGcc                                              >  4:676943/1‑149 (MQ=255)
                                                                                                  tCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCa                                             <  4:329764/149‑1 (MQ=255)
                                                                                                    gCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCACTGCGGCGACCAGAAGCCATCGCGTCCACGCAACCAGTGTGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCTTATGCCCGATATGCCCAGCCAGACCGCCAgg                                           >  3:427449/1‑149 (MQ=255)
                                                                                                    gCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAgg                                           >  4:331021/1‑149 (MQ=255)
                                                                                                      aTCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGTAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGc                                         <  1:250069/149‑1 (MQ=255)
                                                                                                                gccaACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAg                               <  2:555819/149‑1 (MQ=255)
                                                                                                                 ccaACTCCCCCTGCGCTGCCGCTGCAATGCCGCGCCCAGAAGCAATCGCTTCGACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGAATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAga                              >  1:577917/1‑148 (MQ=255)
                                                                                                                  caACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGtt                             <  4:421725/149‑1 (MQ=255)
                                                                                                                  caACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGtt                             <  3:331021/149‑1 (MQ=255)
                                                                                                                  caACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGtt                             >  1:410177/1‑149 (MQ=255)
                                                                                                                  caACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGtt                             >  1:507100/1‑149 (MQ=255)
                                                                                                                  caACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGtt                             <  3:520745/149‑1 (MQ=255)
                                                                                                                   aaCTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGcc                                                      <  2:471076/123‑1 (MQ=255)
                                                                                                                   aaCTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGcc                                                      >  1:471076/1‑123 (MQ=255)
                                                                                                                   aaCTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGttt                            <  1:286923/149‑1 (MQ=255)
                                                                                                                        ccccTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGc                       >  2:655637/1‑149 (MQ=255)
                                                                                                                              cgctgccgctgcAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCa                 <  1:655637/149‑1 (MQ=255)
                                                                                                                                ctgccgctgcAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACt               >  3:524697/1‑149 (MQ=255)
                                                                                                                                 tgccgctgcAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTa              <  3:452291/149‑1 (MQ=255)
                                                                                                                                  gtcgcggcATTGCAGCGGCCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGttt                            >  2:661649/3‑134 (MQ=255)
                                                                                                                                  gtcgcggcATTGCAGCGGCCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGttt                            <  1:661649/132‑1 (MQ=255)
                                                                                                                                  gtcgcggcATTGCAGCGGCCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGttt                            <  3:634331/132‑1 (MQ=255)
                                                                                                                                  gtcgcggcATTGCAGCGGCCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGttt                            >  4:634331/7‑134 (MQ=255)
                                                                                                                                    cgctgcAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTacac           <  2:394121/149‑1 (MQ=255)
                                                                                                                                    cgctgcAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTacac           <  4:171030/149‑1 (MQ=255)
                                                                                                                                           aTGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTACACCgccgcc    >  1:511846/1‑149 (MQ=255)
                                                                                                                                           aTGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTACACCgccgcc    >  3:556327/1‑149 (MQ=255)
                                                                                                                                             gCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTACACCGCCGCCaa  >  2:258623/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TTTAATATCAGCGATCAGCCTTGCAAGCGTACGTGCGGAGCGGTGAATCAGCTGCTGCGCCTGCTCGTCACCCTGCCCGGCGCGCGTGAAAATAGTTTTCGCATCCGCGCCAGCCAACTCCCCCTGCGCTGCCGCTGCAATGCCGCGACCAGAAGCAATCGCTTCCACGCAACCTGTGCGTCCACAGCCGCAGACTGGGCCGTGTGGATCGGCAAGCGTATGCCCGATATGCCCCGCCAGACCGCCAGGGCCGGTAAGCAGTTTGCAGCCGCTCACTACACCGCCGCCAA  >  NC_000913/3369664‑3369953

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