Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,900 G→A G113G (GGC→GGT yieO ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,9000GA99.6% 946.3 / ‑5.5 271G113G (GGC→GGTyieOputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (135/134);  total (136/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 294 total aligned reads displayed.

CGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACT  >  NC_000913/3937810‑3937989
                                                                                          |                                                                                         
cGGCATGTCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcacc                                                                                         <  1:1783518/93‑1 (MQ=255)
 ggCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcaccg                                                                                        >  2:3639921/1‑93 (MQ=255)
 ggCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcaccg                                                                                        <  1:1830999/93‑1 (MQ=255)
  gCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGcaccgc                                                                                       <  2:7751229/93‑1 (MQ=255)
   cATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCaccgcc                                                                                      <  2:7463886/93‑1 (MQ=255)
    aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCt                                                                                     <  1:6218602/93‑1 (MQ=255)
    aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCt                                                                                     >  1:4058917/1‑93 (MQ=255)
    aTGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCt                                                                                     <  2:4814402/93‑1 (MQ=255)
      ggCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTAt                                                                                   >  2:4891982/1‑93 (MQ=255)
      ggCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTAt                                                                                   <  1:8028005/93‑1 (MQ=255)
      ggCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTAt                                                                                   <  1:774367/93‑1 (MQ=255)
       gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATc                                                                                  >  2:5824675/1‑93 (MQ=255)
       gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATc                                                                                  <  2:3320105/93‑1 (MQ=255)
        cGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcc                                                                                 >  1:8540678/1‑93 (MQ=255)
         gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATccc                                                                                >  1:694237/1‑93 (MQ=255)
         gACAAAATTCAATACTGGAAGAAGTTCATTACGAGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATcccc                                                                               <  1:793696/93‑1 (MQ=255)
           cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCt                                                                              >  2:3848762/1‑93 (MQ=255)
           cAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCt                                                                              <  1:6413718/93‑1 (MQ=255)
             aaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGa                                                                            >  2:4245425/1‑93 (MQ=255)
               aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAAt                                                                          <  2:1243566/93‑1 (MQ=255)
                ttCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATa                                                                         >  1:1211343/1‑93 (MQ=255)
                 tCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATaa                                                                        >  1:6338799/1‑93 (MQ=255)
                   aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAcc                                                                      >  1:4748202/1‑93 (MQ=255)
                   aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAcc                                                                      <  1:2610803/93‑1 (MQ=255)
                    aTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAccc                                                                     >  2:2833880/1‑93 (MQ=255)
                    aTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAAccc                                                                     <  2:4373619/93‑1 (MQ=255)
                       cTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGa                                                                  <  1:1013020/93‑1 (MQ=255)
                       cTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGa                                                                  <  2:7224552/93‑1 (MQ=255)
                        tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGaa                                                                 <  2:4450973/93‑1 (MQ=255)
                        tGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGaa                                                                 >  2:5281419/1‑93 (MQ=255)
                           aagaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGAc                                                              >  1:8687571/1‑93 (MQ=255)
                            agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACa                                                             >  1:1270891/1‑93 (MQ=255)
                             gaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACaa                                                            <  1:878907/93‑1 (MQ=255)
                             gaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACaa                                                            >  1:8379939/1‑93 (MQ=255)
                                gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCa                                                         >  1:3817711/1‑93 (MQ=255)
                                gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCa                                                         >  2:2667778/1‑93 (MQ=255)
                                gTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGAATCATCATTGCAACGCCTATCCCCTGAATAACCCGGAAGACAACCa                                                         >  1:5391022/1‑93 (MQ=255)
                                 ttCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAg                                                        <  2:2413805/93‑1 (MQ=255)
                                  tCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGc                                                       >  1:3917937/1‑93 (MQ=255)
                                  tCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGc                                                       >  2:7256798/1‑93 (MQ=255)
                                      tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtgtg                                                   >  2:3912168/1‑93 (MQ=255)
                                      tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCtgtg                                                   <  2:6750630/93‑1 (MQ=255)
                                       aCGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAATCCGGAAGACAACCAGCTGTgg                                                  >  2:233064/1‑93 (MQ=255)
                                       aCGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTgg                                                  >  1:6627741/1‑93 (MQ=255)
                                       aCGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTgg                                                  >  1:8655174/1‑93 (MQ=255)
                                        cGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGt                                                 >  1:2450432/1‑93 (MQ=255)
                                        cGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGt                                                 <  2:399856/93‑1 (MQ=255)
                                        cGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGt                                                 <  1:5634328/93‑1 (MQ=255)
                                         gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTa                                                >  1:227202/1‑93 (MQ=255)
                                          aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAg                                               <  1:4138370/93‑1 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  2:8344730/1‑93 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  1:7339130/1‑93 (MQ=255)
                                           ggATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGc                                              >  1:7596519/1‑93 (MQ=255)
                                             aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGa                                            <  2:4749969/93‑1 (MQ=255)
                                              tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGat                                           >  2:7858227/1‑92 (MQ=255)
                                              tAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGaa                                           >  2:601977/1‑93 (MQ=255)
                                               aaGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAt                                          <  2:322980/93‑1 (MQ=255)
                                                aGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCCGCTGTGGTAGCGAAtt                                         >  1:4436016/1‑93 (MQ=255)
                                                    cgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGaa                                     <  2:3014809/93‑1 (MQ=255)
                                                       aGTAAGGCCAGCCGAGCAACAGGCATGATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                  <  2:1036954/93‑1 (MQ=255)
                                                       aGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                  >  2:1599554/1‑93 (MQ=255)
                                                       aGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                  <  1:5722098/93‑1 (MQ=255)
                                                       aGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                  <  2:3725191/93‑1 (MQ=255)
                                                       aGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGt                                  <  1:7865271/93‑1 (MQ=255)
                                                        gTAAGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTg                                 >  1:4898890/1‑93 (MQ=255)
                                                          aaGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcg                               <  2:2089093/93‑1 (MQ=255)
                                                          aaGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcg                               >  1:2763585/1‑93 (MQ=255)
                                                           aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              >  2:8167516/1‑93 (MQ=255)
                                                           aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              >  1:5775416/1‑93 (MQ=255)
                                                           aGGCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTgcgc                              <  1:5104114/93‑1 (MQ=255)
                                                            ggCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCa                             <  2:617957/93‑1 (MQ=255)
                                                            ggCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCa                             <  2:535744/93‑1 (MQ=255)
                                                            ggCCAGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCa                             <  1:2074839/93‑1 (MQ=255)
                                                                aGCCGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGcc                         <  2:7970891/93‑1 (MQ=255)
                                                                   cGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                      <  1:6010983/93‑1 (MQ=255)
                                                                   cGAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                      >  1:1439426/1‑93 (MQ=255)
                                                                    gAGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCagag                     >  1:4728280/1‑93 (MQ=255)
                                                                     aGCAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAgaga                    <  2:5057791/93‑1 (MQ=255)
                                                                     aGCAACAGGCATCATAATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGCGCGCAGGCCAgaga                    >  1:7188711/1‑93 (MQ=255)
                                                                       cAACAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAc                  <  2:1968961/93‑1 (MQ=255)
                                                                         aCAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                >  2:5504380/1‑93 (MQ=255)
                                                                         aCAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                >  2:5725236/1‑93 (MQ=255)
                                                                         aCAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                >  1:2953753/1‑93 (MQ=255)
                                                                          cAGGCATCATCATTGCACCGCCTATCCCCTGAATACCCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCa               <  2:1532648/93‑1 (MQ=255)
                                                                          cAGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCa               <  2:8732748/93‑1 (MQ=255)
                                                                           aGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCaa              <  2:6808090/93‑1 (MQ=255)
                                                                           aGGCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCaa              <  1:7500512/93‑1 (MQ=255)
                                                                            ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  2:4533814/1‑93 (MQ=255)
                                                                            ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             >  2:6897165/1‑93 (MQ=255)
                                                                            ggCATCATCATTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAt             <  1:8423770/93‑1 (MQ=255)
                                                                              catcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgt           <  1:1313212/93‑1 (MQ=255)
                                                                              catcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgt           <  1:3200092/93‑1 (MQ=255)
                                                                               atcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgtg          <  2:6611226/93‑1 (MQ=255)
                                                                                tcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGTAAGTGCGCAGGCCAGAGAACCCAATGTGa         >  1:3748913/1‑93 (MQ=255)
                                                                                tcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGa         >  1:5479090/1‑93 (MQ=255)
                                                                                tcatcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGa         >  1:7077692/1‑93 (MQ=255)
                                                                                 catcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa        <  1:6526753/93‑1 (MQ=255)
                                                                                 catcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCACCGGGGGCAGCGAATTAGAAAGTGCGCAGGCAAGAGAACCAAATGTGaa        >  2:1091680/1‑93 (MQ=255)
                                                                                 catcatTGCACCGCCTAGCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAAGTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGaa        <  1:1421561/93‑1 (MQ=255)
                                                                                   tcatTGCACCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACa      <  2:7627000/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACT  >  NC_000913/3937810‑3937989

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