Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,090,226 T→C 6.9% C158R (TGC→CGC)  yggI → Zn‑dependent metalloprotease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,2260TC 6.9% 376.0 / 13.9 131C158R (TGC→CGC) yggIZn‑dependent metalloprotease‑related protein
Reads supporting (aligned to +/- strand):  ref base T (63/58);  new base C (5/4);  total (69/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 132 total aligned reads displayed.

TGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGT  >  NC_000913/3090078‑3090375
                                                                                                                                                    |                                                                                                                                                     
   tCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCgg                                                                                                                                                   >  1:330697/1‑150 (MQ=255)
      cgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGCGGCCGTCTATCGCTGTGTTCACTGCGGTGa                                                                                                                                                >  2:44692/1‑150 (MQ=255)
      cgcccgTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGa                                                                                                                                                >  2:499522/1‑150 (MQ=255)
            tcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACAGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATAGCTGTGTTCACTGCGGTGAACAGCt                                                                                                                                          >  1:898599/1‑150 (MQ=255)
            gcgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCt                                                                                                                                          <  1:447439/149‑1 (MQ=255)
             cgtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACacagc                                                                                                                                         <  1:977622/150‑5 (MQ=255)
              gtACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTCCCGCTGCAAGTGCCAGGAGCAGCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACCGCAGTGAACacagcg                                                                                                                                        <  1:843792/150‑6 (MQ=255)
                  gCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTCCCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGc                                                                                                                                    <  2:1331733/150‑1 (MQ=255)
                  gCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGc                                                                                                                                    >  2:122907/1‑150 (MQ=255)
                   cATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCg                                                                                                                                   <  1:1534769/150‑1 (MQ=255)
                    aTCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGa                                                                                                                                  <  1:772665/150‑1 (MQ=255)
                    aTCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGa                                                                                                                                  <  1:1566286/150‑1 (MQ=255)
                      cAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGaaa                                                                                                                                <  1:1407571/150‑1 (MQ=255)
                       aGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAAt                                                                                                                               >  2:1390800/1‑150 (MQ=255)
                       aGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAAt                                                                                                                               >  2:286952/1‑150 (MQ=255)
                        gTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACGGCGGTGAACCGCTGGTTGCGAAATa                                                                                                                              >  1:1060939/1‑150 (MQ=255)
                          tCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATagc                                                                                                                            >  2:1616604/1‑148 (MQ=255)
                          tCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAc                                                                                                                            >  2:557288/1‑150 (MQ=255)
                           cGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAAcc                                                                                                                           >  1:1061742/1‑150 (MQ=255)
                            gAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          <  2:1449105/150‑1 (MQ=255)
                            gAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          <  1:528676/150‑1 (MQ=255)
                            gAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          >  1:1673488/1‑150 (MQ=255)
                            gAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCa                                                                                                                          >  1:1691538/1‑150 (MQ=255)
                                tGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTg                                                                                                                      >  1:1185403/1‑150 (MQ=255)
                                       cGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATc                                                                                                               >  2:1447812/1‑150 (MQ=255)
                                        gTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCa                                                                                                              >  1:1628460/1‑150 (MQ=255)
                                         tGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAg                                                                                                             <  2:938719/150‑1 (MQ=255)
                                         tGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAg                                                                                                             <  1:1326212/150‑1 (MQ=255)
                                           cGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGa                                                                                                           >  2:1231742/1‑150 (MQ=255)
                                            gTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCCCTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGaa                                                                                                          >  1:1271144/1‑150 (MQ=255)
                                            gTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGaa                                                                                                          >  1:1683998/1‑150 (MQ=255)
                                            gTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGaa                                                                                                          >  1:1049056/1‑150 (MQ=255)
                                             gcGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAAc                                                                                                         <  1:122907/149‑1 (MQ=255)
                                               gCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACtt                                                                                                       >  1:965013/1‑150 (MQ=255)
                                                cAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACttt                                                                                                      >  2:1195494/1‑150 (MQ=255)
                                                  acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                    <  2:1476814/150‑1 (MQ=255)
                                                  acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                    <  2:928269/150‑1 (MQ=255)
                                                  acacCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTcc                                                                                                    <  1:330969/150‑1 (MQ=255)
                                                  acacCTACCACTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCGGACATACTCGCGTAGTTCGGGGCGAGGACGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATTACCATCTGAACTATCAGGAACTTTcc                                                                                                    <  1:294474/150‑1 (MQ=255)
                                                    acCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCCCTGCGGTGAACAGCTGGTTGCGAAATCACCATCTGAACTATCAGGAACTTTCCTg                                                                                                  >  1:451761/1‑150 (MQ=255)
                                                    acCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTg                                                                                                  <  1:894520/150‑1 (MQ=255)
                                                     ccTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGa                                                                                                 >  1:636642/1‑150 (MQ=255)
                                                        tCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCt                                                                                              <  2:1472104/150‑1 (MQ=255)
                                                         ccccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTg                                                                                             <  1:691515/150‑1 (MQ=255)
                                                           ccTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAAAGTTCCTGAtagttc                                                                                           <  1:717987/150‑6 (MQ=255)
                                                               ccGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGGTCGTGGCGAGGCCGTCTATCGCTGTGTTCCCTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGAGCTGGCTGAt                                                                                       >  2:229456/1‑150 (MQ=255)
                                                                 gctgctAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTg                                                                                     <  1:368259/150‑1 (MQ=255)
                                                                 gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACATCTGGTTGCGAAATAACCATCTGAACTATCAGGAAATTTCGAGATCTGTCTGATTg                                                                                     >  2:190702/1‑150 (MQ=255)
                                                                 gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTg                                                                                     <  2:612954/150‑1 (MQ=255)
                                                                 gCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTg                                                                                     <  2:1628460/150‑1 (MQ=255)
                                                                  cTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGCACTTTCCTGATCTGGCTGATTGc                                                                                    >  2:14619/1‑150 (MQ=255)
                                                                  cTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    >  2:373291/1‑150 (MQ=255)
                                                                       aGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATAcc                                                                               <  2:330697/150‑1 (MQ=255)
                                                                         tGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCaa                                                                             <  1:1459384/150‑1 (MQ=255)
                                                                            cAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAAc                                                                          <  2:1061742/150‑1 (MQ=255)
                                                                             aGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACa                                                                         >  2:652379/1‑150 (MQ=255)
                                                                              ggAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAg                                                                        <  2:19021/150‑1 (MQ=255)
                                                                              ggAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAg                                                                        >  1:772335/1‑150 (MQ=255)
                                                                              ggAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAg                                                                        >  2:493950/1‑150 (MQ=255)
                                                                              ggAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCATGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    <  1:705987/138‑1 (MQ=255)
                                                                              ggAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCATGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGc                                                                                    >  2:705987/1‑138 (MQ=255)
                                                                              ggAGCATCAGCTTACCGTACGACGCCTTAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTTCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAg                                                                        <  1:821985/150‑1 (MQ=255)
                                                                                  cATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCttt                                                                    <  2:9642/150‑1 (MQ=255)
                                                                                  cATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCttt                                                                    >  2:1433804/1‑150 (MQ=255)
                                                                                    tCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCg                                                                  <  2:175647/150‑1 (MQ=255)
                                                                                         ttACCTTACGCCGCCAACAACGCGTAGTTCGGGGCGAGGCCGGCTAGCGCTGTGTTCACTGCGGTGAACAGCGGGTTGCGAAATACCCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACCGCTTTCGCTACg                                                             <  2:254604/150‑1 (MQ=255)
                                                                                          tACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGGCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt                                                            >  2:1375317/1‑150 (MQ=255)
                                                                                           aCCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt                                                           <  2:965013/150‑1 (MQ=255)
                                                                                            ccGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg                                                          >  1:27229/1‑150 (MQ=255)
                                                                                             cGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc                                                         >  2:539144/1‑150 (MQ=255)
                                                                                               tACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCTCTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg                                                       <  1:1390800/150‑1 (MQ=255)
                                                                                               tACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg                                                       <  2:1204150/150‑1 (MQ=255)
                                                                                               tACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg                                                       <  1:286952/150‑1 (MQ=255)
                                                                                                 cgccgcCATTATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCGGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTACTGGc                                                     <  1:1150875/150‑1 (MQ=255)
                                                                                                  gccgccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt                                                    <  1:1440047/150‑1 (MQ=255)
                                                                                                  gccgccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt                                                    <  2:1201522/150‑1 (MQ=255)
                                                                                                  gccgccATAATCGCGTAGGTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt                                                    <  1:88864/150‑1 (MQ=255)
                                                                                                   ccgccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc                                                   >  2:1535719/1‑150 (MQ=255)
                                                                                                    cgccATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg                                                  >  1:1591495/1‑150 (MQ=255)
                                                                                                         tAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTa                                             <  2:675264/150‑1 (MQ=255)
                                                                                                          aaTCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGAGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTaa                                            >  1:151289/1‑150 (MQ=255)
                                                                                                            tCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGTACTTTCCTGATCTTGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAaca                                          <  2:1112262/150‑1 (MQ=255)
                                                                                                                 tAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAg                                     <  1:1169309/150‑1 (MQ=255)
                                                                                                                     tCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAg                                 <  2:1591495/150‑1 (MQ=255)
                                                                                                                         ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGTTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGAt                             <  2:154312/150‑1 (MQ=255)
                                                                                                                         ggCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGAt                             >  1:1346397/1‑150 (MQ=255)
                                                                                                                            gAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTa                          >  1:1409919/1‑150 (MQ=255)
                                                                                                                             aGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTAc                         <  2:1549541/150‑1 (MQ=255)
                                                                                                                                ccGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGGTGCGAAATAGCCATCGGCACTATCAGGAACTTTCCTGCTATGGCTGATTGCATACCAAAACAGCTTACGCTACGTGGCTGGCTCGTTTTAACACGGAGTAAGGGATGAACCGt                      >  1:831289/1‑150 (MQ=255)
                                                                                                                                    cTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGGAATAACCATCTGAACTATCAGTAACTTTCCTGATCTGGCTGATTGAATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAtt                  <  1:866452/150‑1 (MQ=255)
                                                                                                                                    cTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCACTACGTTGCTGGCGCGTTTTAACACGGAGTAAGTGATGTACAGTTAtt                  >  2:871431/1‑150 (MQ=255)
                                                                                                                                    atATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTAtt                  <  1:1198003/149‑1 (MQ=255)
                                                                                                                                      aTCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTg                <  1:1282865/150‑1 (MQ=255)
                                                                                                                                        ccagctgtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTc              >  1:228812/6‑150 (MQ=255)
                                                                                                                                        ccagctgtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTc              >  1:410479/6‑150 (MQ=255)
                                                                                                                                         gCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt             >  2:1228219/1‑150 (MQ=255)
                                                                                                                                         cagctgtTCACCGCAGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCt             >  1:510411/5‑150 (MQ=255)
                                                                                                                                            gtgtTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTAtt          >  2:980508/1‑150 (MQ=255)
                                                                                                                                                  aCTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTACTGATCTGGCTGATTGCATACCAAAACAGATTTCGCTACGTAGCTGGCTATTTTTAACACCGAGTAAGTGATGAACCGTTATTTATCTATttctgag    >  2:302690/1‑148 (MQ=255)
                                                                                                                                                    tGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCggt  <  2:1683998/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTTCCCGCCCGTCGTACGCATCAGTTCGAACTGCAATCCGTGCGTCGCAACACCTTCCCCTACCGCTGCAAGTGCCAGGAGCATCAGCTTACCGTACGCCGCCATAATCGCGTAGTTCGTGGCGAGGCCGTCTATCGCTGTGTTCACTGCGGTGAACAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGT  >  NC_000913/3090078‑3090375

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