Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,672,847 A→G 5.2% A101A (GCT→GCC)  mdtI ← multidrug efflux system transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,8470AG 5.2% 238.5 / 2.6 78A101A (GCT→GCC) mdtImultidrug efflux system transporter
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCTTGA  >  NC_000913/1672701‑1672991
                                                                                                                                                  |                                                                                                                                                
aaaGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGc                                                                                                                                                <  4:47254/149‑1 (MQ=255)
               gcCGATAACCTGGCCATTTCGGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAATTAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAg                                                                                                                                 <  3:271227/149‑1 (MQ=255)
               gcCGATAACCGTGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCGTCATCAGGCAAGTATCACCATGAGCATACAAGCCAACAGCAAGAcact                                                                                                                                 >  4:59774/1‑146 (MQ=255)
               gcCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAg                                                                                                                                 >  3:201636/1‑149 (MQ=255)
                   aTAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGcc                                                                                                                              >  8:101631/1‑148 (MQ=255)
                       ccGGGCCATTTCTGTTACCGGTTTTCGCGACATGCTGGATCAGCTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTGCAGCAGGCAAGTTTCACCATGATCATTCCAGCCACCAGCAAGACCAGGCCAATcc                                                                                                                         >  1:178139/1‑149 (MQ=255)
                       ccGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGc                                                                                                                                         <  8:268906/133‑1 (MQ=255)
                       ccGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGc                                                                                                                                         >  7:268906/1‑133 (MQ=255)
                        cGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGCCAGCGCGGGCAGCGTCTTCATCAGGCAATTTTCACCATGATCATTCCAGCCACCTGCAAGACCAGGCCAATACa                                                                                                                        >  8:375689/1‑149 (MQ=255)
                        cGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCCGCGTCTTCATCAGGCAAGTTTCACCATGAGCATTCCAGCCAACCGCAAGACCAGCCCACTc                                                                                                                          >  5:346150/1‑147 (MQ=255)
                           gCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGcc                                                                                                                     <  5:193650/149‑1 (MQ=255)
                             cATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCtt                                                                                                                   <  1:153203/149‑1 (MQ=255)
                              aTTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCttt                                                                                                                  <  6:10030/149‑1 (MQ=255)
                               tttCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTa                                                                                                                 >  6:112515/1‑149 (MQ=255)
                               tttCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTa                                                                                                                 <  8:24530/149‑1 (MQ=255)
                                ttCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTAc                                                                                                                >  7:59066/1‑149 (MQ=255)
                                  cTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGa                                                                                                              >  1:390229/1‑149 (MQ=255)
                                     ttACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGAtt                                                                                                            >  7:360348/1‑148 (MQ=255)
                                     ttACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGAt                                                                                                             <  7:188206/147‑1 (MQ=255)
                                        ccGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTaa                                                                                                         <  1:59089/148‑1 (MQ=255)
                                        ccGGTTTTCGCGACAAGCGGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTaa                                                                                                         <  2:309015/148‑1 (MQ=255)
                                         cGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCGTCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAAcc                                                                                                       >  2:189591/1‑149 (MQ=255)
                                            ttttCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCg                                                                                                      <  8:241268/147‑1 (MQ=255)
                                                   gACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCaaa                                                                                             >  1:158895/1‑149 (MQ=255)
                                                   gACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCaaa                                                                                             >  1:263140/1‑149 (MQ=255)
                                                   gACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCaaa                                                                                             <  3:108944/149‑1 (MQ=255)
                                                         cTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAAc                                                                                            >  3:342039/1‑144 (MQ=255)
                                                         cTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGAt                                                                                       <  2:263140/149‑1 (MQ=255)
                                                            gATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATcc                                                                                     >  2:244845/1‑148 (MQ=255)
                                                            gATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATcc                                                                                     >  3:296545/1‑148 (MQ=255)
                                                              tCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAAc                                                                                  <  8:360348/149‑1 (MQ=255)
                                                                aaCTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATAAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCt                                                                                >  5:112056/1‑149 (MQ=255)
                                                                 aCTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTg                                                                               <  8:59066/149‑1 (MQ=255)
                                                                   tGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCGGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCg                                                                             >  8:169971/1‑149 (MQ=255)
                                                                   tGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCGGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCg                                                                             >  8:164700/1‑149 (MQ=255)
                                                                       gCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTAAAGCAGGCAATTTTCACCATGATAATACACGCGAAGATAAAGAACAGGCCAATCCAGCAATTACGATTTAAACGATGACAAAACAAGAACCAAACTGCGTCt                                                                          >  6:71268/1‑148 (MQ=255)
                                                                        cAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGCCCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTaa                                                                        >  3:314424/1‑149 (MQ=255)
                                                                               aatGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTggcg                                                                 >  4:180798/1‑149 (MQ=255)
                                                                                atGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTggcgg                                                                <  4:152059/149‑1 (MQ=255)
                                                                                   ttCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCaa                                                             <  2:390229/149‑1 (MQ=255)
                                                                                      aaaCGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATAc                                                          >  5:96870/1‑149 (MQ=255)
                                                                                       aaCGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATAc                                                          >  7:13998/1‑148 (MQ=255)
                                                                                          gCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCaaa                                                      <  6:112056/149‑1 (MQ=255)
                                                                                            tGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCACTCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATCCCAAccc                                                    >  6:194574/1‑149 (MQ=255)
                                                                                              cccGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCg                                                  >  3:338493/1‑149 (MQ=255)
                                                                                              cccGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCg                                                  <  1:244845/149‑1 (MQ=255)
                                                                                               ccGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGc                                                 >  7:353459/1‑149 (MQ=255)
                                                                                                cGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTa                                                                         <  4:161384/124‑1 (MQ=255)
                                                                                                cGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTa                                                                         >  3:161384/1‑124 (MQ=255)
                                                                                                cGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTATCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGcc                                                <  2:114091/149‑1 (MQ=255)
                                                                                                    aGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAAc                                                                                                        <  6:370989/89‑1 (MQ=255)
                                                                                                    aGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAAc                                                                                                        >  5:370989/1‑89 (MQ=255)
                                                                                                     gcgcGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCcaca                                           <  4:338493/149‑1 (MQ=255)
                                                                                                       gcgGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAAt                                         <  4:201636/149‑1 (MQ=255)
                                                                                                         gggCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGc                                       <  3:180798/149‑1 (MQ=255)
                                                                                                              gCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTTCCCAAACAAGATCCAACATGCGGCTAACGTGGAGGCAATACCAACCCCGCCCCACAATGc                                       >  8:381186/1‑144 (MQ=255)
                                                                                                              gCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGc                                       <  7:381186/144‑1 (MQ=255)
                                                                                                              gCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAg                                  <  6:96870/149‑1 (MQ=255)
                                                                                                              gCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAg                                  <  5:112515/149‑1 (MQ=255)
                                                                                                                   ttCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCga                               >  8:174542/1‑147 (MQ=255)
                                                                                                                   ttCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCga                               <  7:174542/147‑1 (MQ=255)
                                                                                                                     catcaGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGTCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCgac                              <  5:153471/146‑1 (MQ=255)
                                                                                                                      atcaGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCCAACAAGATCCAACATGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCgacagac                          <  3:59774/149‑1 (MQ=255)
                                                                                                                      atcaGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCgacagac                          <  4:296545/149‑1 (MQ=255)
                                                                                                                             aaGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCACCCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGa                   >  4:204169/1‑149 (MQ=255)
                                                                                                                                   cACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCt              <  6:277339/148‑1 (MQ=255)
                                                                                                                                    aCCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCtt             <  4:342039/148‑1 (MQ=255)
                                                                                                                                       aTGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCtt             <  6:301142/145‑1 (MQ=255)
                                                                                                                                       aTGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCtt             >  5:301142/1‑145 (MQ=255)
                                                                                                                                           tCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCt     >  8:267172/1‑149 (MQ=255)
                                                                                                                                           tCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACCATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGcccc                                              >  3:156413/1‑108 (MQ=255)
                                                                                                                                           tCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACCATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGcccc                                              <  4:156413/108‑1 (MQ=255)
                                                                                                                                            cATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAAGACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCtt    >  6:142913/1‑149 (MQ=255)
                                                                                                                                             aTTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCtt    >  6:127135/1‑148 (MQ=255)
                                                                                                                                              ttCCGGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCTTga  <  7:169971/149‑1 (MQ=255)
                                                                                                                                              ttCCGGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCTTga  <  7:164700/149‑1 (MQ=255)
                                                                                                                                              ttCCAGCCAACAGCAACACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAATATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCTTga  <  3:290227/149‑1 (MQ=255)
                                                                                                                                               tCCAGCCAACAGCAAGAACAGGCCAAACAAGCCTTAACGATTTAACCGTTGACCAAACCAGCTCCAACCTGCGGCGCACGTGGCGGACATACCACAccc                                                   >  8:273809/1‑99 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGACAAGCTGGATCAACTGTCGCAAAAATAATGTTCAAACGCTGCCCGACAGCGCGGGCAGCGTCTTCATCAGGCAAGTTTCACCATGATCATTCCAGCCAACAGCAAGACCAGGCCAATCCAGCCTTTACGATTTAACCGTTGACCAAACAAGATCCAACCTGCGGCTAACGTGGCGGCAATACCAAACCCGCCCCACAATGCATAAGCGACAGACAAGTCGATCCCTTTAACGGCTTGA  >  NC_000913/1672701‑1672991

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: ‑