Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,333,863 A→G R91R (CGT→CGC murA ← UDP‑N‑glucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,8630AG99.0% 713.9 / ‑2.8 207R91R (CGT→CGCmurAUDP‑N‑glucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (97/108);  total (99/108)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 214 total aligned reads displayed.

GGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCG  >  NC_000913/3333774‑3333951
                                                                                         |                                                                                        
 gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCa                                                                                      >  1:963989/1‑93 (MQ=255)
 gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCa                                                                                      <  1:1001378/93‑1 (MQ=255)
 gACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCa                                                                                      <  2:2373751/93‑1 (MQ=255)
    gcgcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATgg                                                                                   <  1:5438117/93‑1 (MQ=255)
     cgcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGt                                                                                  <  1:3715860/93‑1 (MQ=255)
      gcACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGtt                                                                                 <  2:273736/93‑1 (MQ=255)
          cGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTaa                                                                             >  2:6026136/1‑93 (MQ=255)
           gATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAc                                                                            >  2:2838565/1‑93 (MQ=255)
            aTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAcc                                                                           >  1:103138/1‑93 (MQ=255)
            aTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAcc                                                                           >  2:1446763/1‑93 (MQ=255)
            aTCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAAcc                                                                           >  1:58156/1‑93 (MQ=255)
             tCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCa                                                                          <  1:3887034/93‑1 (MQ=255)
                tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAt                                                                       >  2:5817912/1‑93 (MQ=255)
                tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAt                                                                       >  1:3849769/1‑93 (MQ=255)
                tACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGAt                                                                       >  2:2437608/1‑93 (MQ=255)
                 aCAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATc                                                                      <  2:1036004/93‑1 (MQ=255)
                   aaCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCgc                                                                    >  2:3787695/1‑92 (MQ=255)
                   aaCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGt                                                                    >  2:503048/1‑93 (MQ=255)
                    aCCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTa                                                                   <  1:1866733/93‑1 (MQ=255)
                     ccGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTaa                                                                  <  2:3497869/93‑1 (MQ=255)
                       gCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAgg                                                                >  2:2105401/1‑93 (MQ=255)
                        ccAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGt                                                               <  2:5756754/93‑1 (MQ=255)
                        ccAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGt                                                               <  1:2474320/93‑1 (MQ=255)
                         cAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTg                                                              >  2:4849654/1‑93 (MQ=255)
                         cAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTg                                                              >  2:3008396/1‑93 (MQ=255)
                          aGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTgc                                                             >  2:3607118/1‑93 (MQ=255)
                            gTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTgcgc                                                           >  2:6052039/1‑93 (MQ=255)
                            gTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTgcgc                                                           >  1:2802113/1‑93 (MQ=255)
                                tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                                       >  2:22876/1‑93 (MQ=255)
                                tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                                       >  2:1994070/1‑93 (MQ=255)
                                tGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGaa                                                       <  1:6147859/93‑1 (MQ=255)
                                 gAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAAt                                                      <  2:829880/93‑1 (MQ=255)
                                  aaaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATa                                                     <  1:1233953/93‑1 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGGGCATGGTTTTAACCAGATCGTAAGGTGCCAGAATACa                                                   <  1:5016049/93‑1 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                                                    <  1:1289283/93‑1 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                                                    >  1:3257492/1‑93 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATAc                                                    >  2:4223538/1‑93 (MQ=255)
                                   aaCTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATAc                                                    >  1:4861024/1‑93 (MQ=255)
                                     cTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACAt                                                  <  2:4564033/93‑1 (MQ=255)
                                      ttGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAtt                                                 >  2:2501931/1‑93 (MQ=255)
                                      ttGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACAtt                                                 >  2:69402/1‑93 (MQ=255)
                                      ttGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACAtt                                                 >  1:822672/1‑93 (MQ=255)
                                       tGCCCCTGACCAAAGCGCGCGACCAGCGGCCCCAGCGCCAAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTa                                                <  2:936742/93‑1 (MQ=255)
                                         ccccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAAc                                              <  2:4790932/93‑1 (MQ=255)
                                          cccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACg                                             >  1:1323963/1‑93 (MQ=255)
                                           ccTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACgc                                            >  2:2557483/1‑92 (MQ=255)
                                            cTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTc                                           <  2:4016695/93‑1 (MQ=255)
                                             tGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcg                                          >  2:4943109/1‑93 (MQ=255)
                                              gACCAAAGCGCGCTGCCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgc                                         <  1:79037/93‑1 (MQ=255)
                                              gACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgc                                         <  2:3259081/93‑1 (MQ=255)
                                               aCCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTcgcg                                        <  1:4463269/93‑1 (MQ=255)
                                                ccAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCgg                                       >  1:5390394/1‑93 (MQ=255)
                                                 cAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCggg                                      <  2:1115118/93‑1 (MQ=255)
                                                  aaaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGGAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGc                                     >  1:6083474/1‑93 (MQ=255)
                                                  aaaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGc                                     >  1:6104097/1‑93 (MQ=255)
                                                   aaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCa                                    <  2:410006/93‑1 (MQ=255)
                                                   aaGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCa                                    >  2:5348827/1‑93 (MQ=255)
                                                    aGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCAt                                   <  2:889067/93‑1 (MQ=255)
                                                     gcgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc                                  <  1:3991161/93‑1 (MQ=255)
                                                     gcgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc                                  >  1:4072685/1‑93 (MQ=255)
                                                     gcgcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATc                                  <  2:3060461/93‑1 (MQ=255)
                                                       gcgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCaa                                >  2:5700866/1‑93 (MQ=255)
                                                        cgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAat                               <  1:4226635/93‑1 (MQ=255)
                                                        cgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAat                               <  2:5656086/93‑1 (MQ=255)
                                                        cgcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAat                               <  2:5214811/93‑1 (MQ=255)
                                                         gcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAata                              >  2:5857858/1‑93 (MQ=255)
                                                         gcTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAata                              <  2:5932285/93‑1 (MQ=255)
                                                          cTACCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAatat                             >  2:1404557/1‑93 (MQ=255)
                                                            aCCAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGc                           <  1:517543/93‑1 (MQ=255)
                                                             ccAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGca                          >  2:3446349/1‑93 (MQ=255)
                                                              cAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcac                         <  1:546763/93‑1 (MQ=255)
                                                              cAGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcac                         <  2:2148568/93‑1 (MQ=255)
                                                               aGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcaca                        <  2:4474385/93‑1 (MQ=255)
                                                               aGCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGcaca                        >  2:1778828/1‑93 (MQ=255)
                                                                gCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAg                       >  1:2863651/1‑93 (MQ=255)
                                                                gCGGCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAg                       >  2:1693633/1‑93 (MQ=255)
                                                                  ggCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGaa                     <  1:3275499/93‑1 (MQ=255)
                                                                   gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAc                    >  2:1984592/1‑93 (MQ=255)
                                                                   gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAc                    >  1:607413/1‑93 (MQ=255)
                                                                   gCCCCAGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGGGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAAc                    >  1:5474168/1‑93 (MQ=255)
                                                                        aGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTa               <  1:3758420/93‑1 (MQ=255)
                                                                        aGCGCCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTa               <  2:2610200/93‑1 (MQ=255)
                                                                         gcgcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTAc              >  1:2637526/1‑93 (MQ=255)
                                                                           gcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGt            <  2:1679811/93‑1 (MQ=255)
                                                                           gcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGt            <  1:3312072/93‑1 (MQ=255)
                                                                           gcCCAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGt            >  1:4721841/1‑93 (MQ=255)
                                                                            cccAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGtt           >  1:1520656/1‑93 (MQ=255)
                                                                            cccAGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACGGAACCATTACGtt           >  1:1823768/1‑93 (MQ=255)
                                                                              cAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCt         >  1:4150114/1‑93 (MQ=255)
                                                                              cAGATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCt         >  1:2595568/1‑93 (MQ=255)
                                                                               aGATAGAAGCGCGCATGGTTTTAACCAGATCGCAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTa        <  2:3622230/93‑1 (MQ=255)
                                                                                gATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTAc       <  2:3273865/93‑1 (MQ=255)
                                                                                gATAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACAATTACGTTCTAc       >  2:1279163/1‑93 (MQ=255)
                                                                                 aTAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACt      <  2:2571072/93‑1 (MQ=255)
                                                                                  tAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     <  1:3029188/93‑1 (MQ=255)
                                                                                  tAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     >  2:4255927/1‑93 (MQ=255)
                                                                                  tAGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     <  1:914719/93‑1 (MQ=255)
                                                                                  aaGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACtt     <  1:2165060/92‑1 (MQ=255)
                                                                                   aGAAGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACttt    <  1:4970690/93‑1 (MQ=255)
                                                                                     aaGCGCGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCg  <  2:507639/93‑1 (MQ=255)
                                                                                         |                                                                                        
GGACGCGCACCGATCGTACAACCGCCAGGTAGTGAAACTTGCCCCTGACCAAAGCGCGCTACCAGCGGCCCCAGCGCCCAGATAGAAGCACGCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCG  >  NC_000913/3333774‑3333951

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