Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,335,317 T→G Y393S (TAT→TCT)  murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,335,3170TG96.4% 337.0 / 2.0 84Y393S (TAT→TCT) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (36/45);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CACTAAACCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGATAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCGATTTCGGC  >  NC_000913/3335169‑3335462
                                                                                                                                                    |                                                                                                                                                 
cacTAAACCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGaa                                                                                                                                                 <  1:1217316/151‑1 (MQ=255)
    aaaCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAAt                                                                                                                                               <  4:1125064/149‑1 (MQ=255)
    aaaCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAAt                                                                                                                                               <  4:462094/149‑1 (MQ=255)
    aaaCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATAc                                                                                                                                             >  1:798738/1‑151 (MQ=255)
      aCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACg                                                                                                                                            >  1:1000952/1‑150 (MQ=255)
      aCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACg                                                                                                                                            >  1:800928/1‑150 (MQ=255)
      aCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGa                                                                                                                                           <  1:1206892/151‑1 (MQ=255)
      aCCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGa                                                                                                                                           <  4:653529/151‑1 (MQ=255)
       ccGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGAt                                                                                                                                          >  1:1232484/1‑151 (MQ=255)
        tgAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATc                                                                                                                                         >  4:291127/2‑151 (MQ=255)
            aaaaTGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAacc                                                                                                                                     >  4:355262/1‑151 (MQ=255)
               aTGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAaccacc                                                                                                                                  >  4:559862/1‑151 (MQ=255)
                 gAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACcgt                                                                                                                                >  1:1113774/1‑151 (MQ=255)
                 gAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACcgt                                                                                                                                <  1:1115752/151‑1 (MQ=255)
                 gAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACcgt                                                                                                                                >  4:1261773/1‑151 (MQ=255)
                  aGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACcgtc                                                                                                                               <  1:638172/151‑1 (MQ=255)
                  aGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACcgtc                                                                                                                               >  1:431628/1‑151 (MQ=255)
                     ggTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCgtcgtc                                                                                                                            >  4:441523/1‑151 (MQ=255)
                          ggCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTc                                                                                                                       <  2:132688/151‑1 (MQ=255)
                           gCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTcc                                                                                                                      <  4:1122705/151‑1 (MQ=255)
                           gCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTcc                                                                                                                      <  1:509541/151‑1 (MQ=255)
                           gCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTcc                                                                                                                      <  1:64288/151‑1 (MQ=255)
                             ggTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGc                                                                                                                    >  1:1256923/1‑151 (MQ=255)
                                aGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAAt                                                                                                                 >  1:1229742/1‑151 (MQ=255)
                                     cgcgAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATAcagc                                                                                                            >  1:1331017/1‑151 (MQ=255)
                                       cgAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATAcagcca                                                                                                          <  4:259737/151‑1 (MQ=255)
                                        gAACGGTTCTGACAGCTCTCACACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCTTTCGTAGCCACGGTCGATGTGAGAAATACGATCAACCTCCGTCGTCCTTACCTCAATAcagcca                                                                                                          <  4:704294/150‑1 (MQ=255)
                                              ggCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCa                                                                                                   <  4:890940/151‑1 (MQ=255)
                                               gCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCAc                                                                                                  <  4:209443/151‑1 (MQ=255)
                                                cTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCAcc                                                                                                 >  4:830353/1‑151 (MQ=255)
                                                  gCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCAcc                                                                                                 >  3:160959/1‑149 (MQ=255)
                                                     aGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCt                                                                                            <  4:551174/151‑1 (MQ=255)
                                                           cAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTg                                                                                       >  1:750453/1‑150 (MQ=255)
                                                             gACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATg                                                                                    >  4:44653/1‑151 (MQ=255)
                                                             gACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATg                                                                                    <  4:881768/151‑1 (MQ=255)
                                                              aCGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGc                                                                                   <  1:181367/151‑1 (MQ=255)
                                                                 attattCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACg                                                                                >  4:505386/1‑151 (MQ=255)
                                                                    attCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAg                                                                             <  3:71189/151‑1 (MQ=255)
                                                                     ttCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGc                                                                               >  3:152340/1‑148 (MQ=255)
                                                                     ttCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGc                                                                               >  1:42628/1‑148 (MQ=255)
                                                                       cGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATc                                                                          <  1:386530/151‑1 (MQ=255)
                                                                       cGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATc                                                                          <  1:503172/151‑1 (MQ=255)
                                                                           tttCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGtt                                                                      >  1:927625/1‑151 (MQ=255)
                                                                            ttCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGtt                                                                      <  1:915455/150‑1 (MQ=255)
                                                                            ttCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGtt                                                                      <  4:902635/150‑1 (MQ=255)
                                                                            ttCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGt                                                                       <  4:785768/149‑1 (MQ=255)
                                                                              cacaCGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCgg                                                                        >  4:182028/1‑146 (MQ=255)
                                                                              cacaCGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGc                                                                    <  4:622761/150‑1 (MQ=255)
                                                                                   gCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAAc                                                              <  1:772490/151‑1 (MQ=255)
                                                                                   gCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAAc                                                              <  4:1113686/151‑1 (MQ=255)
                                                                                   gCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAAc                                                              >  4:1276844/1‑151 (MQ=255)
                                                                                      cAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCtg                                                           <  1:1208534/151‑1 (MQ=255)
                                                                                       aaTATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCtgt                                                          >  1:461265/1‑151 (MQ=255)
                                                                                        atatTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGATAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCtgtg                                                         <  1:170895/151‑1 (MQ=255)
                                                                                        atatTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCtgtg                                                         <  4:1046312/151‑1 (MQ=255)
                                                                                        atatTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCtgtg                                                         >  2:76448/1‑151 (MQ=255)
                                                                                          atTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTgcg                                                       >  1:1124700/1‑151 (MQ=255)
                                                                                             tGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGcc                                                     >  1:1205018/1‑150 (MQ=255)
                                                                                                  cTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGaa                                                 >  4:62911/1‑149 (MQ=255)
                                                                                                   tAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGt                                              >  4:429323/1‑151 (MQ=255)
                                                                                                   tAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGt                                              >  4:744912/1‑151 (MQ=255)
                                                                                                    aaaGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGtt                                             >  3:27543/1‑151 (MQ=255)
                                                                                                    aaaGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGt                                              <  2:78810/150‑1 (MQ=255)
                                                                                                     aaGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGttt                                            <  1:152218/151‑1 (MQ=255)
                                                                                                      aGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGtttt                                           <  3:8123/151‑1 (MQ=255)
                                                                                                        cgcgcAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGtttt                                           <  4:40319/149‑1 (MQ=255)
                                                                                                          cgcAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGATAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCa                                        >  1:451227/1‑150 (MQ=255)
                                                                                                            cAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGATAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAaca                                     <  4:211755/151‑1 (MQ=255)
                                                                                                              gTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACAcc                                   >  1:1060867/1‑151 (MQ=255)
                                                                                                                ttGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCg                                  <  4:956968/150‑1 (MQ=255)
                                                                                                                 tGTCTTCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGt                                 <  1:238086/150‑1 (MQ=255)
                                                                                                                     ttCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACa                             >  4:950554/1‑150 (MQ=255)
                                                                                                                      tCAATGCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACa                             <  1:622139/149‑1 (MQ=255)
                                                                                                                           gCGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACg                      <  4:842557/151‑1 (MQ=255)
                                                                                                                            cGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACgg                     >  4:1110116/1‑151 (MQ=255)
                                                                                                                            cGTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACgg                     <  4:985158/151‑1 (MQ=255)
                                                                                                                             gTTCGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGt                    <  2:10971/151‑1 (MQ=255)
                                                                                                                                cGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTAtt                 <  4:135982/151‑1 (MQ=255)
                                                                                                                                cGTAGCCACGATCGATGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTAtt                 <  4:798692/151‑1 (MQ=255)
                                                                                                                                        cgatcgatGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCg          >  1:45855/1‑150 (MQ=255)
                                                                                                                                        cgatcgatGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCGa         <  1:927687/151‑1 (MQ=255)
                                                                                                                                             gatGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCGATTTCg    <  2:133997/151‑1 (MQ=255)
                                                                                                                                             gatGTGAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCGATTTCg    <  4:826295/151‑1 (MQ=255)
                                                                                                                                                gtgAGAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCGATTTCGGc  <  4:1420836/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CACTAAACCGAAAAAATGAGTGGTCTGGCGGTAGCCCCGCGAACGGGGCTGCCAGCTCTCAGACGATTATTCGCCTTTCACACGCTCAATATTTGCACCTAAAGCGCGCAGTTTGTCTTCAATGCGTTCGTAGCCACGATCGATGTGATAAATACGATCAACCACCGTCGTCCCTTCCGCAATACAGCCAGCCAGCACCAGGCTTGCTGATGCACGCAGATCGGTTGCCATAACCTGTGCGCCAGAAAGTTTTTCAACACCGTGACAAATAACGGTATTGCTTTCGATTTCGGC  >  NC_000913/3335169‑3335462

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