Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,156,030 C→T 6.1% A256A (GCC→GCT mdtB → multidrug efflux system, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,156,0300CT 6.1% 207.8 / 2.4 65A256A (GCC→GCTmdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (2/2);  total (33/32)
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.44e-01

GACCAGCGAAACCGTGCGCACCGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTC  >  NC_000913/2155886‑2156173
                                                                                                                                                |                                                                                                                                               
gACCAGCGAAACCGTGCGCACCGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTAc                                                                                                                                              >  3:321812/1‑148 (MQ=255)
     gCGAAACCGTGCGCACCGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGaa                                                                                                                                         <  7:310945/148‑1 (MQ=255)
             gTGCGCACCGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGgcgcgc                                                                                                                                 <  6:29369/148‑1 (MQ=255)
                    ccGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGCTGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCg                                                                                                                          >  4:327920/1‑148 (MQ=255)
                     cGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGt                                                                                                                         <  4:321812/148‑1 (MQ=255)
                      gCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGgcg                                                                                                                                    >  6:227827/1‑136 (MQ=255)
                      gCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGgcg                                                                                                                                    <  5:227827/136‑1 (MQ=255)
                          ttACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTggg                                                                                                                    <  5:71178/148‑1 (MQ=255)
                           tACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGc                                                                                                                   >  3:370906/1‑148 (MQ=255)
                           tACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGc                                                                                                                   >  3:94578/1‑148 (MQ=255)
                             ccGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCg                                                                                                                  <  6:132598/147‑1 (MQ=255)
                                    tAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGCACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGc                                                                                                           >  1:234429/1‑147 (MQ=255)
                                         ttAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTg                                                                                                      >  5:180249/1‑147 (MQ=255)
                                         ttAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGt                                                                                                     <  8:279957/148‑1 (MQ=255)
                                          tAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTg                                                                                                      >  5:199938/1‑146 (MQ=255)
                                           aaCTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTa                                                                                                    <  3:234521/147‑1 (MQ=255)
                                             cTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTagag                                                                                                 >  7:318860/1‑148 (MQ=255)
                                                ggCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCa                                                                                               <  4:168030/147‑1 (MQ=255)
                                                 gCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAg                                                                                             <  4:370906/148‑1 (MQ=255)
                                                           gCCGCGACGGCCCTTCCCTTGCGGTCACGCTTTCAGCTAACGACCAGATGCCATCCGACGAAGAGTATCGCCAGCTAATCATCGCCTCCCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGaaa                                                                                    <  4:318915/147‑1 (MQ=255)
                                                                      ccTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCg                                                                        <  8:362600/148‑1 (MQ=255)
                                                                         tCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCg                                                                                                                  >  6:387209/1‑103 (MQ=255)
                                                                         tCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCg                                                                                                                  <  5:387209/103‑1 (MQ=255)
                                                                         tCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGaaa                                                                                    >  1:230887/1‑133 (MQ=255)
                                                                         tCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGaaa                                                                                    <  2:230887/133‑1 (MQ=255)
                                                                         tCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGgcg                                                                     >  2:197254/1‑148 (MQ=255)
                                                                            cGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTg                                                                                          >  7:340384/1‑124 (MQ=255)
                                                                            cGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTg                                                                                          <  8:340384/124‑1 (MQ=255)
                                                                             gTGCGGTCACGCTTGCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTg                                                                 <  1:296350/148‑1 (MQ=255)
                                                                              tGCGGTCACGCTTTCCGCGAACGCCCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTgg                                                                >  5:116772/1‑148 (MQ=255)
                                                                                   tCACGCCTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGaa                                                           >  6:28306/1‑148 (MQ=255)
                                                                                    cACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTCGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGaa                                                           >  7:223998/1‑147 (MQ=255)
                                                                                      cGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACaa                                                        <  8:318860/148‑1 (MQ=255)
                                                                                      cGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACa                                                         >  7:390610/1‑147 (MQ=255)
                                                                                        cTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAg                                                      <  4:94578/148‑1 (MQ=255)
                                                                                          ttCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACaaa                                                       <  5:283369/145‑1 (MQ=255)
                                                                                          ttCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACaaa                                                       >  6:283369/1‑145 (MQ=255)
                                                                                                   cGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGc                                                                           >  7:192411/1‑116 (MQ=255)
                                                                                                   cGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGc                                                                           <  8:192411/116‑1 (MQ=255)
                                                                                                     aCCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGt                                                                  >  2:250776/1‑123 (MQ=255)
                                                                                                     aCCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGt                                                                  <  1:250776/123‑1 (MQ=255)
                                                                                                     aCCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATtgt                                         >  2:278399/1‑148 (MQ=255)
                                                                                                      ccAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATtgtg                                        <  2:234429/148‑1 (MQ=255)
                                                                                                          aTGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGtgatga                                    <  1:197254/148‑1 (MQ=255)
                                                                                                            gCAATCCGCCGAAGAGTATCGCCAGATAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGCGCGATGACGCAACTGTAGAGCAACGTGCAGAAAACAGCTGGATCGGCCCGTGGGCGAACAAAGAACAGGCCATTGTCATGAAt                                  >  3:166757/1‑148 (MQ=255)
                                                                                                                tCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGtt                               >  1:74801/1‑147 (MQ=255)
                                                                                                                 ccgccgAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGtga                                       <  5:90546/138‑1 (MQ=255)
                                                                                                                 ccgccgAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGtga                                       >  6:90546/1‑138 (MQ=255)
                                                                                                                      gAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAgcgc                         >  1:89076/1‑147 (MQ=255)
                                                                                                                        agagTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAAGGTTCAGCGCCAg                      >  5:317643/1‑148 (MQ=255)
                                                                                                                               cGCCAGCTTATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTg               >  1:203094/1‑148 (MQ=255)
                                                                                                                               cGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTgg                                                                >  5:92505/1‑99 (MQ=255)
                                                                                                                               cGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTgg                                                                <  6:92505/99‑1 (MQ=255)
                                                                                                                                   aGCTAATCATCGCTTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAAc                                                   >  8:322645/1‑108 (MQ=255)
                                                                                                                                   aGCTAATCATCGCTTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAAc                                                   <  5:14229/108‑1 (MQ=255)
                                                                                                                                   aGCTAATCATCGCTTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAAc                                                   <  7:322645/108‑1 (MQ=255)
                                                                                                                                   aGCTAATCATCGCTTACCAGAACGGCGCGCCAATTCGTCTGGGCGAAGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAAc                                                   >  6:14229/1‑108 (MQ=255)
                                                                                                                                    gCTACTGATCGCCTACCAGAACGGCGCGCAAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAAc          <  2:352855/148‑1 (MQ=255)
                                                                                                                                    gCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAAc          <  5:28306/148‑1 (MQ=255)
                                                                                                                                    gCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAAc          <  6:180249/148‑1 (MQ=255)
                                                                                                                                    gCTAATCATCGCCTACCAGAACGGCGAGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAAc          <  6:116772/148‑1 (MQ=255)
                                                                                                                                        atcatcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTGCAGCGCCAGCCCGGTGCTAACATTa      >  2:379204/1‑148 (MQ=255)
                                                                                                                                          catcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCCCCAGCCCGGTGCTAACAtt       <  1:150414/145‑1 (MQ=255)
                                                                                                                                          catcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCACCAGCCCGGTGCTAACAtt       >  2:150414/1‑145 (MQ=255)
                                                                                                                                            tcGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTAtctc  >  1:153561/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GACCAGCGAAACCGTGCGCACCGCCATTACCGGCGCTAACGTTAACTCGGCAAAAGGTAGCCTCGACGGCCCTTCCCGTGCGGTCACGCTTTCCGCGAACGACCAGATGCAATCCGCCGAAGAGTATCGCCAGCTAATCATCGCCTACCAGAACGGCGCGCCAATTCGTCTGGGCGATGTCGCAACTGTAGAGCAAGGTGCAGAAAACAGCTGGCTCGGCGCGTGGGCGAACAAAGAACAGGCCATTGTGATGAATGTTCAGCGCCAGCCCGGTGCTAACATTATCTC  >  NC_000913/2155886‑2156173

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