Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 853,967 G→A 5.5% P107P (CCG→CCA ybiR → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913853,9670GA 5.5% 230.6 / 2.4 72P107P (CCG→CCAybiRputative transporter
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base A (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

GAGTTAAGCGGTTATTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCG  >  NC_000913/853821‑854111
                                                                                                                                                  |                                                                                                                                                
gAGTTAAGCGGTTATTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCg                                                                                                                                                  <  1:90877/147‑1 (MQ=255)
    tAAGCGGTTATTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGAc                                                                                                                                             >  2:392396/1‑148 (MQ=255)
       gCGGTTATTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTAt                                                                                                                                          <  1:392396/148‑1 (MQ=255)
        cGGTTATTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTAt                                                                                                                                          >  5:106891/1‑147 (MQ=255)
          gTTATTTTTATGTTCTGGGGCGCAAAATTGTGCGCCGCTTTGCTACGGATCGTCTGCTGGCGATGTTTATTGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCAc                                                                                                                                       <  5:348472/148‑1 (MQ=255)
             aTTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGc                                                                                                                                     <  1:195874/147‑1 (MQ=255)
                  gATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTaaaa                                                                                                                                <  2:334015/147‑1 (MQ=255)
                  gATGTGCTGGGGCGCAAAATGGGGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTaaaa                                                                                                                                <  1:348095/147‑1 (MQ=255)
                        cTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACtgt                                                                                                                         >  3:12122/1‑148 (MQ=255)
                           gggCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACtgtg                                                                                                                        <  1:139237/146‑1 (MQ=255)
                           gggCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACtgtg                                                                                                                        <  2:106078/146‑1 (MQ=255)
                           gggCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGAGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACtgt                                                                                                                         <  4:12122/145‑1 (MQ=255)
                            ggCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTga                                                                                                                     <  6:106891/148‑1 (MQ=255)
                             gcgcAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCAGCGTTCGACCTTTCTGACCAACGATGTCGCGCTGTTTAATGGTGTTCCGCTGCCATTCACGCTAAAAAGACtgtgtgc                                                                                                                     >  6:96435/1‑146 (MQ=255)
                              cgcAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACtgtgt                                                                                                                       >  3:348959/1‑144 (MQ=255)
                                   aaTGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCGAAAAAGACTGTGTGAGATCCCg                                                                                                              >  5:133039/1‑148 (MQ=255)
                                    aTGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTGTAAGGTGCTGGCGGCGGCGCTGCGTTCTACCTTTCGGACCAAAGAAGTCGCGATTTTTATTGTTGTTCCGCTTACTAAGACGCTAAAAAGAATGTGTGAGATCCCgg                                                                                                             >  2:141747/1‑148 (MQ=255)
                                               cTTTGCTACGGAGCGTCGGCTGGCTATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTg                                                                                                  <  3:248502/148‑1 (MQ=255)
                                               cTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGGTGGCGGCGGCGGTGCTTTCTAACTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATACCGGTTAATCGGCt                                                                                                   >  5:420084/1‑147 (MQ=255)
                                                       cGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTAtttt                                                                                          <  7:435156/148‑1 (MQ=255)
                                                       cGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTAtttt                                                                                          <  5:145391/148‑1 (MQ=255)
                                                         gAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTg                                                                                        <  6:65688/148‑1 (MQ=255)
                                                         gAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTg                                                                                        <  4:348959/148‑1 (MQ=255)
                                                             gTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGatt                                                                                               >  2:288367/1‑137 (MQ=255)
                                                             gTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGatt                                                                                               <  1:288367/137‑1 (MQ=255)
                                                                ggCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCt                                                                                 <  4:93014/148‑1 (MQ=255)
                                                                 gCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTg                                                                                >  3:373633/1‑148 (MQ=255)
                                                                     gCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGc                                                                              <  4:373633/146‑1 (MQ=255)
                                                                           tttATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGgcgc                                                                                  >  1:75105/1‑136 (MQ=255)
                                                                           tttATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGgcgc                                                                                  <  2:75105/136‑1 (MQ=255)
                                                                            ttATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGc                                                                     >  6:273979/1‑148 (MQ=255)
                                                                               tGGGGCTGGCGGCGGCGATGCTTTAAACCTTTCTCACCAACGATGTCGCGCTGTTTATTGTTGATCCGCTGACTATAACGCTAAAAAGACTGTGTGAGAGCCCGGTTACTCTGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTgg                                                                  <  6:420084/148‑1 (MQ=255)
                                                                                 gTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGCGGCGCTGGCGGTCAACGGTGGtt                                                                >  7:71752/1‑148 (MQ=255)
                                                                                  tGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGatt                                                                                               <  1:286415/116‑1 (MQ=255)
                                                                                  tGCTGGCGGCCGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGatt                                                                                               >  2:286415/1‑116 (MQ=255)
                                                                                   gCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCg                                                              >  4:126505/1‑148 (MQ=255)
                                                                                   gCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCg                                                              >  1:226584/1‑148 (MQ=255)
                                                                                       gcggcggcgCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            <  6:133039/146‑1 (MQ=255)
                                                                                       gcggcggcgCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTa                                                           <  2:226584/147‑1 (MQ=255)
                                                                                             gcgcTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTg                                                        <  8:433307/144‑1 (MQ=255)
                                                                                             gcgcTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTg                                                        >  7:433307/1‑144 (MQ=255)
                                                                                                   ccttCTACCATTCAGCCCCACGTTGTCGCGGTGGTTATTGGTGTTCCGCTGTCTATCAAGCTAAAAAGACTTGGTGAGCTCCCGGGTAATCGGCTGCTTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTg                                              <  8:71752/146‑1 (MQ=255)
                                                                                                          ccTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAAccc                                       >  3:180195/1‑148 (MQ=255)
                                                                                                             ttCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACCCCAATTGGCACCCCGc                                     >  5:187874/1‑147 (MQ=255)
                                                                                                              tCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGAGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            >  2:219192/1‑123 (MQ=255)
                                                                                                              tCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGAGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            <  1:219192/123‑1 (MQ=255)
                                                                                                                 gACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCaaaa                                 >  3:165437/1‑147 (MQ=255)
                                                                                                                 gACCAACGATGTCGCGCTGTTTATTGTTGTTCCACTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCaaaa                                 >  2:80526/1‑147 (MQ=255)
                                                                                                                   ccAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAg                                                                                                                              <  8:274738/52‑1 (MQ=255)
                                                                                                                   ccAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAg                                                                                                                              >  7:274738/1‑52 (MQ=255)
                                                                                                                   ccAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAatat                              >  3:151270/1‑148 (MQ=255)
                                                                                                                    cAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATAtt                             <  4:180195/148‑1 (MQ=255)
                                                                                                                     aaCGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTc                            >  7:300748/1‑148 (MQ=255)
                                                                                                                     aaCGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGACAATTGGCAACCCGCAAAATATTc                            >  3:98576/1‑148 (MQ=255)
                                                                                                                      aCGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCaa                                                 <  7:344398/126‑1 (MQ=255)
                                                                                                                      aCGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCaa                                                 >  8:344398/1‑126 (MQ=255)
                                                                                                                      aCGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCt                           <  4:165437/148‑1 (MQ=255)
                                                                                                                       cGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCtt                          <  3:126505/148‑1 (MQ=255)
                                                                                                                           gTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            <  4:50201/110‑1 (MQ=255)
                                                                                                                           gTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            >  3:50201/1‑110 (MQ=255)
                                                                                                                             cgcgCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTgg                    >  5:179564/1‑148 (MQ=255)
                                                                                                                               cgcTGTTTATTGTTGTTCCGCTGCCTATCACCCTAAAAAGACTGTGTCAGATCCCGGTTAAGCGGATGAATATTTTAGACCCGCAGGCAGTCAAAGATGGTTCTATAATTACGCCAATTGGCAACCCGCAAAATATTCTTATCTggg                   >  3:297702/1‑147 (MQ=255)
                                                                                                                                 cTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGAc                >  6:110990/1‑148 (MQ=255)
                                                                                                                                  tGTTTATTGTTGTTCCACTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAAcc                                        <  6:90378/123‑1 (MQ=255)
                                                                                                                                  tGTTTATTGTTGTTCCACTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAAcc                                        >  5:90378/1‑123 (MQ=255)
                                                                                                                                    tttATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            <  3:192765/101‑1 (MQ=255)
                                                                                                                                    tttATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCt                                                            >  4:192765/1‑101 (MQ=255)
                                                                                                                                     ttATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTc            >  1:59822/1‑148 (MQ=255)
                                                                                                                                        ttgttgTTCCACTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTg          <  1:80526/147‑1 (MQ=255)
                                                                                                                                           ttgttCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCaa                                   <  7:127470/119‑1 (MQ=255)
                                                                                                                                           ttgttCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCaa                                   >  8:127470/1‑119 (MQ=255)
                                                                                                                                              ttCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCttt    <  2:59822/147‑1 (MQ=255)
                                                                                                                                                ccGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCg  <  4:98576/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GAGTTAAGCGGTTATTTTGATGTGCTGGGGCGCAAAATGGTGCGCCGCTTTGCTACGGAGCGTCGGCTGGCGATGTTTATGGTGCTGGCGGCGGCGCTGCTTTCTACCTTTCTGACCAACGATGTCGCGCTGTTTATTGTTGTTCCGCTGACTATCACGCTAAAAAGACTGTGTGAGATCCCGGTTAATCGGCTGATTATTTTTGAGGCGCTGGCAGTCAACGCTGGTTCGCTACTGACGCCAATTGGCAACCCGCAAAATATTCTTATCTGGGGACGTTCTGGTCTTTCG  >  NC_000913/853821‑854111

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