Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,241 A→G intergenic (+19/‑164) cstA → / → ybdD carbon starvation protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,2410AG97.7% 668.6 / 7.4 213intergenic (+19/‑164)cstA/ybdDcarbon starvation protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (117/91);  minor base . (3/2);  total (120/93)
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.87e-01

Only 100 of 243 total aligned reads displayed.

AAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAT  >  NC_000913/631151‑631330
                                                                                          |                                                                                         
aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTg                                                                                         >  1:3974090/1‑93 (MQ=255)
aaagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTg                                                                                         >  1:2638126/1‑93 (MQ=255)
 aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGgtgt                                                                                        <  1:2536263/93‑1 (MQ=255)
  agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGCGAGGGGGGGGTGtt                                                                                       >  2:2552468/1‑93 (MQ=255)
   gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGGGCAGGCAAAAGGCGGACACTAAAGTCAGAGTGAGGAGGGGGTGTTg                                                                                      >  2:1678756/1‑93 (MQ=255)
   gaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTg                                                                                      <  1:1972719/93‑1 (MQ=255)
    aaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTgg                                                                                     <  1:6265433/93‑1 (MQ=255)
     aaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGc                                                                                    >  2:5847894/1‑93 (MQ=255)
      aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCg                                                                                   >  1:223303/1‑93 (MQ=255)
         ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAAt                                                                                <  1:834633/93‑1 (MQ=255)
         ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAAt                                                                                <  1:3656886/93‑1 (MQ=255)
         ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAAt                                                                                <  2:4666722/93‑1 (MQ=255)
           gTATGAACGAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAATTCAGAGTGAGGGGGGGGTGTTGGCGAATGt                                                                              <  2:3457883/93‑1 (MQ=255)
           gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGt                                                                              >  1:6697989/1‑93 (MQ=255)
           gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGGTGGCGAATGt                                                                              >  2:418282/1‑93 (MQ=255)
           gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGt                                                                              >  1:4572613/1‑93 (MQ=255)
            tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtt                                                                             >  2:4593802/1‑93 (MQ=255)
            tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtt                                                                             >  2:2474296/1‑93 (MQ=255)
             aTGAACCAATTCCGGAAAATGTCGAGGAGATCGTAGCGCAGGCAAAAGGCGCACACTACAGTCAGAGTGAAGGCGGAGAGCTGGCGAATGTTg                                                                            >  1:6007980/1‑93 (MQ=255)
             aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGAGAGGGGGGGGTGTTGGCGAATGCTg                                                                            >  2:3958023/1‑93 (MQ=255)
               gAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGc                                                                          >  2:6181283/1‑93 (MQ=255)
                aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCt                                                                         >  1:6519244/1‑93 (MQ=255)
                aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCt                                                                         <  2:6383382/93‑1 (MQ=255)
                 aCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCtt                                                                        <  2:739974/93‑1 (MQ=255)
                  ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTa                                                                       >  2:1568895/1‑93 (MQ=255)
                  ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTa                                                                       <  2:783784/93‑1 (MQ=255)
                    aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGt                                                                     <  1:3341029/93‑1 (MQ=255)
                    aaTGCCGGAAAATGTCGAGGAGATCCTGGCGCAGCCAAAAGGCGCACACTAAATTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGt                                                                     <  1:5541281/93‑1 (MQ=255)
                     aTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTg                                                                    >  1:4141033/1‑93 (MQ=255)
                       gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTTAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGcc                                                                  <  2:4323354/93‑1 (MQ=255)
                       gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGcc                                                                  <  2:2997317/93‑1 (MQ=255)
                        ccGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTAccc                                                                 <  1:6498314/93‑1 (MQ=255)
                         cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCa                                                                >  1:4804353/1‑93 (MQ=255)
                          ggAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAg                                                               <  2:5192284/93‑1 (MQ=255)
                          ggAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAg                                                               <  1:4908817/93‑1 (MQ=255)
                          ggAAAATGTCGAGGAGATCATGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAg                                                               <  2:2731685/93‑1 (MQ=255)
                           gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGATTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                              <  2:280013/93‑1 (MQ=255)
                           gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                              <  2:2246870/93‑1 (MQ=255)
                            aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                             >  1:928003/1‑93 (MQ=255)
                              aaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGtt                                                           <  2:5900409/93‑1 (MQ=255)
                                tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTcc                                                         <  1:2869021/93‑1 (MQ=255)
                                tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTcc                                                         >  1:3098516/1‑93 (MQ=255)
                                 gTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                        >  2:3270871/1‑93 (MQ=255)
                                 gTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                        <  2:1735228/93‑1 (MQ=255)
                                   cGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                                      >  1:1079036/1‑93 (MQ=255)
                                     aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                    <  1:6168753/93‑1 (MQ=255)
                                      ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                                   >  1:6883/1‑93 (MQ=255)
                                      ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                                   >  2:2970523/1‑93 (MQ=255)
                                        agaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAcc                                                 >  1:2896105/1‑93 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                <  2:2133286/93‑1 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                <  2:184714/93‑1 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                <  2:92511/93‑1 (MQ=255)
                                         gaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCAccc                                                >  1:6346099/1‑93 (MQ=255)
                                           tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                              >  1:3062574/1‑93 (MQ=255)
                                             gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc                                            <  1:5229413/93‑1 (MQ=255)
                                                    aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcc                                     <  2:1876615/93‑1 (MQ=255)
                                                    aGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTcc                                     >  1:4760112/1‑93 (MQ=255)
                                                     ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                                    >  1:395509/1‑93 (MQ=255)
                                                     ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                                    >  2:2447400/1‑93 (MQ=255)
                                                     ggCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTccc                                    >  1:4879973/1‑93 (MQ=255)
                                                        aaaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCgg                                 >  2:1443222/1‑93 (MQ=255)
                                                         aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGt                                >  2:3881170/1‑93 (MQ=255)
                                                          aaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTg                               >  2:861676/1‑93 (MQ=255)
                                                          aaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTg                               >  1:3466640/1‑93 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              <  2:1970920/93‑1 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              >  2:1387832/1‑93 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              <  2:4873611/93‑1 (MQ=255)
                                                           aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                              <  1:916168/93‑1 (MQ=255)
                                                            ggCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTggg                             >  2:5965034/1‑93 (MQ=255)
                                                             gcgcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgggg                            <  2:5433553/93‑1 (MQ=255)
                                                             gcgcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgggg                            <  2:3065176/93‑1 (MQ=255)
                                                               gcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACACTCTCCCCGGTGGGGCg                          >  1:1873884/1‑93 (MQ=255)
                                                                 acacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                        <  1:1359501/93‑1 (MQ=255)
                                                                  cacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                       <  2:1108237/93‑1 (MQ=255)
                                                                   acTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCCGTGGGGCGAggg                      >  1:1565273/1‑93 (MQ=255)
                                                                    cTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgggg                     >  2:3468179/1‑93 (MQ=255)
                                                                    cTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgggg                     >  2:581930/1‑93 (MQ=255)
                                                                        aGTCAGAGTGAGGGGGGGGTTTTGGCGAATGTTGGCTTAGTGCCCCGGGCTTCCTCTCACCCTCATCCTCTCCCCGGTGGGGTGACAGGACta                 >  1:3696454/1‑92 (MQ=255)
                                                                        aGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg                 >  2:2597895/1‑93 (MQ=255)
                                                                        aGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTg                 >  2:5749468/1‑93 (MQ=255)
                                                                          tCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc               <  1:2405846/93‑1 (MQ=255)
                                                                           cAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAcc              >  2:2903529/1‑93 (MQ=255)
                                                                           cAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAcc              <  2:5679020/93‑1 (MQ=255)
                                                                            agagTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg             >  1:1663582/1‑93 (MQ=255)
                                                                            agagTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCg             >  1:1399476/1‑93 (MQ=255)
                                                                                tGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcg         <  2:4487863/93‑1 (MQ=255)
                                                                                 gAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcgc        >  1:5178231/1‑93 (MQ=255)
                                                                                  aGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGcgcg       >  1:1813234/1‑93 (MQ=255)
                                                                                   ggggggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt      >  1:4456419/1‑93 (MQ=255)
                                                                                   ggggggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt      <  2:6721595/93‑1 (MQ=255)
                                                                                   ggggggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt      <  1:3447915/93‑1 (MQ=255)
                                                                                   ggggggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGt      <  1:5441599/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt     <  2:1227915/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt     <  2:4939506/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt     <  1:4659866/93‑1 (MQ=255)
                                                                                     ggggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTg    >  2:1936162/1‑93 (MQ=255)
                                                                                      ggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa   <  2:2608939/93‑1 (MQ=255)
                                                                                      ggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa   >  2:5310272/1‑93 (MQ=255)
                                                                                       ggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt  <  1:5855988/93‑1 (MQ=255)
                                                                                       ggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt  <  2:2309262/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAT  >  NC_000913/631151‑631330

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