Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,895,605 A→G E267E (GAA→GAG pabB → aminodeoxychorismate synthase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,895,6050AG99.0% 376.7 / ‑4.6 97E267E (GAA→GAGpabBaminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (41/55);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

ATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAA  >  NC_000913/1895469‑1895739
                                                                                                                                        |                                                                                                                                      
aTTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATc                                                                                                                                     <  2:621497/140‑1 (MQ=255)
    tGGCGATGAATGGCAGGCATTCCCTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGa                                                                                                                                 <  2:1490463/140‑1 (MQ=255)
      gCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcc                                                                                                                               <  2:1301129/140‑1 (MQ=255)
      gCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcc                                                                                                                               >  1:1389193/1‑140 (MQ=255)
      gCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTCGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGAcc                                                                                                                               >  2:1171802/1‑140 (MQ=255)
       cGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAccc                                                                                                                              >  2:1445261/1‑140 (MQ=255)
       cGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAccc                                                                                                                              >  2:427971/1‑140 (MQ=255)
         aTGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgc                                                                                                                            <  2:531629/140‑1 (MQ=255)
         aTGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgc                                                                                                                            <  2:528696/140‑1 (MQ=255)
         aTGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgc                                                                                                                            >  2:233271/1‑140 (MQ=255)
          tGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgcc                                                                                                                           <  1:371417/140‑1 (MQ=255)
          tGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgcc                                                                                                                           <  1:1973892/140‑1 (MQ=255)
           gAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgccc                                                                                                                          >  2:417111/1‑140 (MQ=255)
           gAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGAcccgccc                                                                                                                          >  1:193029/1‑140 (MQ=255)
                gcaggcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTa                                                                                                                     >  1:648372/1‑140 (MQ=255)
                  aggcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTaaa                                                                                                                   <  1:1836140/140‑1 (MQ=255)
                   ggcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAg                                                                                                                  <  2:781975/140‑1 (MQ=255)
                    gcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAgg                                                                                                                 <  1:858959/140‑1 (MQ=255)
                    gcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAgg                                                                                                                 <  1:953432/140‑1 (MQ=255)
                    gcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAgg                                                                                                                 <  2:1389193/140‑1 (MQ=255)
                    gcaTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAgg                                                                                                                 <  1:996409/140‑1 (MQ=255)
                     caTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGc                                                                                                                <  1:815220/140‑1 (MQ=255)
                      aTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCa                                                                                                               <  2:1087677/140‑1 (MQ=255)
                      aTTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCa                                                                                                               <  1:1535514/140‑1 (MQ=255)
                       ttCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCAc                                                                                                              <  2:448988/140‑1 (MQ=255)
                        tCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCTGACCCGCCCGATTAAAGGCACg                                                                                                             >  2:1171141/1‑140 (MQ=255)
                          cTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCa                                                                                                               <  1:1040858/136‑1 (MQ=255)
                          cTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCa                                                                                                               >  2:1040858/1‑136 (MQ=255)
                              aGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCa                                                                                                       >  2:1723041/1‑140 (MQ=255)
                              aGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCa                                                                                                       >  2:1729394/1‑140 (MQ=255)
                                cTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACg                                                                                                     <  2:193029/140‑1 (MQ=255)
                                 ttAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGc                                                                                                    <  2:731943/140‑1 (MQ=255)
                                      cAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGcc                                                                                               <  1:527961/140‑1 (MQ=255)
                                         gCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGa                                                                                            >  1:1037305/1‑140 (MQ=255)
                                          ccAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGAt                                                                                           >  2:859092/1‑140 (MQ=255)
                                          ccAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGAt                                                                                           >  1:804723/1‑140 (MQ=255)
                                           cAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATc                                                                                          <  1:485812/140‑1 (MQ=255)
                                              ccGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTc                                                                                       >  1:1126359/1‑140 (MQ=255)
                                                   cgcCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCt                                                                                        <  1:7287/134‑1 (MQ=255)
                                                   cgcCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCt                                                                                        >  2:7287/1‑134 (MQ=255)
                                                   cgcCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGaa                                                                                  >  1:1887656/1‑140 (MQ=255)
                                                    gcCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAg                                                                                 <  2:337141/140‑1 (MQ=255)
                                                      cATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGAt                                                                               <  1:617294/140‑1 (MQ=255)
                                                              gcTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACa                                                                       <  2:1718903/140‑1 (MQ=255)
                                                              gcTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACa                                                                       <  2:1062349/140‑1 (MQ=255)
                                                              gcTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACa                                                                       <  1:427971/140‑1 (MQ=255)
                                                              gcTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACa                                                                       <  2:1027153/140‑1 (MQ=255)
                                                                ttttttACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAg                                                                     <  2:1025108/140‑1 (MQ=255)
                                                                  ttttACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGATTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATTCTCAGGAAGATAGCAAACAAGCa                                                                   >  1:1104801/1‑140 (MQ=255)
                                                                          cTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACt                                                           >  2:1265358/1‑140 (MQ=255)
                                                                            tGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACa                                                                       >  1:244565/1‑126 (MQ=255)
                                                                            tGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACa                                                                       <  2:244565/126‑1 (MQ=255)
                                                                            tGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACACCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTgg                                                         <  1:417111/140‑1 (MQ=255)
                                                                                cAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGaa                                                     <  2:1126359/140‑1 (MQ=255)
                                                                                  gggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACt                                                   <  2:104662/140‑1 (MQ=255)
                                                                                  gggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACt                                                   <  2:1887656/140‑1 (MQ=255)
                                                                                  gggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACt                                                   <  2:804723/140‑1 (MQ=255)
                                                                                  gggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACt                                                   <  1:1723041/140‑1 (MQ=255)
                                                                                  gggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACt                                                   <  2:1415425/140‑1 (MQ=255)
                                                                                  gggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACt                                                   <  1:1729394/140‑1 (MQ=255)
                                                                                   ggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTc                                                  >  2:319328/1‑140 (MQ=255)
                                                                                   ggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTc                                                  >  2:673602/1‑140 (MQ=255)
                                                                                   ggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTc                                                  >  1:1834113/1‑140 (MQ=255)
                                                                                   ggTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTc                                                  >  1:446397/1‑140 (MQ=255)
                                                                                    gTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCa                                                 >  2:1236267/1‑140 (MQ=255)
                                                                                     tGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAg                                                >  1:1531308/1‑140 (MQ=255)
                                                                                     tGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAg                                                >  2:281524/1‑140 (MQ=255)
                                                                                        aaTTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGa                                             >  1:735688/1‑140 (MQ=255)
                                                                                           tttAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAg                                          <  1:1918733/140‑1 (MQ=255)
                                                                                             tAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGAt                                        >  2:1159972/1‑140 (MQ=255)
                                                                                             tAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGAt                                        >  1:946477/1‑140 (MQ=255)
                                                                                               aGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCg                                      <  1:1159972/140‑1 (MQ=255)
                                                                                               aGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCg                                      <  1:319328/140‑1 (MQ=255)
                                                                                               aGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCg                                      >  2:1774483/1‑140 (MQ=255)
                                                                                                     tCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGa                                <  1:673602/140‑1 (MQ=255)
                                                                                                           gagCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCt                          >  2:413213/1‑140 (MQ=255)
                                                                                                            agCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGaa                               <  1:1239694/134‑1 (MQ=255)
                                                                                                            agCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGaa                               >  2:1239694/1‑134 (MQ=255)
                                                                                                            agCGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCtg                         >  2:758655/1‑140 (MQ=255)
                                                                                                              cGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAATTCAGCGAAAGATCGTGCCGAAAATCtgat                       >  1:634868/1‑140 (MQ=255)
                                                                                                              cGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAATTCAGCGAAAGATCGTGCCGAAAATCtgat                       <  2:634868/140‑1 (MQ=255)
                                                                                                              cGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCtgat                       >  2:159610/1‑140 (MQ=255)
                                                                                                              cGGTTTATTCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCtgat                       <  2:812216/140‑1 (MQ=255)
                                                                                                                     ttCTTTGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTc                >  2:1493263/1‑140 (MQ=255)
                                                                                                                         ttGTGATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGAtt            <  1:159610/140‑1 (MQ=255)
                                                                                                                          tgtgATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGAttt           >  1:1731830/1‑140 (MQ=255)
                                                                                                                          tgtgATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGAttt           <  1:1445261/140‑1 (MQ=255)
                                                                                                                           gtgATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTa          >  2:795724/1‑140 (MQ=255)
                                                                                                                           gtgATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTa          >  2:1552888/1‑140 (MQ=255)
                                                                                                                           gtgATAATAGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTa          <  2:1037305/140‑1 (MQ=255)
                                                                                                                              ataataGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATg       <  2:873235/140‑1 (MQ=255)
                                                                                                                              ataataGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATg       <  1:1265358/140‑1 (MQ=255)
                                                                                                                              ataataGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATg       <  2:1834113/140‑1 (MQ=255)
                                                                                                                              ataataGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATg       <  2:1531308/140‑1 (MQ=255)
                                                                                                                              ataataGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATg       <  1:758655/140‑1 (MQ=255)
                                                                                                                                  taGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTa   >  2:327557/1‑140 (MQ=255)
                                                                                                                                  taGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTa   <  1:696422/140‑1 (MQ=255)
                                                                                                                                   aGTGAGATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCTATTTAATGCGTaa  <  1:413213/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATTCTGGCGATGAATGGCAGGCATTCCTTCAGCTTAATCAGGCCAACCGCGCGCCATTTAGCGCTTTTTTACGTCTTGAACAGGGTGCAATTTTAAGCCTTTCGCCAGAGCGGTTTATTCTTTGTGATAATAGTGAAATCCAGACCCGCCCGATTAAAGGCACGCTACCACGCCTGCCCGATCCTCAGGAAGATAGCAAACAAGCAGTAAAACTGGCGAACTCAGCGAAAGATCGTGCCGAAAATCTGATGATTGTCGATTTAATGCGTAA  >  NC_000913/1895469‑1895739

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