Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 548,239 C→A 5.8% S520R (AGC→AGA fdrA → putative NAD(P)‑binding acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913548,2390CA 5.8% 243.9 / 3.6 77S520R (AGC→AGAfdrAputative NAD(P)‑binding acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base A (3/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAA  >  NC_000913/548091‑548382
                                                                                                                                                    |                                                                                                                                               
tGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGc                                                                                                                                                 >  1:590830/1‑149 (MQ=255)
 ggAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCt                                                                                                                                                <  1:123269/149‑1 (MQ=255)
  gAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTATCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCtt                                                                                                                                               >  1:925538/1‑149 (MQ=255)
    aGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTg                                                                                                                                             >  1:655469/1‑149 (MQ=255)
         cGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTg                                                                                                                                        <  2:1260227/149‑1 (MQ=255)
          ggggATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTgg                                                                                                                                       <  1:471994/149‑1 (MQ=255)
           gggATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGa                                                                                                                                      <  2:1408520/149‑1 (MQ=255)
           gggATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGa                                                                                                                                      <  2:376885/149‑1 (MQ=255)
             gATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGc                                                                                                                                    >  1:240892/1‑149 (MQ=255)
             gATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGc                                                                                                                                    <  1:1465452/149‑1 (MQ=255)
             gATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGc                                                                                                                                    >  1:838631/1‑149 (MQ=255)
                    gTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCAGCACCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCCTCATTACTCGAAAACGTCGCCGTGAATAACATCGGATTACGCAGCTTTGCGATGGAGCTACaacc                                                                                                                             >  2:273075/1‑147 (MQ=255)
                     tCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAACgc                                                                                                                            >  1:13178/1‑149 (MQ=255)
                        tGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGcc                                                                                                                         <  1:37821/149‑1 (MQ=255)
                         gAGTTCGCTACCGGAAGCCACCTTGCGGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCa                                                                                                                        <  1:902416/149‑1 (MQ=255)
                            ttCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGATTTGCGCTGGAGCTACAAAGCGCCAGCa                                                                                                                     >  2:1176645/1‑149 (MQ=255)
                            ttCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGATTTGCGCTGGAGCTACAAAGCGCCAGCa                                                                                                                     >  2:387584/1‑149 (MQ=255)
                            ttCGCTACAGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATACGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCa                                                                                                                     >  1:945412/1‑149 (MQ=255)
                               gCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAAc                                                                                                                  <  2:954446/149‑1 (MQ=255)
                                    cGGAAGCCACCTTGCTGGAGGCAGGGTTAATTCATCCTCTCTCGCATGCCGCACTGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGAGTAACATCGGATTACGCTGCGTTGCGCTGGAGCTACAAAGCGCCAGCAAACCggat                                                                                                             >  1:99452/1‑147 (MQ=255)
                                         gCCACCTTGCTGGCGGCAGCGTTTATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCa                                                                                                        >  1:1378736/1‑149 (MQ=255)
                                         gCCACCTTGCTGGCGGCAGCGTTAATTCATCTGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCa                                                                                                        <  1:492235/149‑1 (MQ=255)
                                           cACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCCCAGCAACACACCCCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCAtt                                                                                                      >  1:450406/1‑149 (MQ=255)
                                            aCCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTa                                                                                                     >  2:1195753/1‑149 (MQ=255)
                                             ccTTGCTTGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGGCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTAc                                                                                                    <  2:1455348/149‑1 (MQ=255)
                                              cTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTAcc                                                                                                   >  2:51128/1‑149 (MQ=255)
                                              cTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTAcc                                                                                                   >  1:176814/1‑149 (MQ=255)
                                              cTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGACGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTAcc                                                                                                   >  2:46235/1‑149 (MQ=255)
                                               ttGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCa                                                                                                  >  1:709881/1‑149 (MQ=255)
                                                   tGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACCGCAACCCACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCGGGCGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATgg                                                                                              >  1:1152566/1‑149 (MQ=255)
                                                     gcggcAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTc                                                                                            <  2:370685/149‑1 (MQ=255)
                                                     gcggcAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTc                                                                                            >  1:1119451/1‑149 (MQ=255)
                                                       ggcAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGc                                                                                          <  2:425507/149‑1 (MQ=255)
                                                             gTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCg                                                                                    >  1:64048/1‑149 (MQ=255)
                                                               tAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCCGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGcc                                                                                  >  1:27728/1‑149 (MQ=255)
                                                                    cATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCTTTACTTGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTgg                                                                             <  1:117488/149‑1 (MQ=255)
                                                                    cATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCCGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTgg                                                                             >  2:614147/1‑149 (MQ=255)
                                                                    cATCCGCGCTCGCCTGCCGCACAGCAACGCACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTgg                                                                             <  1:1042463/149‑1 (MQ=255)
                                                                      tCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCa                                                                           <  2:1507104/149‑1 (MQ=255)
                                                                      tCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCa                                                                           >  2:166935/1‑149 (MQ=255)
                                                                      tCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGAGCCAGCAAACCGGTAGTGCATTACCAATGGTCGCCAGTCGCCGGTggaa                                                                           >  1:469823/1‑147 (MQ=255)
                                                                       ccGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATTGTCGCCAGTCGCCGGTGGCaa                                                                          <  2:945412/149‑1 (MQ=255)
                                                                        cGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAAt                                                                         <  1:497602/149‑1 (MQ=255)
                                                                                 cTGCCGCACCGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCAACAAAGCGCCAGCAAACCGGTTGTGCATTAACAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTg                                                                >  2:1068351/1‑149 (MQ=255)
                                                                                 cTGCCCCAAAGCAACACACACAGTCATTTATTGAAAAAGTCGCCGTGATTAACATCGGATTAAGCAGATTTGCGCAGGAGCTAAAAATCGGCAGCAAACCGGTTGTGCATTACCAATGGTAGACAGTCGCCGGTGGCAATTAAAAACTg                                                                >  2:1367256/1‑149 (MQ=255)
                                                                                  tGCCGCACAGCAACTCACACCGTGATTACTGGATAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTGGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTgg                                                               <  2:539691/149‑1 (MQ=255)
                                                                                  tGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTgg                                                               <  2:1378736/149‑1 (MQ=255)
                                                                                  tGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTgg                                                               <  2:655469/149‑1 (MQ=255)
                                                                                     cGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTc                                                            >  2:1069340/1‑149 (MQ=255)
                                                                                        acaGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGtt                                                         <  2:1198368/149‑1 (MQ=255)
                                                                                         caGCAACACACACCGTCATTACTGGAACACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGttt                                                        <  1:85747/149‑1 (MQ=255)
                                                                                          aGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTtta                                                       <  2:13178/149‑1 (MQ=255)
                                                                                           gCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTttat                                                      <  2:1381044/149‑1 (MQ=255)
                                                                                             aaCACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTttatta                                                    <  2:1281396/149‑1 (MQ=255)
                                                                                             aaCACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGAGTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTttatta                                                    <  1:760104/149‑1 (MQ=255)
                                                                                                  acacCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACg                                               >  1:1366347/1‑149 (MQ=255)
                                                                                                    acCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGtt                                             >  2:383429/1‑149 (MQ=255)
                                                                                                     ccGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGttt                                            <  1:202613/149‑1 (MQ=255)
                                                                                                      cGTCATTACTGGAAAACGTCTCCGTGATTATCATCGGATTACGCTGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTg                                           <  2:838631/149‑1 (MQ=255)
                                                                                                      cGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTg                                           <  2:64048/149‑1 (MQ=255)
                                                                                                      cGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTg                                           >  1:1167221/1‑149 (MQ=255)
                                                                                                       gTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGc                                          >  1:1510030/1‑149 (MQ=255)
                                                                                                           ttACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATa                                      >  2:1266418/1‑149 (MQ=255)
                                                                                                             aCTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAg                                    <  1:217260/149‑1 (MQ=255)
                                                                                                                 gAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGt                                >  2:965848/1‑149 (MQ=255)
                                                                                                                     aCGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTccc                            <  2:925538/149‑1 (MQ=255)
                                                                                                                     aCGTCGCCGTGATTAACATCGGATTACGCAGATTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTccc                            <  1:387584/149‑1 (MQ=255)
                                                                                                                     aCGTCGCCGTGATTAACATCGGATTACGCAGATTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTccc                            <  1:1176645/149‑1 (MQ=255)
                                                                                                                       gTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCAt                          <  2:176814/149‑1 (MQ=255)
                                                                                                                       gTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCAt                          <  2:159629/149‑1 (MQ=255)
                                                                                                                        tCGCCGTGATTAACATCTGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATg                         <  1:1115790/149‑1 (MQ=255)
                                                                                                                               gATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACAt                  <  2:235632/149‑1 (MQ=255)
                                                                                                                                     cATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTgg            >  1:620134/1‑149 (MQ=255)
                                                                                                                                      aTCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAACACTGGCTCGTTTATTAGAACGTTTGAAATAAGGGGTTACCATGTTTACATCAGTgga           >  1:1240990/1‑148 (MQ=255)
                                                                                                                                          gATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCaa       >  1:677566/1‑149 (MQ=255)
                                                                                                                                          gATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCaa       >  1:671672/1‑149 (MQ=255)
                                                                                                                                          gATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCAGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCAGGTTTACATCAGTGGCGCaa       >  1:107339/1‑149 (MQ=255)
                                                                                                                                               cGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCaa  <  2:590830/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTTGTGCATTACCAATGGTCGCCAGTCGCCGGTGGCAATAAAAAACTGGCTCGTTTATTAGAACGTTTGCAATAAGGGGTTCCCATGTTTACATCAGTGGCGCAAGCCAA  >  NC_000913/548091‑548382

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