Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,160,692 C→T 5.7% I769I (ATC→ATT mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,160,6920CT 5.7% 410.3 / 8.9 120I769I (ATC→ATTmdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (56/57);  new base T (3/4);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

ACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTAACCTGCCGACCGGAAAATCGCTCTCGGACGC  >  NC_000913/2160545‑2160835
                                                                                                                                                   |                                                                                                                                               
   aGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAAt                                                                                                                                              >  4:894145/1‑148 (MQ=255)
        ccAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         >  1:1060862/1‑148 (MQ=255)
        ccAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         <  3:894145/148‑1 (MQ=255)
        ccAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         <  2:1060862/148‑1 (MQ=255)
         cAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGaa                                                                                                                                        <  2:293991/148‑1 (MQ=255)
              gTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCaa                                                                                                                                   >  4:871139/1‑148 (MQ=255)
              gTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCaa                                                                                                                                   <  1:477792/148‑1 (MQ=255)
               tCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCaaa                                                                                                                                  <  3:871139/148‑1 (MQ=255)
                    tAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGAt                                                                                                                             >  3:272757/1‑148 (MQ=255)
                         aaCGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         >  1:392521/1‑131 (MQ=255)
                         aaCGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         <  2:392521/131‑1 (MQ=255)
                           cGCCTTCGGTCAGCGGCAAATCTCGTCCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGaa                                                                                                                                        >  4:56625/1‑130 (MQ=255)
                           cGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGaa                                                                                                                                        <  3:56625/130‑1 (MQ=255)
                                 cGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTAt                                                                                                                <  4:272757/148‑1 (MQ=255)
                                             aaTCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATgg                                                                                                    >  4:655238/1‑148 (MQ=255)
                                                ctcGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCaa                                                                                                 >  1:555463/1‑148 (MQ=255)
                                                ctcGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCaa                                                                                                 >  3:149208/1‑148 (MQ=255)
                                                    aCCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                                                                                           >  4:129128/1‑134 (MQ=255)
                                                    aCCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCg                                                                                              >  2:171297/1‑147 (MQ=255)
                                                    aCCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCg                                                                                              <  1:171297/147‑1 (MQ=255)
                                                      cATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAg                                                                                                                                 <  2:636575/110‑1 (MQ=255)
                                                      cATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCg                                                                                           >  2:623965/1‑148 (MQ=255)
                                                       aTTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                                                                                           >  3:398731/1‑131 (MQ=255)
                                                       aTTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGc                                                                                                           <  4:398731/131‑1 (MQ=255)
                                                        tttACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGaa                                                                                         <  3:43888/148‑1 (MQ=255)
                                                              aGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAg                                                                                                                                 >  3:1023498/1‑102 (MQ=255)
                                                              aGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAg                                                                                                                                 <  4:1023498/102‑1 (MQ=255)
                                                                      aaCCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCgg                                                                           <  1:623965/148‑1 (MQ=255)
                                                                        ccAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCgg                                                                                             >  2:951924/1‑128 (MQ=255)
                                                                         cAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCaaa                                                                                                                                  >  3:1143534/1‑90 (MQ=255)
                                                                             ataAAGTGGTGATGGAAGTGGATCCGCGCTATGCCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCgg                                                                           >  3:721466/1‑141 (MQ=255)
                                                                             ataAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCgg                                                                           <  3:714929/141‑1 (MQ=255)
                                                                             ataAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCgg                                                                           >  4:714929/1‑141 (MQ=255)
                                                                             ataAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAatca                                                                    <  1:986445/148‑1 (MQ=255)
                                                                             ataAAGCGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATAgg                                                                           <  4:721466/141‑1 (MQ=255)
                                                                              taAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAatcat                                                                   <  4:149208/148‑1 (MQ=255)
                                                                                aaGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTAGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAAtcatca                                                                 <  3:655238/148‑1 (MQ=255)
                                                                                 aGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAg                                                                <  4:591182/148‑1 (MQ=255)
                                                                                 aGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAg                                                                <  4:52864/148‑1 (MQ=255)
                                                                                    ggtgATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGc                                                                                                                        <  2:392200/89‑1 (MQ=255)
                                                                                         tGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCg                                                                            <  3:465091/128‑1 (MQ=255)
                                                                                         tGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCg                                                                            >  2:520724/1‑128 (MQ=255)
                                                                                           gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATTGATAACGaaca                                                                                                                                      >  1:322801/1‑66 (MQ=38)
                                                                                           gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATTGATAACGaaca                                                                                                                                      <  2:322801/68‑3 (MQ=38)
                                                                                           gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         <  3:374752/65‑1 (MQ=255)
                                                                                           gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCgg                                                                                             <  1:343773/109‑1 (MQ=255)
                                                                                           gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCgg                                                                                             >  2:343773/1‑109 (MQ=255)
                                                                                           gAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATcgg                                                      >  4:939915/1‑148 (MQ=255)
                                                                                            aaGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGa                                                                        >  1:24115/1‑129 (MQ=255)
                                                                                            aaGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGa                                                                        <  2:24115/129‑1 (MQ=255)
                                                                                             aGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATTGATAACGaa                                                                                                                                        <  4:391325/64‑1 (MQ=39)
                                                                                             aGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGaa                                                                                                                                        <  4:464930/64‑1 (MQ=255)
                                                                                             aGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGaa                                                                                                                                        >  3:464930/1‑64 (MQ=255)
                                                                                             aGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTa                                                                                                         >  4:741402/1‑94 (MQ=255)
                                                                                               tGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCaa                                                                                                 <  1:49978/101‑1 (MQ=255)
                                                                                               tGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCaa                                                                                                 >  2:49978/1‑101 (MQ=255)
                                                                                                 gATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAgg                                                               <  3:815293/133‑1 (MQ=255)
                                                                                                 gATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAgg                                                               >  4:815293/1‑133 (MQ=255)
                                                                                                  aTCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCg                                                                            <  3:572469/119‑1 (MQ=255)
                                                                                                  aTCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCg                                                                            >  4:572469/1‑119 (MQ=255)
                                                                                                   tCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAg                                                                >  3:621372/1‑130 (MQ=255)
                                                                                                   tCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAg                                                                <  4:621372/130‑1 (MQ=255)
                                                                                                    ccGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTa                                                          <  2:714610/135‑1 (MQ=255)
                                                                                                    ccGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTa                                                          >  1:714610/1‑135 (MQ=255)
                                                                                                    ccGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGa                                             >  2:1084219/1‑148 (MQ=255)
                                                                                                              cccAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCgg                                                                           <  3:158521/108‑1 (MQ=255)
                                                                                                              cccAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCgg                                                                           >  4:158521/1‑108 (MQ=255)
                                                                                                              cccAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGtt                                   <  1:1084219/148‑1 (MQ=255)
                                                                                                                cAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTa                                                          <  1:542197/123‑1 (MQ=255)
                                                                                                                cAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTa                                                          >  2:542197/1‑123 (MQ=255)
                                                                                                                 aagACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGt                                                                          >  4:760896/3‑106 (MQ=255)
                                                                                                                 aGGACATCAGTGCGCTGGAATAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTaa                                <  3:939915/148‑1 (MQ=255)
                                                                                                                 aGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGt                                                                          <  3:760896/106‑1 (MQ=255)
                                                                                                                 aGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTaa                                >  4:457580/1‑148 (MQ=255)
                                                                                                                  ggACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         <  4:577829/42‑1 (MQ=255)
                                                                                                                  ggACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGa                                                                                                                                         >  3:577829/1‑42 (MQ=255)
                                                                                                                    aCATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAatcat                                                                   >  3:952242/1‑110 (MQ=255)
                                                                                                                    aCATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTaa                                <  3:943680/145‑1 (MQ=255)
                                                                                                                    aCATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTaa                                >  4:943680/1‑145 (MQ=255)
                                                                                                                       tCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATg                                                                                       <  2:527825/87‑1 (MQ=255)
                                                                                                                       tCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATg                                                                                       >  1:527825/1‑87 (MQ=255)
                                                                                                                                   aaaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATcggcgg                                                   <  4:849016/111‑1 (MQ=255)
                                                                                                                                   aaaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATcggcgg                                                   >  1:719588/1‑111 (MQ=255)
                                                                                                                                    aaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAAtcatca                                                                 >  4:878065/1‑96 (MQ=255)
                                                                                                                                    aaaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAAtcatca                                                                 <  3:878065/96‑1 (MQ=255)
                                                                                                                                     aaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCg                                     >  2:647770/1‑123 (MQ=255)
                                                                                                                                     aaaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCg                                     <  1:647770/123‑1 (MQ=255)
                                                                                                                                      aaaTGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGa                                             <  3:886079/114‑1 (MQ=255)
                                                                                                                                           ttCGTTATTGATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCACCAGGGATTATCGGCGGCCTc                                               >  4:756021/1‑107 (MQ=255)
                                                                                                                                           ttCGTTATTGATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCACCAGGGATTATCGGCGGCCTc                                               <  3:756021/107‑1 (MQ=255)
                                                                                                                                           ttCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCg                                     >  2:232339/1‑117 (MQ=255)
                                                                                                                                            tCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTAt                                                                              <  1:1125184/75‑1 (MQ=255)
                                                                                                                                               ttATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGaa                                                                       <  1:316029/79‑1 (MQ=255)
                                                                                                                                               ttATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTAACCTGCCGACCGGAAAATCGCTCTCGGACGc  >  1:839198/1‑148 (MQ=255)
                                                                                                                                                tATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCg                                     <  3:583629/112‑1 (MQ=255)
                                                                                                                                                tATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCg                                     >  4:898083/1‑112 (MQ=255)
                                                                                                                                                tATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTAACCTGc                           >  1:444110/1‑122 (MQ=255)
                                                                                                                                                   cAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAatcat                                                                   >  1:958005/1‑79 (MQ=255)
                                                                                                                                                   cAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAatcat                                                                   <  2:958005/79‑1 (MQ=255)
                                                                                                                                                   cAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTa                                                          <  3:807719/88‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTTATCAATAACGAAGGCAAAGCGATCCCGCTGTCGTATTTCGCTAAATGGCAACCGGCGAATGCCCCACTATCGGTGAATCATCAGGGATTATCGGCGGCCTCGACCATTTCGTTTAACCTGCCGACCGGAAAATCGCTCTCGGACGC  >  NC_000913/2160545‑2160835

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