Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,335,681 T→C 5.1% intergenic (+223/‑150) yciX → / → acnA uncharacterized protein/aconitate hydratase 1; aconitase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,6810TC 5.1% 369.6 / 6.8 117intergenic (+223/‑150)yciX/acnAuncharacterized protein/aconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base T (53/58);  new base C (3/3);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 117 total aligned reads displayed.

TCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGA  >  NC_000913/1335535‑1335827
                                                                                                                                                  |                                                                                                                                                  
   cAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggc                                                                                                                                               >  8:300333/1‑149 (MQ=255)
    aaCACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                                                                              >  4:426374/1‑149 (MQ=255)
    aaCACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                                                                              >  6:450173/1‑149 (MQ=255)
    aaCACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCAGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                                                                              >  3:252676/1‑149 (MQ=255)
    aaCACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGGTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                                                                              <  3:68358/149‑1 (MQ=255)
     acacAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGcggcgg                                                                                                                                             <  8:334930/149‑1 (MQ=255)
     acacAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATAGGTCCTGGATGGGCGCGTTTTCTCTCTCTATGCTCTTCTGGTCCGAACTCCTCTGGATCATTAATCACTGATCCTTGACACCGCTGcggcgg                                                                                                                                             >  1:228855/1‑149 (MQ=255)
     acaaAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCTCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGcggcgg                                                                                                                                             <  2:225962/149‑1 (MQ=255)
       acaAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCgggg                                                                                                                                           >  6:354219/1‑149 (MQ=255)
       acaAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCgggg                                                                                                                                           <  7:305453/149‑1 (MQ=255)
         aaaaTAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGtt                                                                                                                                         <  7:115640/149‑1 (MQ=255)
         aaaaTAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGtt                                                                                                                                         <  4:127140/149‑1 (MQ=255)
         aaaaTAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGtt                                                                                                                                         >  1:68026/1‑149 (MQ=255)
          aaaTAATAACCGTTAAGGGTGTAGCCTATGATCAACACATATATTAAATATTTGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTg                                                                                                                                        <  5:136206/149‑1 (MQ=255)
           aataatAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGCAGCGGGGtcaa                                                                                                                                       <  4:228125/149‑4 (MQ=255)
                     gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgcttt                                                                                                                             <  2:68026/149‑1 (MQ=255)
                         aGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTATCTGTCGATGCTATTCTGGGCAGAACTCCTCTGGATCATTACTAACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCca                                                                                                                         >  2:369762/1‑149 (MQ=255)
                            gtgtAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACaa                                                                                                                      <  2:61140/149‑1 (MQ=255)
                              gtAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAgg                                                                                                                    <  8:457806/149‑1 (MQ=255)
                                 gCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGttt                                                                                                                 <  5:450173/149‑1 (MQ=255)
                                  ccTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtc                                                                                                                <  2:456647/149‑1 (MQ=255)
                                   cTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtct                                                                                                               >  3:358881/1‑149 (MQ=255)
                                     cgagcaaaCACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtctc                                                                                                              <  1:384840/142‑1 (MQ=255)
                                        aTCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTctcctc                                                                                                           <  8:343327/148‑1 (MQ=255)
                                         tCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCtt                                                                                                         <  5:338087/149‑1 (MQ=255)
                                         tCAACACAAATATGAAATATTGGTCCTGGAGGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCtt                                                                                                         <  2:89492/149‑1 (MQ=255)
                                           aaCACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGttt                                                                                                                 <  6:241297/139‑1 (MQ=255)
                                            acacAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGCTCCTTGACCCCGCTGCGGCGGGGGTGTCAGTTGGTTTGCCCCAAGGTTTgtccggatcg                                                                                                       >  2:384840/1‑142 (MQ=255)
                                            acacAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTa                                                                                                      >  2:52265/1‑149 (MQ=255)
                                               caAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCACTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTc                                                                                                                                      >  1:97716/1‑114 (MQ=255)
                                                aaaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCaa                                                                                                  >  2:122637/1‑149 (MQ=255)
                                                   tatGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAttt                                                                                               >  4:93795/1‑149 (MQ=255)
                                                     tGAATTATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGTGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTgg                                                                                             <  2:228855/149‑1 (MQ=255)
                                                      gAAATATTGGTCATGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTggg                                                                                            >  6:388190/1‑149 (MQ=255)
                                                              ggTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgcttt                                                                                                                             >  2:26438/1‑108 (MQ=255)
                                                              ggTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTAt                                                                                     <  4:252676/148‑1 (MQ=255)
                                                               gTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATc                                                                                    <  4:349508/148‑1 (MQ=255)
                                                                tCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCaa                                                                                                  <  1:282953/133‑1 (MQ=255)
                                                                tCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCaa                                                                                                  >  2:282953/1‑133 (MQ=255)
                                                                    ggATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCg                                                                              >  6:258812/1‑149 (MQ=255)
                                                                     gATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCACGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATc                                                                               <  5:226525/147‑1 (MQ=255)
                                                                            gcgTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCTGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGa                                                                      >  4:6424/1‑149 (MQ=255)
                                                                                      tGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt                                                             >  5:30573/1‑148 (MQ=255)
                                                                                        tCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGtt                                                            <  5:258812/147‑1 (MQ=255)
                                                                                        tCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa                                                           <  8:244350/148‑1 (MQ=255)
                                                                                            tGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCttt                                                      <  8:224008/149‑1 (MQ=255)
                                                                                              ctctTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                    <  3:335623/149‑1 (MQ=255)
                                                                                              ctctTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTa                                                     <  1:122637/148‑1 (MQ=255)
                                                                                               tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                    >  8:303689/1‑148 (MQ=255)
                                                                                               tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTa                                                     >  4:25160/1‑147 (MQ=255)
                                                                                               tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTa                                                     <  3:25160/147‑1 (MQ=255)
                                                                                                 ttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAat                                                   <  2:391606/147‑1 (MQ=255)
                                                                                                      ggCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCa                                              <  4:358881/147‑1 (MQ=255)
                                                                                                         cGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCa                                              >  8:381764/1‑144 (MQ=255)
                                                                                                         cGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCa                                              <  7:381764/144‑1 (MQ=255)
                                                                                                           aaCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                       <  5:354219/149‑1 (MQ=255)
                                                                                                           aaCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                       <  7:64768/149‑1 (MQ=255)
                                                                                                           aaCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                       <  3:6424/149‑1 (MQ=255)
                                                                                                           aaCTACTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTGTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                       <  3:78864/149‑1 (MQ=255)
                                                                                                              tcctcTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                       >  5:378767/1‑146 (MQ=255)
                                                                                                                ctctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCa                                                                                   >  1:144482/1‑100 (MQ=255)
                                                                                                                ctctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCa                                                                                   <  2:144482/100‑1 (MQ=255)
                                                                                                                ctctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGaa                                  >  1:463530/1‑149 (MQ=255)
                                                                                                                ctctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATATTTAATATTCACACTGAAgaga                                   >  4:134128/1‑148 (MQ=255)
                                                                                                                 tctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGaa                                  >  5:289349/1‑148 (MQ=255)
                                                                                                                  ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAAt                                 >  7:37905/1‑148 (MQ=255)
                                                                                                                   tGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCAccc                                           <  7:330329/137‑1 (MQ=255)
                                                                                                                   tGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAAt                                 >  7:120028/1‑147 (MQ=255)
                                                                                                                     gATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAgg                             <  7:303689/149‑1 (MQ=255)
                                                                                                                     gATCATTACTAAATGATCCTTTACCCCGCTTGGGCGGGGTTGTAATTTGCTTTGCCACAATGTTTCTCCTCTTTTATCAAATTGGGGTTTTATCAAATCGTTACGCGATGTTTTTGTTATCTTTAATATTCACCCTGAAGAGAATCAg                              <  3:134128/148‑1 (MQ=255)
                                                                                                                          ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATAtt                                                <  1:370957/125‑1 (MQ=255)
                                                                                                                          ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATAtt                                                >  2:370957/1‑125 (MQ=255)
                                                                                                                          ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTc                        >  4:45737/1‑149 (MQ=255)
                                                                                                                          ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATAAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTc                        <  3:93795/149‑1 (MQ=255)
                                                                                                                           tACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAGGAGAATCAGGGCTTCg                       >  1:365355/1‑149 (MQ=255)
                                                                                                                            aCTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCg                       >  2:453901/1‑148 (MQ=255)
                                                                                                                            aCTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGGTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTc                        >  1:375924/1‑147 (MQ=255)
                                                                                                                            aCTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGCTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTc                        <  2:375924/147‑1 (MQ=255)
                                                                                                                                aCTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGaa                                  >  1:168261/1‑133 (MQ=255)
                                                                                                                                aCTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAg                              <  3:117342/137‑1 (MQ=255)
                                                                                                                                aCTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAAcc                  >  7:169400/1‑149 (MQ=255)
                                                                                                                                aCTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAAcc                  >  8:257288/1‑149 (MQ=255)
                                                                                                                                   gATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTg               >  2:270530/1‑149 (MQ=255)
                                                                                                                                    aTCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCt                <  1:193163/147‑1 (MQ=255)
                                                                                                                                    aTCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCt                <  3:287651/147‑1 (MQ=255)
                                                                                                                                    aTCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGCTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTg               <  7:300333/148‑1 (MQ=255)
                                                                                                                                    aTCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGCTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGt              >  8:340467/1‑149 (MQ=255)
                                                                                                                                     tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGt              >  3:78343/1‑148 (MQ=255)
                                                                                                                                      tgacaaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAaga                                     >  6:119438/6‑124 (MQ=255)
                                                                                                                                      tgacaaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAaga                                     <  5:119438/119‑1 (MQ=255)
                                                                                                                                       cTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCa            <  1:453901/148‑1 (MQ=255)
                                                                                                                                        ttGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCAt           >  4:84105/1‑148 (MQ=255)
                                                                                                                                        ttGACCCCGCTGCGGCGGGCGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTAACAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTAACCCTGAAGAGAATCACGGATTCGCAACCCAGTCAt           >  2:177787/1‑148 (MQ=255)
                                                                                                                                        acaaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCAt           >  4:62308/4‑148 (MQ=255)
                                                                                                                                          aaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCt                          >  5:200306/2‑131 (MQ=255)
                                                                                                                                          aaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCt                          <  6:200306/130‑1 (MQ=255)
                                                                                                                                           aCCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAaga                                     <  2:347126/119‑1 (MQ=255)
                                                                                                                                           aCCCCGCTGCGGCGGGGTTGTCATTTGCGTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAaga                                     >  1:347126/1‑119 (MQ=255)
                                                                                                                                             cccGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAag       >  1:25419/1‑147 (MQ=255)
                                                                                                                                                gCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAggagga  <  6:289349/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGA  >  NC_000913/1335535‑1335827

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