Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,717,169 T→G 7.4% L139R (CTC→CGC) ‡ ung → uracil‑DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,717,1690TG 7.4% 243.5 / 11.1 81L139R (CTC→CGC) ‡unguracil‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (38/37);  new base G (2/4);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

TGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCATCCGT  >  NC_000913/2717022‑2717318
                                                                                                                                                   |                                                                                                                                                     
tGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCg                                                                                                                                                     <  2:1264318/150‑1 (MQ=255)
      tGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACACCAGCCTCGGCTggg                                                                                                                                               <  2:1617303/150‑1 (MQ=255)
       gTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGa                                                                                                                                              >  1:1312555/1‑150 (MQ=255)
         ataAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtg                                                                                                                                            <  2:907476/150‑5 (MQ=255)
          taAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtgg                                                                                                                                           <  1:1644453/150‑6 (MQ=255)
           aaaGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGAGGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtggg                                                                                                                                          <  1:1829155/150‑7 (MQ=255)
           aaaGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAcc                                                                                                                                          <  1:1198270/150‑1 (MQ=255)
           aaaGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAAcc                                                                                                                                          <  1:1436317/150‑1 (MQ=255)
           aaaGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCGAGGCTggcgtggg                                                                                                                                          <  2:1207002/150‑7 (MQ=255)
                 cTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCAcc                                                                                                                                    >  2:142580/1‑150 (MQ=255)
                 cTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCAcc                                                                                                                                    >  1:1827936/1‑150 (MQ=255)
                   ggAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGa                                                                                                                                  >  1:35607/1‑150 (MQ=255)
                    gAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGAt                                                                                                                                 >  1:367751/1‑150 (MQ=255)
                         tACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGt                                                                                                                            <  1:1783666/150‑1 (MQ=255)
                          aCTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTg                                                                                                                           <  2:35607/150‑1 (MQ=255)
                          aCTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTg                                                                                                                           <  2:1900455/150‑1 (MQ=255)
                           cTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGa                                                                                                                          >  1:273400/1‑150 (MQ=255)
                                ccGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGc                                                                                                                     >  1:1170041/1‑150 (MQ=255)
                                   ggCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCCCGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTg                                                                                                                  >  1:1825768/1‑150 (MQ=255)
                                    gCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAACCCTTCACCGATAAAGTGATCAGCCTGa                                                                                                                 >  2:1366491/1‑150 (MQ=255)
                                     cTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGAt                                                                                                                >  1:144027/1‑150 (MQ=255)
                                     cTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGAt                                                                                                                <  1:831613/150‑1 (MQ=255)
                                        cACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTaa                                                                                                             >  1:1176149/1‑150 (MQ=255)
                                          cccGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAcc                                                                                                           <  1:1651361/150‑1 (MQ=255)
                                            cGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAg                                                                                                         <  2:25838/150‑1 (MQ=255)
                                             gCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGc                                                                                                        >  1:1786697/1‑150 (MQ=255)
                                             gCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGc                                                                                                        >  2:1041012/1‑150 (MQ=255)
                                              cccTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCa                                                                                                       <  2:1825768/150‑1 (MQ=255)
                                                  aaTCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcgc                                                                                                   <  1:997639/150‑1 (MQ=255)
                                                   aTCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATcgcg                                                                                                  >  2:556094/1‑150 (MQ=255)
                                                       tGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAgg                                                                                              >  2:1231982/1‑150 (MQ=255)
                                                         gTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCg                                                                                            >  1:1087314/1‑150 (MQ=255)
                                                         gTTATCTTGAAAGCTCGGCGCGTCAGGGCGTTATGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCACCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAAACAGCATCGCGAAcgcg                                                                                            >  1:1760352/1‑150 (MQ=255)
                                                          ttATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgt                                                                                           <  2:1176149/150‑1 (MQ=255)
                                                              cTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCgtggtg                                                                                       <  2:699737/150‑1 (MQ=255)
                                                               ttGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGgtgt                                                                                      >  1:1634650/1‑150 (MQ=255)
                                                                    aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttg                                                                                 <  2:1827936/150‑1 (MQ=255)
                                                                    aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttg                                                                                 <  2:667709/150‑1 (MQ=255)
                                                                    aGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttg                                                                                 <  1:592154/150‑1 (MQ=255)
                                                                     gCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttgt                                                                                >  2:1953439/1‑150 (MQ=255)
                                                                       tGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTttgttg                                                                              <  2:1918926/150‑1 (MQ=255)
                                                                           cgcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTggg                                                                          <  1:685040/150‑1 (MQ=255)
                                                                           cgcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTggg                                                                          <  2:1485388/150‑1 (MQ=255)
                                                                            gcgTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTgggg                                                                         >  2:1999184/1‑150 (MQ=255)
                                                                              gTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGAt                                                                       <  2:367751/150‑1 (MQ=255)
                                                                               tCATGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATc                                                                      <  1:405089/150‑1 (MQ=255)
                                                                                cAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCg                                                                     >  2:406719/1‑150 (MQ=255)
                                                                                  gggCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCa                                                                   >  2:1152533/1‑150 (MQ=255)
                                                                                  gggCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCa                                                                   <  2:634431/150‑1 (MQ=255)
                                                                                   ggCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCAt                                                                  >  2:340602/1‑150 (MQ=255)
                                                                                    gCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATg                                                                 <  2:1414970/150‑1 (MQ=255)
                                                                                      gtacTGCTACTCAAAACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATgcg                                                               >  2:1766197/4‑150 (MQ=255)
                                                                                           gCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCaaaa                                                          >  1:1548215/1‑150 (MQ=255)
                                                                                           aatACTCAAAACTGTGTTGACGGTACGCGCAGGTCAGACGAATTCCAACGCCAGCATCGGATGAGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCaaaa                                                          >  2:216960/3‑150 (MQ=255)
                                                                                                  aaTACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAggg                                                   <  1:406719/150‑1 (MQ=255)
                                                                                                     aCTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCg                                                >  1:1578527/1‑150 (MQ=255)
                                                                                                     aCTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCg                                                >  1:891297/1‑150 (MQ=255)
                                                                                                     aCTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTCTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCg                                                >  2:472961/1‑150 (MQ=255)
                                                                                                          gtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGttttt                                                                                  <  2:1222622/111‑1 (MQ=255)
                                                                                                          gtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTggg                                                                          >  1:1895315/1‑119 (MQ=255)
                                                                                                          gtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTggg                                                                          <  2:1895315/119‑1 (MQ=255)
                                                                                                          gtTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTAACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATtat                                           >  2:921912/1‑150 (MQ=255)
                                                                                                          attGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGttttt                                                                                  >  1:1222622/2‑111 (MQ=255)
                                                                                                           ttGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATtata                                          <  1:1319587/150‑1 (MQ=255)
                                                                                                            tGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAg                                         <  1:1766197/150‑1 (MQ=255)
                                                                                                             gACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTCTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGa                                        <  1:472961/150‑1 (MQ=255)
                                                                                                              aCGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGAt                                       >  1:505909/1‑150 (MQ=255)
                                                                                                               cGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATa                                      <  1:1327794/150‑1 (MQ=255)
                                                                                                                 gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAg                                    <  2:1634650/150‑1 (MQ=255)
                                                                                                                 gTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAg                                    <  1:1041012/150‑1 (MQ=255)
                                                                                                                  tACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGc                                   >  1:1609863/1‑150 (MQ=255)
                                                                                                                  tACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGc                                   <  2:144027/150‑1 (MQ=255)
                                                                                                                                gcgcATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTa                     >  1:803686/1‑150 (MQ=255)
                                                                                                                                     ttCCCACGCCCGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCAGGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGaa                >  1:501608/1‑150 (MQ=255)
                                                                                                                                     ttCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGTGGCGATTATAGATAAGCAACGCCATCATGTACTGaa                <  2:1312555/150‑1 (MQ=255)
                                                                                                                                      tCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATCTACTGaaa               <  1:353070/150‑1 (MQ=255)
                                                                                                                                       cccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAg              >  2:128630/1‑150 (MQ=255)
                                                                                                                                        gtttcccAGCCGAGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGc             >  2:609446/6‑150 (MQ=255)
                                                                                                                                        ccACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGc             >  1:1554798/1‑150 (MQ=255)
                                                                                                                                         tttcccAGCCGAGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCa            >  2:1984734/5‑150 (MQ=255)
                                                                                                                                                  cTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCATCCg   <  2:1311215/150‑1 (MQ=255)
                                                                                                                                                   tCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCATCCGt  >  1:882370/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGAATATGTATAAAGAGCTGGAAAATACTATTCCGGGCTTCACCCGCCCTAATCATGGTTATCTTGAAAGCTGGGCGCGTCAGGGCGTTCTGCTACTCAATACTGTGTTGACGGTACGCGCAGGTCAGGCGCATTCCCACGCCAGCCTCGGCTGGGAAACCTTCACCGATAAAGTGATCAGCCTGATTAACCAGCATCGCGAAGGCGTGGTGTTTTTGTTGTGGGGATCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCATCCGT  >  NC_000913/2717022‑2717318

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