Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,154,916 G→A 45.4% intergenic (‑69/‑85) argE ← / → argC acetylornithine deacetylase/N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,154,9160GA45.4% 30.1 / 133.4 99intergenic (‑69/‑85)argE/argCacetylornithine deacetylase/N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base A (19/26);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

CTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTAT  >  NC_000913/4154770‑4155060
                                                                                                                                                  |                                                                                                                                                
ctctTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGt                                                                                                                                                 >  1:910221/1‑148 (MQ=255)
  ctTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTg                                                                                                                                               <  2:421232/148‑1 (MQ=255)
  ctTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTg                                                                                                                                               >  2:35397/1‑148 (MQ=255)
  ctTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGa                                                                                                                                              <  4:470214/149‑1 (MQ=255)
  ctTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGa                                                                                                                                              >  4:198438/1‑149 (MQ=255)
   ttCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGa                                                                                                                                              >  1:774796/1‑148 (MQ=255)
   ttCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAg                                                                                                                                             >  3:358914/1‑149 (MQ=255)
    tCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAg                                                                                                                                             <  2:291149/148‑1 (MQ=255)
             cgtATTGAGGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAAGTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTAGCAATATTCATGCAGTATGTATGAATAAAAATACACTAACATTGAGCGTAATaaaa                                                                                                                                   <  4:15319/147‑1 (MQ=255)
             cTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATaaaa                                                                                                                                   <  3:842144/149‑1 (MQ=255)
             cTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATaaaa                                                                                                                                   <  3:112159/149‑1 (MQ=255)
             cTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATaaaa                                                                                                                                   >  4:129034/1‑149 (MQ=255)
             cTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATaaa                                                                                                                                    >  1:9781/1‑148 (MQ=255)
              ttATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAAc                                                                                                                                  >  3:601612/1‑149 (MQ=255)
                aTTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAAccc                                                                                                                                >  4:946303/1‑149 (MQ=255)
                     aGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAg                                                                                                                           <  3:198438/149‑1 (MQ=255)
                       gtgtGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGcc                                                                                                                         <  3:289943/149‑1 (MQ=255)
                           ggCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTaa                                                                                                                     <  3:582682/149‑1 (MQ=255)
                           ggCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTa                                                                                                                      <  1:643102/148‑1 (MQ=255)
                            gCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACAAGCCGTAAg                                                                                                                    >  3:204467/1‑149 (MQ=255)
                             cAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAg                                                                                                                    >  2:594491/1‑148 (MQ=255)
                             cAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAg                                                                                                                    >  1:170724/1‑148 (MQ=255)
                                tCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTg                                                                                                                 >  2:528588/1‑148 (MQ=255)
                                tCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTg                                                                                                                 >  1:221863/1‑148 (MQ=255)
                                tCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGAAGGTGa                                                                                                                <  2:48569/148‑1 (MQ=255)
                                 cAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGa                                                                                                                >  1:71628/1‑148 (MQ=255)
                                 cAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGa                                                                                                                >  1:596816/1‑148 (MQ=255)
                                  agagCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAAt                                                                                                              >  3:109702/1‑149 (MQ=255)
                                     gcgcGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGGTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGt                                                                                                            >  1:416445/1‑148 (MQ=255)
                                      cgcgGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGttt                                                                                                          <  3:95711/149‑1 (MQ=255)
                                      cgcgGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGtt                                                                                                           <  1:69282/148‑1 (MQ=255)
                                        cgGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGtttt                                                                                                         <  1:495777/148‑1 (MQ=255)
                                         ggTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTAc                                                                                                       <  3:668243/149‑1 (MQ=255)
                                          gTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACg                                                                                                      >  4:414915/1‑149 (MQ=255)
                                           tAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGt                                                                                                     >  3:98170/1‑149 (MQ=255)
                                                 tcGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGc                                                                                                                                            <  3:230411/104‑1 (MQ=255)
                                                 tcGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGc                                                                                                                                            >  4:230411/1‑104 (MQ=255)
                                                 tcGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAAcc                                                                                               >  4:33431/1‑149 (MQ=255)
                                                  cGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCt                                                                                              >  4:493391/1‑149 (MQ=255)
                                                    aTAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGtttt                                                                                                         >  2:122800/1‑136 (MQ=255)
                                                    aTAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGtttt                                                                                                         <  1:122800/136‑1 (MQ=255)
                                                              gTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGAc                                                                                   >  2:677209/1‑148 (MQ=255)
                                                                 aTTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaa                                                                               >  3:480209/1‑149 (MQ=255)
                                                                  tttGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaag                                                                              >  3:904236/1‑149 (MQ=255)
                                                                   ttGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaag                                                                              >  2:231508/1‑148 (MQ=255)
                                                                     gTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaagaag                                                                           <  4:371717/149‑1 (MQ=255)
                                                                      tttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaagaag                                                                           <  2:573579/148‑1 (MQ=255)
                                                                       ttttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCGACCAGACATAAGAAGGt                                                                         <  4:358914/149‑1 (MQ=255)
                                                                         ttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTg                                                                        >  1:341620/1‑148 (MQ=255)
                                                                         ttCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAACACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTg                                                                        >  1:888409/1‑148 (MQ=255)
                                                                           cATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGaa                                                                      <  1:594491/148‑1 (MQ=255)
                                                                              tgttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAg                                                                   >  1:402372/1‑148 (MQ=255)
                                                                              tgttgACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAg                                                                   <  2:7438/148‑1 (MQ=255)
                                                                                   acacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCg                                                              >  1:165227/1‑148 (MQ=255)
                                                                                   acacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGa                                                             <  4:177927/149‑1 (MQ=255)
                                                                                    cacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGa                                                             <  2:774796/148‑1 (MQ=255)
                                                                                    cacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGa                                                             <  2:416445/148‑1 (MQ=255)
                                                                                    cacacCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGa                                                             <  2:966836/148‑1 (MQ=255)
                                                                                       acCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAgg                                                                                                                   >  4:112072/1‑91 (MQ=255)
                                                                                       acCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAgg                                                                                                                   <  3:112072/91‑1 (MQ=255)
                                                                                          tctGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGa                                                       <  1:256124/148‑1 (MQ=255)
                                                                                            tGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATa                                                    <  3:384873/149‑1 (MQ=255)
                                                                                            tGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATa                                                    <  4:601612/149‑1 (MQ=255)
                                                                                              gTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACg                                                  <  4:399142/149‑1 (MQ=255)
                                                                                              gTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACg                                                  <  3:493391/149‑1 (MQ=255)
                                                                                               tCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACg                                                  >  1:674673/1‑148 (MQ=255)
                                                                                                cATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGc                                                 >  2:221833/1‑148 (MQ=255)
                                                                                                   gATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGaa                                                      >  4:384231/1‑140 (MQ=255)
                                                                                                   gATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGaa                                                      <  3:384231/140‑1 (MQ=255)
                                                                                                     tAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGAtt                                            >  2:716826/1‑148 (MQ=255)
                                                                                                       gTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATtgt                                          <  1:325856/148‑1 (MQ=255)
                                                                                                         aTCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTgg                                        <  1:528588/148‑1 (MQ=255)
                                                                                                         aTCAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTggg                                       >  4:951361/1‑149 (MQ=255)
                                                                                                          tCAGTATTCATGCAGAATTTATGAATAAAAATACACTAACATTGAGCGTAATAAATCCCCCCAGCCCTAAGATGACTGTTTTACGTTTAACCTGTCACCCAGACATGGCAAGGTGAATAGACCCGATGTTGAATACGCGGATTGTGGGt                                      <  4:842444/149‑1 (MQ=255)
                                                                                                          tCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTggg                                       <  2:910221/148‑1 (MQ=255)
                                                                                                          tCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTggg                                       <  2:402372/148‑1 (MQ=255)
                                                                                                           cAATATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGt                                      >  1:238553/1‑148 (MQ=255)
                                                                                                            aaTATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAAc                                                                                                <  4:402430/89‑1 (MQ=255)
                                                                                                            aaTATTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAAc                                                                                                >  3:402430/1‑89 (MQ=255)
                                                                                                              tatTCATGCAGTATTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGcc                                   >  2:316095/1‑148 (MQ=255)
                                                                                                                ttCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAg                                 >  2:53261/1‑148 (MQ=255)
                                                                                                                          aTTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATtgt                                          <  2:723205/129‑1 (MQ=255)
                                                                                                                          aTTTATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATtgt                                          >  1:723205/1‑129 (MQ=255)
                                                                                                                           tttATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTg                      <  1:677209/148‑1 (MQ=255)
                                                                                                                            ttATGAATAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTgg                     <  2:159423/148‑1 (MQ=255)
                                                                                                                                 aaTAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCag                <  2:221863/148‑1 (MQ=255)
                                                                                                                                 aaTAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCaga               <  3:951361/149‑1 (MQ=255)
                                                                                                                                  aTAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCt                                                                                              <  2:315892/69‑1 (MQ=255)
                                                                                                                                  aTAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCt                                                                                              >  1:315892/1‑69 (MQ=255)
                                                                                                                                  aTAAAAATACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCaga               >  1:24214/1‑148 (MQ=255)
                                                                                                                                   tAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCagag              >  2:541728/1‑148 (MQ=255)
                                                                                                                                       aaTACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAg          <  1:22609/148‑1 (MQ=255)
                                                                                                                                       aaTACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAg          >  2:923298/1‑148 (MQ=255)
                                                                                                                                       aaTACACTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGt         <  4:204467/149‑1 (MQ=255)
                                                                                                                                            acTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGAcc     <  2:238553/148‑1 (MQ=255)
                                                                                                                                             cTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGa                                                             <  1:554386/91‑1 (MQ=255)
                                                                                                                                             cTAACATTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGa                                                             >  2:554386/1‑91 (MQ=255)
                                                                                                                                              tAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTAt  <  4:40565/149‑1 (MQ=255)
                                                                                                                                              tAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTAt  >  4:457805/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTAT  >  NC_000913/4154770‑4155060

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