Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 397,560 A→T 5.5% N308Y (AAC→TAC)  sbmA → peptide antibiotic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913397,5600AT 5.5% 280.0 / 2.2 94N308Y (AAC→TAC) sbmApeptide antibiotic transporter
Reads supporting (aligned to +/- strand):  ref base A (42/46);  new base T (4/2);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GGCAGTGGTAGGGATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTACCAA  >  NC_000913/397412‑397708
                                                                                                                                                    |                                                                                                                                                    
ggCAGTGGTAGGGATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGAAGAGGCTGCCAACCGTAAAGAGCTGGTTTAAGGTGAAGACGATGCCCCGCGCGCGACGCCGCCTACGGTACGCGAGGTGTTTAGCGCCGTACGGaaaaa                                                                                                                                                     >  2:104178/1‑150 (MQ=255)
 gCAGTGGTAGGGATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAAc                                                                                                                                                    <  1:1590775/150‑1 (MQ=255)
   aGTGGTAGGGATCAAACTGCCGGGGCTGGAGTTTAAAACCCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTa                                                                                                                                                  >  2:832361/1‑150 (MQ=255)
       gTAGGGATCAAACTGCCGGGGCTGGAGATTAAAAACCAGCGTGTAGATGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTAAGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTc                                                                                                                                              >  2:921000/1‑150 (MQ=255)
         aGGGATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCg                                                                                                                                            >  1:951915/1‑150 (MQ=255)
         aGGGATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCg                                                                                                                                            <  1:1696839/150‑1 (MQ=255)
         aGGGATCAAACTGCAGGGGCTGGAGTTTAAAAACCAGCTTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCg                                                                                                                                            <  1:252290/150‑1 (MQ=255)
            gATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTCCGGTACGCGAGCTGTTTAGCGCCGTACGGAAAACCTATTTCCGCct                                                                                                                                         >  1:74365/1‑150 (MQ=255)
            gATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCct                                                                                                                                         >  2:723686/1‑150 (MQ=255)
             aTCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCctc                                                                                                                                        <  2:54884/150‑1 (MQ=255)
               cAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTATATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTa                                                                                                                                      >  1:1009733/1‑150 (MQ=255)
                  aCTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTAttt                                                                                                                                   <  2:3491/150‑1 (MQ=255)
                  aCTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTAttt                                                                                                                                   >  1:173630/1‑150 (MQ=255)
                     gCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCa                                                                                                                                >  1:571156/1‑150 (MQ=255)
                       cGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACt                                                                                                                              >  2:244984/1‑150 (MQ=255)
                          ggCTGGAGTTTAAAAACCAGCGGGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACtata                                                                                                                           <  2:960276/150‑1 (MQ=255)
                               gAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTAt                                                                                                                      <  1:7120/150‑1 (MQ=255)
                                aGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTAtt                                                                                                                     >  1:1704616/1‑150 (MQ=255)
                                 gTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTAttt                                                                                                                    <  1:503883/150‑1 (MQ=255)
                                     aaaaaCCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAAc                                                                                                                <  2:859822/150‑1 (MQ=255)
                                     aaaaaCCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAAc                                                                                                                >  2:176531/1‑150 (MQ=255)
                                      aaaaCCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTAt                                                                                                                                     <  1:1362300/128‑1 (MQ=255)
                                      aaaaCCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTAt                                                                                                                                     >  2:1362300/1‑128 (MQ=255)
                                       aaaCCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACAt                                                                                                              >  1:987907/1‑150 (MQ=255)
                                          ccAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGc                                                                                                           >  2:845188/1‑150 (MQ=255)
                                           cAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACta                                                                                                                             >  2:750026/1‑131 (MQ=255)
                                           cAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACta                                                                                                                             <  1:750026/131‑1 (MQ=255)
                                           cAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGcc                                                                                                          <  1:1172970/150‑1 (MQ=255)
                                             gCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCGCTATTTTCACTATATGTATTTCAACATCGCCCg                                                                                                        <  2:457190/150‑1 (MQ=255)
                                                 gtAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCAtt                                                                                                    <  1:1151653/150‑1 (MQ=255)
                                                 gtAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCAtt                                                                                                    <  2:951915/150‑1 (MQ=255)
                                                    gagGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATtctc                                                                                                 >  2:1277319/1‑150 (MQ=255)
                                                     tggCTGCCTACCGTAAAGAGCTTGTTTATGGTGAAGACGATGCCACGCGCGCGACTCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTctct                                                                                                <  2:1009733/149‑1 (MQ=255)
                                                      ggcggcCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTa                                                                                               <  2:1063878/150‑1 (MQ=255)
                                                      ggCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTa                                                                                               >  2:242409/1‑150 (MQ=255)
                                                      ggCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTa                                                                                               >  2:1031348/1‑150 (MQ=255)
                                                          gCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTg                                                                                           >  2:172133/1‑150 (MQ=255)
                                                            cTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCa                                                                                         >  2:1361071/1‑150 (MQ=255)
                                                              aCCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAgg                                                                                       <  2:1116330/150‑1 (MQ=255)
                                                               ccGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGt                                                                                      <  2:173630/150‑1 (MQ=255)
                                                                 gTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCg                                                                                    >  1:1022033/1‑150 (MQ=255)
                                                                    aaGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATa                                                                                 <  1:1276684/150‑1 (MQ=255)
                                                                        gCTGGTTTATGGTGAGGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTATATTTGCAGGTCGATAACGt                                                                             >  1:608788/1‑150 (MQ=255)
                                                                        gCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGt                                                                             >  2:95280/1‑150 (MQ=255)
                                                                        gCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGt                                                                             >  2:415497/1‑150 (MQ=255)
                                                                        gCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGt                                                                             >  1:380722/1‑150 (MQ=255)
                                                                         cTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGtt                                                                            <  2:226178/150‑1 (MQ=255)
                                                                           ggTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGtttt                                                                          >  2:461576/1‑150 (MQ=255)
                                                                             tttATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCg                                                                        <  2:763776/150‑1 (MQ=255)
                                                                             tttATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCg                                                                        <  2:739781/150‑1 (MQ=255)
                                                                               tATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGt                                                                      <  1:677603/150‑1 (MQ=255)
                                                                                aTGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGtt                                                                     <  1:176531/150‑1 (MQ=255)
                                                                                 tGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGttt                                                                    >  1:164178/1‑150 (MQ=255)
                                                                                   gTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAACCTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCTGGTCGATAACGTTTTCGGTTTGt                                                                  >  1:235017/1‑150 (MQ=255)
                                                                                     gAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTc                                                                >  2:1015990/1‑150 (MQ=255)
                                                                                        gACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTg                                                             >  2:1265884/1‑150 (MQ=255)
                                                                                           gATGCCACGCGCGCGATGCAGCCTACTGTACGCGAGCTGTTTAGCGCCGTACGGAAATACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGGATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTg                                                          <  1:654705/150‑1 (MQ=255)
                                                                                           gATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTg                                                          <  1:415497/150‑1 (MQ=255)
                                                                                             tGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGtt                                                        <  1:1387387/150‑1 (MQ=255)
                                                                                             tGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTAGCTGtt                                                        >  1:154234/1‑150 (MQ=255)
                                                                                                  cgcgcgcgACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCTTTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGt                                                   <  1:845188/150‑1 (MQ=255)
                                                                                                   gcgcgcgACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTc                                                  <  1:778605/150‑1 (MQ=255)
                                                                                                     gcgcgACGCCGCCTACGGTACGCGAGCTTTTTAGCGCCGTACGGAAAACTATTTCCGCCTCTAATTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAt                                               <  1:1314804/150‑1 (MQ=255)
                                                                                                     gcgcgACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCaa                                                >  2:884773/1‑150 (MQ=255)
                                                                                                     gcgcgACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCaa                                                >  2:691337/1‑150 (MQ=255)
                                                                                                      cgcgACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAt                                               <  1:723686/150‑1 (MQ=255)
                                                                                                       gcgACGCCGCCTACGGTACGCTAGCTGTTTAGCGCCGTACGTAAAAACTATTTCCTCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTTTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAtt                                              <  1:244984/150‑1 (MQ=255)
                                                                                                       gcgACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAtt                                              >  2:1668237/1‑150 (MQ=255)
                                                                                                           cgcctcCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAATACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACAACGCCCGCATTCTCTATTTGCAGGTCGATAACTTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAttgttg                                          <  2:1400703/150‑1 (MQ=255)
                                                                                                           cgccgcCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAttgttg                                          <  2:1441287/150‑1 (MQ=255)
                                                                                                           cgccgcCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAttgttg                                          <  1:172133/150‑1 (MQ=255)
                                                                                                           cgccgcCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCATGCTGTTTCCGTCAAttcttg                                          >  1:1453621/1‑150 (MQ=255)
                                                                                                              cgccTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACCATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCg                                       <  2:1024969/150‑1 (MQ=255)
                                                                                                              cgccTACGGTACGCGAGCTGTATAGCGCCGTACGGAAAAACTATATCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCATGATGTTTCCGTCAATTGTTGCCg                                       >  1:1623685/1‑150 (MQ=255)
                                                                                                               gccTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGGTGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCgg                                      <  1:877200/150‑1 (MQ=255)
                                                                                                                ccTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGt                                     >  2:1380567/1‑150 (MQ=255)
                                                                                                                ccTACGGTACGCGAGCTGTATAGCGCCGTACGGAAAAACTATTTCCGCCGCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACTTTTTCGGTTTGTTCTTGCTGTTTCCGTCATTAGTAGCCGGt                                     >  2:889042/1‑150 (MQ=255)
                                                                                                                  tACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTAc                                   >  2:159944/1‑150 (MQ=255)
                                                                                                                  tACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTAc                                   >  1:89460/1‑150 (MQ=255)
                                                                                                                       tACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGGCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTa                              <  1:1017943/150‑1 (MQ=255)
                                                                                                                         cgcgAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACg                            <  1:1361071/150‑1 (MQ=255)
                                                                                                                           cgAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCt                          >  1:1216100/1‑150 (MQ=255)
                                                                                                                               cTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGc                      <  2:164178/150‑1 (MQ=255)
                                                                                                                               cTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGc                      <  2:1383539/150‑1 (MQ=255)
                                                                                                                                tGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGcc                     <  1:921000/150‑1 (MQ=255)
                                                                                                                                     aGCGCCGTACGGAAAATCTATTTACGCCTCTTTTTTCTCTATTTGTATTTCAACTTCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAAATGTTGCCGGTACGATTACGCTCGGCCtgttg                <  2:779714/150‑1 (MQ=255)
                                                                                                                                          gaggcGGAAATAGTTTTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCa           >  2:1025588/5‑150 (MQ=255)
                                                                                                                                          gaggcGGAAATAGTTTTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCa           >  2:789159/5‑150 (MQ=255)
                                                                                                                                           gTACGGAAAAACTATATCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCAGTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTGCCGGCAATTGGTGCCGGTACGATTACGCTCGGCCTGATGACGCAg          <  2:955350/150‑1 (MQ=255)
                                                                                                                                            ggcGGAAATAGTTTTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGa         >  1:947864/3‑150 (MQ=255)
                                                                                                                                            ggcGGAAATAGTTTTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGa         >  1:54235/3‑150 (MQ=255)
                                                                                                                                             aCGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGAt        <  1:1179562/150‑1 (MQ=255)
                                                                                                                                               ggAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTa      <  2:1397779/150‑1 (MQ=255)
                                                                                                                                               ggAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTa      >  2:230883/1‑150 (MQ=255)
                                                                                                                                                 taaaaCTATTTACGCCTCTATTTTCACTATTTGTAATTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTAcc    <  2:235017/149‑1 (MQ=255)
                                                                                                                                                 aaaaaCTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTAcc    >  1:1186690/1‑150 (MQ=255)
                                                                                                                                                  aaaaCTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTACCa   >  2:282518/1‑150 (MQ=255)
                                                                                                                                                  aaaaCTATTTCCGCCTCTATTTTCACTATATGTATATCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTACCa   <  1:435792/150‑1 (MQ=255)
                                                                                                                                                   aaaCTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGTTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTACCaa  <  1:242409/150‑1 (MQ=255)
                                                                                                                                                   aaaCTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTACCaa  >  2:118146/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCAGTGGTAGGGATCAAACTGCCGGGGCTGGAGTTTAAAAACCAGCGTGTAGAGGCTGCCTACCGTAAAGAGCTGGTTTATGGTGAAGACGATGCCACGCGCGCGACGCCGCCTACGGTACGCGAGCTGTTTAGCGCCGTACGGAAAAACTATTTCCGCCTCTATTTTCACTATATGTATTTCAACATCGCCCGCATTCTCTATTTGCAGGTCGATAACGTTTTCGGTTTGTTCTTGCTGTTTCCGTCAATTGTTGCCGGTACGATTACGCTCGGCCTGATGACGCAGATTACCAA  >  NC_000913/397412‑397708

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: