Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,541,791 G→A 100% G462D (GGT→GAT)  feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,541,7910GA98.2% 380.3 / ‑2.2 111G462D (GGT→GAT) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (48/61);  total (50/61)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 112 total aligned reads displayed.

TGTCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGGTATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTATCAT  >  NC_000913/3541695‑3541885
                                                                                                |                                                                                              
tGTCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGatat                                                                                             <  2:2827466/100‑1 (MQ=255)
 gTCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATAtt                                                                                            >  1:385472/1‑100 (MQ=255)
  tCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATtg                                                                                           >  2:2619864/1‑100 (MQ=255)
  tCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATtg                                                                                           <  1:3153971/100‑1 (MQ=255)
  tCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATtg                                                                                           >  2:2818584/1‑100 (MQ=255)
   ccTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATtgt                                                                                          <  1:2022888/100‑1 (MQ=255)
    cTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATtgtg                                                                                         >  1:906634/1‑100 (MQ=255)
    cTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATtgtg                                                                                         <  1:1417627/100‑1 (MQ=255)
     tGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGa                                                                                        <  2:1646355/100‑1 (MQ=255)
      gcggcgCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGa                                                                                        <  1:1482760/99‑1 (MQ=255)
       cggcgCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGTCGGTCTTCTCTCTGTACATGCTGGATATTGTGAt                                                                                       >  1:2320245/1‑99 (MQ=255)
        ggcgCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATgg                                                                                     >  1:2175695/1‑100 (MQ=255)
          cgcgcgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCg                                                                                   >  2:742818/1‑100 (MQ=255)
          cgcgcgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCg                                                                                   >  2:360940/1‑100 (MQ=255)
          cgcgcgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCg                                                                                   >  2:2746411/1‑100 (MQ=255)
           gcgcgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCgg                                                                                  <  1:3053601/100‑1 (MQ=255)
            cgcgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGt                                                                                 <  1:571684/100‑1 (MQ=255)
            cgcgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGt                                                                                 <  2:2815575/100‑1 (MQ=255)
              cgTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGc                                                                               <  2:2995222/100‑1 (MQ=255)
               gTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCt                                                                              <  2:942040/100‑1 (MQ=255)
                   ggCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACt                                                                          >  1:2278905/1‑100 (MQ=255)
                   ggCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACt                                                                          <  1:3088622/100‑1 (MQ=255)
                    gatATCTTCGCAGGATTCGCGGCTGCCTTCTACGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTg                                                                         >  2:2255877/3‑100 (MQ=255)
                     cTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTgg                                                                        <  2:1509573/100‑1 (MQ=255)
                     cTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTgg                                                                        <  1:3154634/100‑1 (MQ=255)
                      tATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGc                                                                       >  2:906557/1‑100 (MQ=255)
                       aTCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGcc                                                                      <  1:766438/100‑1 (MQ=255)
                       aTCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGcc                                                                      <  1:716332/100‑1 (MQ=255)
                       aTCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGcc                                                                      <  1:1542908/100‑1 (MQ=255)
                       aTCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGcc                                                                      <  1:1556217/100‑1 (MQ=255)
                        tCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCt                                                                     <  2:2220681/100‑1 (MQ=255)
                         cTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTg                                                                    >  2:945557/1‑100 (MQ=255)
                         cTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTg                                                                    >  2:718204/1‑100 (MQ=255)
                          ttCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGa                                                                   >  2:33986/1‑100 (MQ=255)
                           tCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGAt                                                                  <  2:45489/100‑1 (MQ=255)
                               cgTAGTCGCGGCTGCCTTCTGCGGGCAGAACGGTGCGCCGGCGGTCGTCTCGCTGTATATGCTGGATATTGTGAGGGCGGTGCTGACTGGCCTGATGCt                                                               <  1:2599634/98‑1 (MQ=255)
                               aGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTc                                                              >  1:3053839/1‑100 (MQ=255)
                                gTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCa                                                             <  2:906637/100‑1 (MQ=255)
                                 tATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCaa                                                            <  2:2649241/100‑1 (MQ=255)
                                  aTTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAg                                                           <  2:2945345/100‑1 (MQ=255)
                                   ttCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGt                                                          <  1:2168645/100‑1 (MQ=255)
                                   ttCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGt                                                          <  2:2486689/100‑1 (MQ=255)
                                    tCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTa                                                         <  1:3021372/100‑1 (MQ=255)
                                    tCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTa                                                         >  2:1936016/1‑100 (MQ=255)
                                    tCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTa                                                         <  1:288017/100‑1 (MQ=255)
                                     cgcgGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTac                                                        >  2:2495950/1‑100 (MQ=255)
                                      gcgGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTaca                                                       <  2:2175711/100‑1 (MQ=255)
                                          cTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACcat                                                   >  2:987955/1‑100 (MQ=255)
                                            gCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACcatca                                                 <  2:860936/100‑1 (MQ=255)
                                            gCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACcatca                                                 <  2:861163/100‑1 (MQ=255)
                                             ccTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACcatcat                                                >  2:2150875/1‑100 (MQ=255)
                                               ttcttcGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATgc                                              >  1:2479613/1‑100 (MQ=255)
                                                tcttcGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATgcg                                             >  1:1250722/1‑100 (MQ=255)
                                                tcttcGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATgcg                                             >  1:2079288/1‑100 (MQ=255)
                                                 cttcGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATgcgc                                            <  2:2432144/100‑1 (MQ=255)
                                                    cGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGTGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGt                                         <  2:2178267/100‑1 (MQ=255)
                                                     gggCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTg                                        >  1:2550841/1‑100 (MQ=255)
                                                     gggCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTACTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTg                                        <  1:2988785/100‑1 (MQ=255)
                                                      ggCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGa                                       <  2:2577986/100‑1 (MQ=255)
                                                      ggCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGa                                       >  1:1966979/1‑100 (MQ=255)
                                                      ggCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGa                                       >  2:1171037/1‑100 (MQ=255)
                                                        cAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAg                                     <  1:1991272/100‑1 (MQ=255)
                                                         aGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGc                                    >  1:2442157/1‑100 (MQ=255)
                                                           aaCGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGa                                  <  1:2969072/100‑1 (MQ=255)
                                                            aCGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGAc                                 <  1:2749411/100‑1 (MQ=255)
                                                                tGCGCTGGCGGTCTTCTCGCTGTATATTCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCg                             <  1:3297054/100‑1 (MQ=255)
                                                                tGCGCTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCg                             >  2:3267171/1‑100 (MQ=255)
                                                                 gcgcTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGt                            >  2:1396524/1‑100 (MQ=255)
                                                                   gcTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGttt                          >  2:1653626/1‑100 (MQ=255)
                                                                   gcTGGCGGTCTTCTCGCTGTATATGCTGGATATAGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGttt                          >  1:2014370/1‑100 (MQ=255)
                                                                    cTGGCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTg                         >  2:373961/1‑100 (MQ=255)
                                                                      ggCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTc                       >  1:1214993/1‑100 (MQ=255)
                                                                      ggCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTc                       >  2:441305/1‑100 (MQ=255)
                                                                       gCGGTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCGGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCa                      <  1:2120996/100‑1 (MQ=255)
                                                                         ggTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATg                    >  2:2708061/1‑100 (MQ=255)
                                                                         ggTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATg                    >  1:3257628/1‑100 (MQ=255)
                                                                          gTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATgg                   <  2:1749247/100‑1 (MQ=255)
                                                                          gTCTTCTCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATgg                   <  1:799813/100‑1 (MQ=255)
                                                                           tcttctCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATgg                   >  2:1887667/1‑99 (MQ=255)
                                                                           tcttctCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGa                  >  1:3268835/1‑100 (MQ=255)
                                                                           tcttctCGCTATATATGCTGGGTATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGa                  >  1:2272308/1‑100 (MQ=255)
                                                                            cttctCGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAg                 >  2:1709249/1‑100 (MQ=255)
                                                                              tctcGCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGtt               >  2:935376/1‑100 (MQ=255)
                                                                                  gCTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGATCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCg           <  2:3322398/100‑1 (MQ=255)
                                                                                   cTGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCgg          <  2:42137/100‑1 (MQ=255)
                                                                                    tGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCgg          <  1:1153250/99‑1 (MQ=255)
                                                                                    tGTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGt         <  2:1750853/100‑1 (MQ=255)
                                                                                     gTATATGCTGGGTATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTc        >  2:1168458/1‑100 (MQ=255)
                                                                                     gTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTc        <  1:373961/100‑1 (MQ=255)
                                                                                     gTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTc        >  1:672929/1‑100 (MQ=255)
                                                                                     gTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTc        <  1:897287/100‑1 (MQ=255)
                                                                                     gTATATGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTc        <  1:1299004/100‑1 (MQ=255)
                                                                                      tataTGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCt       <  2:1483033/100‑1 (MQ=255)
                                                                                      tataTGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCt       <  2:1737172/100‑1 (MQ=255)
                                                                                       atatGCTGGATATTGTGATGGCGGTGCTGACTTGCCTGATTCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTa      <  1:2818560/100‑1 (MQ=255)
                                                                                       atatGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTa      <  2:597108/100‑1 (MQ=255)
                                                                                        tatGCTGGATATTGTGATGGCGGTGCTGCCTGGCCTGATGCTCAATTCCGCCATCATGCGCGGTTAAGCGACGCCGTTTGTCATGGCGTTGCCGGTCTAt     <  1:1505110/100‑1 (MQ=255)
                                                                                        tatGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTAt     >  2:2582577/1‑100 (MQ=255)
                                                                                         atGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTATc    <  2:1721597/100‑1 (MQ=255)
                                                                                          tGCTGGATATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTATCa   <  1:301906/100‑1 (MQ=255)
                                                                                                |                                                                                              
TGTCCTGCGGCGCGCGTCTGGCTATCTTCGCAGTATTCGCGGCTGCCTTCTTCGGGCAGAACGGTGCGCTGGCGGTCTTCTCGCTGTATATGCTGGGTATTGTGATGGCGGTGCTGACTGGCCTGATGCTCAAGTACACCATCATGCGCGGTGAAGCGACGCCGTTTGTCATGGAGTTGCCGGTCTATCAT  >  NC_000913/3541695‑3541885

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