Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 513,465 T→G 5.5% L297R (CTG→CGG) ‡ ybaT → putative amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913513,4650TG 5.5% 472.8 / 16.2 128L297R (CTG→CGG) ‡ybaTputative amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (60/61);  new base G (3/4);  total (63/65)
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.17e-01

Only 100 of 132 total aligned reads displayed.

TTAGCGATGTATCGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGG  >  NC_000913/513321‑513612
                                                                                                                                                |                                                                                                                                                   
    cGATGTATCGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACAGGTTCGcg                                                                                                                                              <  2:329170/148‑2 (MQ=255)
           tCGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTa                                                                                                                                       >  1:454256/1‑148 (MQ=255)
            cGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTaa                                                                                                                                      >  3:816554/1‑148 (MQ=255)
            cGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTaa                                                                                                                                      >  1:852475/1‑148 (MQ=255)
            cGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTaa                                                                                                                                      >  1:842669/1‑148 (MQ=255)
            cGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACAGGTTCGcgt                                                                                                                                             >  1:1014560/1‑139 (MQ=255)
            cGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACAGGTTCGcgt                                                                                                                                             <  2:1014560/141‑3 (MQ=255)
              gCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAca                                                                                                                                    <  2:842669/148‑1 (MQ=255)
               cATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAcat                                                                                                                                   >  1:476321/1‑148 (MQ=255)
               cATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAcat                                                                                                                                   <  4:816554/148‑1 (MQ=255)
                 ttAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGc                                                                                                                                               >  2:748610/1‑134 (MQ=255)
                    gagTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATgg                                                                                                                              <  2:454256/148‑1 (MQ=255)
                     agTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTa                                                                                                                                       >  1:1014415/1‑138 (MQ=255)
                     agTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGa                                                                                                                             <  2:476321/148‑1 (MQ=255)
                        tAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAc                                                                                                                                     >  3:1051077/1‑137 (MQ=255)
                        tAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGacaa                                                                                                                          >  3:14453/1‑148 (MQ=255)
                              aaTATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCgtg                                                                                                                                           >  2:160392/1‑125 (MQ=255)
                              aaTATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCgtg                                                                                                                                           <  1:160392/125‑1 (MQ=255)
                                     cGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAcatc                                                                                                                                  >  2:666693/1‑127 (MQ=255)
                                      gATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTaa                                                                                                                                      <  3:1085889/122‑1 (MQ=255)
                                      gATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTaa                                                                                                                                      >  4:1085889/1‑122 (MQ=255)
                                            gccgTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaacgcgaac                                                                                                      >  3:828982/1‑148 (MQ=255)
                                             ccgTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACt                                                                                                     <  4:828982/148‑1 (MQ=255)
                                              cgTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCgtg                                                                                                                                           >  3:226452/1‑109 (MQ=255)
                                                     cAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGc                                                                                             <  3:179010/148‑1 (MQ=255)
                                                            ctTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGacaac                                                                                                                         <  4:750768/113‑1 (MQ=255)
                                                            ctTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGaa                                                                                      <  2:852475/148‑1 (MQ=255)
                                                                 ccGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAAt                                                                                 >  3:453777/1‑148 (MQ=255)
                                                                      gctCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCtg                                                                            <  4:453777/148‑1 (MQ=255)
                                                                        tCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGa                                                                                                                             <  3:709059/97‑1 (MQ=255)
                                                                         cGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTaa                                                                                          <  4:622198/131‑1 (MQ=255)
                                                                               tgtgGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATaaa                                                                                  >  1:380833/1‑133 (MQ=255)
                                                                               tgtgGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATaaa                                                                                  <  2:380833/133‑1 (MQ=255)
                                                                                gtgGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCaga                                                                  <  1:852129/148‑1 (MQ=255)
                                                                                   ggTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGttt                                                                                                                                        >  2:985291/1‑75 (MQ=255)
                                                                                   ggTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGttt                                                                                                                                        <  1:985291/75‑1 (MQ=255)
                                                                                   ggTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                               <  4:399236/148‑1 (MQ=255)
                                                                                          tgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAcatcat                                                                                                                                >  1:746943/1‑76 (MQ=255)
                                                                                          tgATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAAcatcat                                                                                                                                <  2:746943/76‑1 (MQ=255)
                                                                                            aTCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATTAATCCCTGTGGCGGCag                                                                   >  1:731871/1‑135 (MQ=255)
                                                                                            aTCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCag                                                                   <  2:731871/135‑1 (MQ=255)
                                                                                            aTCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCa                                                      <  2:853814/148‑1 (MQ=255)
                                                                                             tcgtcgTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTCCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTa                                                                                           <  1:1090661/110‑1 (MQ=255)
                                                                                             tcgtcgTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTCCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTa                                                                                           >  2:1090661/1‑110 (MQ=255)
                                                                                              cgtcgtCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATgg                                                                                                                     >  1:611329/1‑83 (MQ=255)
                                                                                                 cgtcATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCtg                                                                            >  4:760493/1‑121 (MQ=255)
                                                                                                 cgtcATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCtg                                                                            <  3:760493/121‑1 (MQ=255)
                                                                                                       cGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACAt                                                                                                                       <  2:500286/72‑1 (MQ=255)
                                                                                                       cGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaacgcga                                                                                                        >  3:664286/1‑87 (MQ=255)
                                                                                                       cGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAgcgaacgcga                                                                                                        <  4:664286/87‑1 (MQ=255)
                                                                                                         gcgcTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATg                                                                                                                      <  4:179150/71‑1 (MQ=255)
                                                                                                         gcgcTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATg                                                                                                                      >  3:179150/1‑71 (MQ=255)
                                                                                                             tttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTa                                                                                           >  2:242283/1‑94 (MQ=255)
                                                                                                             tttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTa                                                                                           <  1:242283/94‑1 (MQ=255)
                                                                                                              ttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCagag                                                                 <  4:393744/119‑1 (MQ=255)
                                                                                                              ttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCagag                                                                 >  3:393744/1‑119 (MQ=255)
                                                                                                              ttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                               <  1:35558/121‑1 (MQ=255)
                                                                                                              ttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTa                                                               >  2:35558/1‑121 (MQ=255)
                                                                                                              ttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTggg                                                         >  1:477737/1‑127 (MQ=255)
                                                                                                              ttACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGgtgt                                    >  2:767998/1‑148 (MQ=255)
                                                                                                                    gCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTa                              <  1:767998/148‑1 (MQ=255)
                                                                                                                       aCGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTg                                  >  2:447159/1‑141 (MQ=255)
                                                                                                                          gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATgg                                                                                                                              >  1:551495/1‑46 (MQ=255)
                                                                                                                          gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATgg                                                                                                                              <  2:551495/46‑1 (MQ=255)
                                                                                                                          gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAAt                                                                                 <  4:601504/91‑1 (MQ=255)
                                                                                                                          gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAAt                                                                                 >  3:601504/1‑91 (MQ=255)
                                                                                                                          gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtggtg                                     >  3:1039327/1‑135 (MQ=255)
                                                                                                                          gCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtggtg                                     <  4:1039327/135‑1 (MQ=255)
                                                                                                                             tCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCgtgg                                       >  3:549017/1‑130 (MQ=255)
                                                                                                                                gCCATTAACGCGAACCTGTTCGCCGTGTTTAAcat                                                                                                                                   <  4:723791/35‑1 (MQ=255)
                                                                                                                                gCCATTAACGCGAACCTGTTCGCCGTGTTTAAcat                                                                                                                                   >  3:723791/1‑35 (MQ=255)
                                                                                                                                gCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGAcgg                  <  4:14453/148‑1 (MQ=255)
                                                                                                                                 ccATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTggcggc                                                                     <  1:1064533/96‑1 (MQ=255)
                                                                                                                                 ccATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTggcggc                                                                     >  2:1064533/1‑96 (MQ=255)
                                                                                                                                 ccATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCag                                                                   >  4:196624/1‑98 (MQ=255)
                                                                                                                                  cATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACat                                                  >  4:1025356/1‑114 (MQ=255)
                                                                                                                                  cATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGAt                                <  1:382228/132‑1 (MQ=255)
                                                                                                                                  cATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGAcggcg                >  2:560197/1‑148 (MQ=255)
                                                                                                                                    ttAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTg                                                                                                    >  3:1003835/1‑62 (MQ=255)
                                                                                                                                    ttAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTg                                                                                                    <  4:1003835/62‑1 (MQ=255)
                                                                                                                                     tAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTaa                                                                                          <  3:784482/71‑1 (MQ=255)
                                                                                                                                       cggcgAACAGGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATaa                                                                                   <  1:706317/74‑1 (MQ=255)
                                                                                                                                       cggcgAACAGGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATaa                                                                                   >  2:706317/3‑76 (MQ=255)
                                                                                                                                       aCGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATg                                                                                        >  1:755284/1‑71 (MQ=255)
                                                                                                                                       aCGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGaa                                                                                      >  1:916814/1‑73 (MQ=255)
                                                                                                                                       aCGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTggcg                                                                       <  1:979829/88‑1 (MQ=255)
                                                                                                                                       aCGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTggcg                                                                       >  2:979829/1‑88 (MQ=255)
                                                                                                                                        cgcgAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCt                                                            <  2:1070230/98‑1 (MQ=255)
                                                                                                                                        cgcgAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCt                                                            >  1:1070230/1‑98 (MQ=255)
                                                                                                                                        cgcgAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGa                                 <  4:2161/125‑1 (MQ=255)
                                                                                                                                        cgcgAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGa                                 >  3:2161/1‑125 (MQ=255)
                                                                                                                                        cgcgAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTAt                             >  2:824653/1‑129 (MQ=255)
                                                                                                                                          cgAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACggcggc              >  4:1104558/1‑142 (MQ=255)
                                                                                                                                            aaCCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTg                                  <  1:816250/120‑1 (MQ=255)
                                                                                                                                            aaCCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTg                                  >  2:816250/1‑120 (MQ=255)
                                                                                                                                              ccTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTa    <  1:560197/148‑1 (MQ=255)
                                                                                                                                                tGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCtga                        <  1:859854/126‑1 (MQ=255)
                                                                                                                                                tGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCtga                        >  2:859854/1‑126 (MQ=255)
                                                                                                                                                tGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAgg  >  2:624491/1‑148 (MQ=255)
                                                                                                                                                tGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAgg  >  4:843264/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TTAGCGATGTATCGGCATTAGAGTTAGAAAAATATGCCGATACCGCCGTAGCGCAGGCTGCTTCTCCGCTGCTCGGGCATGTGGGTTATGTGATCGTCGTCATCGGCGCTTTACTGGCGACGGCTTCAGCCATTAACGCGAACCTGTTCGCCGTGTTTAACATCATGGACAACATGGGCAGCGAACGCGAACTGCCGAAGCTAATGAATAAATCCCTGTGGCGGCAGAGTACCTGGGGCAACATTATTGTCGTGGTGTTGATTATGCTGATGACGGCGGCACTGAATTTAGG  >  NC_000913/513321‑513612

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