Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,331,998 T→G D499E (GAT→GAG)  proP → proline/glycine betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,9980TG93.2% 354.7 / 13.8 119D499E (GAT→GAG) proPproline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (49/61);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 121 total aligned reads displayed.

ACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCC  >  NC_000913/4331901‑4332094
                                                                                                 |                                                                                                
 cGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGaa                                                                                               >  2:3970022/1‑100 (MQ=255)
 cGATAAGATCGAGCAGAAAAGCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGaa                                                                                               <  2:2666077/100‑1 (MQ=255)
   ataataTCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATa                                                                                             <  1:1424836/100‑1 (MQ=255)
   ataataTCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGAATTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATa                                                                                             >  2:1696436/1‑100 (MQ=255)
     aataTCGAGCAGAAAACCGATGATATTGACCACGAGATTGCCGATTTGCAGTCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAg                                                                                           <  2:1868726/100‑1 (MQ=255)
      atatCGGGCTTAAAATCGATGAGATTTACCACGAGATTGCTGCTTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATACGCGAATTGAGGAATGAGc                                                                                          <  2:259868/100‑1 (MQ=255)
      atatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGc                                                                                          <  1:197380/100‑1 (MQ=255)
       tatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCt                                                                                         <  2:2273281/100‑1 (MQ=255)
        atCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTg                                                                                        <  1:984958/100‑1 (MQ=255)
         tCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGAGCAGCAACATCCGCGAATTGAGGAATAAGCTGa                                                                                       >  2:149956/1‑100 (MQ=255)
         tCGAGCAGAAAATAGATGATATTGACCACGAGATTGCCGATTTGCACGCGAAACGTACACGCCTGTTGCAGCAACATCAGCGAATTGAGGAACAAgagga                                                                                       >  2:3718784/1‑96 (MQ=255)
          cGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGaa                                                                                      <  1:1648519/100‑1 (MQ=255)
          cGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGaa                                                                                      >  1:3034224/1‑100 (MQ=255)
           gAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTCCCCGCCTGGTCCAGCAAGATCCGCGACTTGACGAATAAGCTGaaa                                                                                     >  2:3819859/1‑100 (MQ=255)
           gAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGaaa                                                                                     >  1:543527/1‑100 (MQ=255)
            aGCAGAAAATCGATGATATTGACCACGAGATTGCCGAATTGCAGGCGAAACGTACACGCCTGGTGCAGCAACCTACGCGAATTGCGGAATAAGCTGacac                                                                                    >  2:1573748/1‑100 (MQ=255)
              cAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACgg                                                                                  >  2:3065718/1‑100 (MQ=255)
               aGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGa                                                                                 <  2:1502698/100‑1 (MQ=255)
                  aaaTCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATgg                                                                              <  1:3928670/100‑1 (MQ=255)
                  aaaTCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATgg                                                                              <  1:526029/100‑1 (MQ=255)
                   aaTCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGc                                                                             <  1:642373/100‑1 (MQ=255)
                    aTCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGcc                                                                            <  2:1691439/100‑1 (MQ=255)
                     tCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCt                                                                           <  1:1359589/100‑1 (MQ=255)
                       gatgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGa                                                                         >  1:4045177/1‑100 (MQ=255)
                       gatgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGa                                                                         <  2:1472752/100‑1 (MQ=255)
                        atgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGAt                                                                        <  2:3753906/100‑1 (MQ=255)
                        atgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGAt                                                                        >  1:2372242/1‑100 (MQ=255)
                         tgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGAtg                                                                       <  1:1788836/100‑1 (MQ=255)
                           atatTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGAtgtg                                                                     <  2:1095057/100‑1 (MQ=255)
                             atTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGAc                                                                   <  1:1604244/100‑1 (MQ=255)
                             atTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGAc                                                                   >  1:2521700/1‑100 (MQ=255)
                               tGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGc                                                                 <  1:2518160/100‑1 (MQ=255)
                                gACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCt                                                                >  1:3625634/1‑100 (MQ=255)
                                  ccACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGt                                                              >  1:1731969/1‑100 (MQ=255)
                                   cACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTc                                                             >  1:3439771/1‑100 (MQ=255)
                                   cACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTc                                                             >  1:3212612/1‑100 (MQ=255)
                                    aCGAGATTGCCGATTTGCAGGCGAAACGTACCCGCATGGTGCAGCAACATCCGCGAATTGAGGAATGAGCTGAAACGGATGGCCTGATGTGACGCTGTCt                                                            <  2:3949040/100‑1 (MQ=255)
                                     cGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCtt                                                           <  2:412471/100‑1 (MQ=255)
                                       agaTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTAt                                                         >  2:531512/1‑100 (MQ=255)
                                        gaTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATc                                                        <  1:515854/100‑1 (MQ=255)
                                        gaTTGCCGAATTGCAGGTGAAACGTACCCGCCTAGTGCCGCAGCATACGCGAATTGGGGAATAAGCTGCTACGCATGGCCTGATGTGACGCCGTCTTacc                                                        <  2:473596/100‑3 (MQ=255)
                                         aTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCa                                                       <  1:2678787/100‑1 (MQ=255)
                                          ttGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAg                                                      <  2:2538313/100‑1 (MQ=255)
                                           tGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAgg                                                     <  2:1915026/100‑1 (MQ=255)
                                           tGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAgg                                                     <  1:1146189/100‑1 (MQ=255)
                                           ggCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAgg                                                     <  2:1857181/99‑1 (MQ=255)
                                            gCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGc                                                    >  2:1665970/1‑100 (MQ=255)
                                            gCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGc                                                    <  2:3437005/100‑1 (MQ=255)
                                            gCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGc                                                    >  1:2503408/1‑100 (MQ=255)
                                             ccGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGcc                                                   <  2:2580739/100‑1 (MQ=255)
                                             ccGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGcc                                                   >  1:2367821/1‑100 (MQ=255)
                                             ccGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGcc                                                   <  1:2595929/100‑1 (MQ=255)
                                              cGATTTGCAGGGGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCa                                                  <  2:3198279/100‑1 (MQ=255)
                                              cGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCa                                                  <  2:2130154/100‑1 (MQ=255)
                                                 tttGCAGGCTAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGATACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAAtt                                               <  1:3561798/100‑1 (MQ=255)
                                                 tttGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAAtt                                               <  2:627407/100‑1 (MQ=255)
                                                  ttGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTg                                              >  2:60148/1‑100 (MQ=255)
                                                  ttGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTg                                              <  2:661051/100‑1 (MQ=255)
                                                   tGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGa                                             >  2:111143/1‑100 (MQ=255)
                                                    gCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAAGTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTCTCAGGCCACTTGaa                                            >  2:3438364/1‑100 (MQ=255)
                                                     cAGGTGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAc                                           <  1:1044720/100‑1 (MQ=255)
                                                     cAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAc                                           <  1:1085160/100‑1 (MQ=255)
                                                     cAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAc                                           >  1:3559570/1‑100 (MQ=255)
                                                     cAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATCGCCTGATGTGACGCTGTCTTCTCAGGCCAATTGAAc                                           >  2:2752602/1‑100 (MQ=255)
                                                      aGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAct                                          <  2:2840322/100‑1 (MQ=255)
                                                       ggCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAActc                                         >  1:2555759/1‑100 (MQ=255)
                                                        gCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGATACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAActct                                        >  2:3694921/1‑100 (MQ=255)
                                                          gAAAGGTACCCGCCTGGGGCAGTAACATCCGCGATGTTATTAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTa                                      <  1:4078299/100‑1 (MQ=255)
                                                           aaaCGTACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTaa                                     <  1:2925700/100‑1 (MQ=255)
                                                            aaCGTCACCGCCTGGTGCAGCAACACCCGCAAATTGAGGAATAAGCTGAAACGGATGGCGTGATGTGACGCTGTCTGATCAGGCCAATTGAACTCTTAAg                                    <  2:1916541/100‑1 (MQ=255)
                                                                tACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTc                                >  2:3667748/1‑100 (MQ=255)
                                                                tACCCGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATGGAACTCTTAAGGTTc                                <  2:3894559/100‑1 (MQ=255)
                                                                  cccGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCAc                              >  2:2427981/1‑100 (MQ=255)
                                                                  cccGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCAc                              >  2:3049020/1‑100 (MQ=255)
                                                                   ccGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACt                             <  1:1108090/100‑1 (MQ=255)
                                                                   ccGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACt                             <  1:328899/100‑1 (MQ=255)
                                                                    cGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACtt                            <  2:3482613/100‑1 (MQ=255)
                                                                    cGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACtt                            >  1:1809001/1‑100 (MQ=255)
                                                                    cGCCTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACtt                            >  1:3400974/1‑100 (MQ=255)
                                                                       cTGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAt                         >  2:1587233/1‑100 (MQ=255)
                                                                        tGGTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtc                        >  2:1049751/1‑100 (MQ=255)
                                                                         ggTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtct                       >  2:1055170/1‑100 (MQ=255)
                                                                         ggTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtct                       >  1:1059604/1‑100 (MQ=255)
                                                                          gTGCAGCAACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtctc                      >  1:2852557/1‑100 (MQ=255)
                                                                             cagcaACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGa                   >  1:3171325/1‑100 (MQ=255)
                                                                             cagcaACATCCGCGAATTGAGGAATAAGCTGAAACAGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGa                   >  1:3944338/1‑100 (MQ=255)
                                                                              agcaACATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGAc                  <  1:513722/100‑1 (MQ=255)
                                                                                 aaCATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGAcgcg               <  1:166770/100‑1 (MQ=255)
                                                                                   cATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCa             >  2:2492516/1‑100 (MQ=255)
                                                                                   cATCCGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCa             <  2:3271524/100‑1 (MQ=255)
                                                                                      ccGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAc          <  1:3276250/100‑1 (MQ=255)
                                                                                      ccGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAc          <  1:1923069/100‑1 (MQ=255)
                                                                                      ccGCGAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAc          <  1:2385547/100‑1 (MQ=255)
                                                                                       cgcgAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTAAACTCTTAAGGTTCAGTTAATCTCTGACGCGCATAct         >  1:3799318/1‑100 (MQ=255)
                                                                                       cgcaAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAct         >  2:1828728/1‑100 (MQ=255)
                                                                                         cgAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATActct       <  1:2744124/100‑1 (MQ=255)
                                                                                          gAATTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACtctc      >  2:3535371/1‑100 (MQ=255)
                                                                                           aaTTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcc     >  2:3162801/1‑100 (MQ=255)
                                                                                           aaTTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcc     <  1:2470697/100‑1 (MQ=255)
                                                                                              tgaggaATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCtcctcc  >  2:2639397/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCC  >  NC_000913/4331901‑4332094

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