Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,540,239 T→A 26.7% L26Q (CTG→CAG)  feoA → ferrous iron transporter, protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,540,2390TA26.7% 243.3 / 106.6 135L26Q (CTG→CAG) feoAferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base T (45/54);  new base A (18/18);  total (63/72)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 136 total aligned reads displayed.

TAATGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGG  >  NC_000913/3540140‑3540337
                                                                                                   |                                                                                                  
tAATGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACt                                                                                                    >  2:923622/1‑100 (MQ=255)
 aaTGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTg                                                                                                   >  2:1234977/1‑100 (MQ=255)
  aTGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGc                                                                                                  <  2:322340/100‑1 (MQ=255)
  aTGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGc                                                                                                  <  1:1327628/100‑1 (MQ=255)
  aTGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGc                                                                                                  >  1:218455/1‑100 (MQ=255)
    gagaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGTCATATCGCCAAAAACTGCtt                                                                                                >  2:24121/1‑100 (MQ=255)
    gagaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCGAAATACAGCtt                                                                                                <  2:1253056/100‑1 (MQ=255)
    gagaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCtt                                                                                                >  2:525236/1‑100 (MQ=255)
    gagaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCtt                                                                                                <  1:1075508/100‑1 (MQ=255)
    gagaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCtt                                                                                                >  1:2120231/1‑100 (MQ=255)
    gagaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATAAGCCCGGCATATGGCAAAAAACTGCtt                                                                                                >  1:462742/1‑100 (MQ=255)
      gaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTtc                                                                                              >  1:898669/1‑100 (MQ=255)
      gaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTtc                                                                                              >  2:1761929/1‑100 (MQ=255)
      gaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTtc                                                                                              <  1:458726/100‑1 (MQ=255)
      gaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTtc                                                                                              >  1:2162627/1‑100 (MQ=255)
        aaCAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTtctc                                                                                            >  2:600273/1‑100 (MQ=255)
        aaCAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTtctc                                                                                            <  1:1035430/100‑1 (MQ=255)
        aaCAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTtctc                                                                                            >  1:384764/1‑100 (MQ=255)
         aCAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTctct                                                                                           <  1:970470/100‑1 (MQ=255)
          cAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAGATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCtt                                                                                          <  1:1688671/100‑1 (MQ=255)
          cAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCtt                                                                                          >  1:1201290/1‑100 (MQ=255)
           aaGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTg                                                                                         <  2:312934/100‑1 (MQ=255)
            aGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTgg                                                                                        >  2:2197179/1‑100 (MQ=255)
              tAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCa                                                                                      <  2:218452/100‑1 (MQ=255)
                  cACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGtt                                                                                  <  1:2365893/100‑1 (MQ=255)
                  cACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGtt                                                                                  >  1:374789/1‑100 (MQ=255)
                  cACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGtt                                                                                  >  1:2166589/1‑100 (MQ=255)
                   aCCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTa                                                                                 <  2:1375742/100‑1 (MQ=255)
                      tATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAATCTGCTTTCTCTTGGCATGTTACCt                                                                              <  2:1273723/100‑1 (MQ=255)
                      tATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCt                                                                              >  1:1083694/1‑100 (MQ=255)
                          cAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGct                                                                          >  1:1577613/1‑100 (MQ=255)
                          cAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGct                                                                          <  2:816669/100‑1 (MQ=255)
                          cAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGct                                                                          <  2:1009299/100‑1 (MQ=255)
                           aaTACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGctc                                                                         >  1:1080634/1‑100 (MQ=255)
                           aaTACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGctc                                                                         <  1:945067/100‑1 (MQ=255)
                           aaTACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGct                                                                          >  2:309970/1‑99 (MQ=255)
                           aaTACACTCCAGATACTGCGTGGAAAATCACTGGCATTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGctc                                                                         >  2:726271/1‑100 (MQ=255)
                            aTACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGctcc                                                                        <  1:463194/100‑1 (MQ=255)
                              acacTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGctcctc                                                                      <  2:1860853/100‑1 (MQ=255)
                               cactccttgtACTGCTTGGAAAATCTCGGGTGTTTCCCGGTAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGctcct                                                                       <  1:1033255/89‑1 (MQ=255)
                               cacTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCctct                                                                     >  2:1613895/1‑100 (MQ=255)
                                 cTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCttt                                                                   <  1:860568/100‑1 (MQ=255)
                                  tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCtttt                                                                  <  1:2014954/100‑1 (MQ=255)
                                  tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCtt                                                                    <  1:28927/98‑1 (MQ=255)
                                  tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTAGCTCCTCtttt                                                                  >  1:306655/1‑100 (MQ=255)
                                  tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCtttt                                                                  <  1:893368/100‑1 (MQ=255)
                                   ccAGATACTGCGTGGAAAATTACTGTCTTTTCCCGTGAAATCAGCCAGAAAAAACAAAAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCttt                                                                   <  2:1516566/98‑1 (MQ=255)
                                    cAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTaa                                                                >  1:1176129/1‑100 (MQ=255)
                                    cAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTaa                                                                <  1:15232/100‑1 (MQ=255)
                                    cAGATACTGCGTGGAAAATCACTGGCTTTCCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTaa                                                                <  2:1439918/100‑1 (MQ=255)
                                     aGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAt                                                               <  1:734572/100‑1 (MQ=255)
                                      gATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAtg                                                              <  2:2306057/100‑1 (MQ=255)
                                      gATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAtg                                                              <  2:258843/100‑1 (MQ=255)
                                      gATACTGCGTCGAAAATCACTGTCTTTTCCCGTGAATTCAGAAAAAAAAAAAAAAAAAAACTGCTTGCTTTTGGCATGTTACCTGGCTCCTCTTTTAAtg                                                              <  2:1576255/100‑1 (MQ=255)
                                       aTACTGCGTGAAAAATCACTAGCTTTTCCCGTTAAATCAGTTCGAAAAAAAAAAAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAtg                                                              <  2:2277956/99‑1 (MQ=255)
                                        tACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAtgt                                                             <  2:1804290/99‑1 (MQ=255)
                                         aCTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATgtgg                                                           <  2:1603614/100‑1 (MQ=255)
                                            gCGTGGAGAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTgc                                                        >  1:1254391/1‑100 (MQ=255)
                                            gCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTgc                                                        <  1:1603828/100‑1 (MQ=255)
                                             cGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTgcg                                                       >  2:2187803/1‑100 (MQ=255)
                                             cGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTgcg                                                       <  1:730963/100‑1 (MQ=255)
                                               tGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGcgcg                                                     <  1:1436222/100‑1 (MQ=255)
                                                ggAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGt                                                    >  1:630344/1‑100 (MQ=255)
                                                 gAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTc                                                   <  2:1289452/100‑1 (MQ=255)
                                                 gAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTc                                                   <  1:1086047/100‑1 (MQ=255)
                                                 gAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCg                                                  <  2:210850/100‑1 (MQ=255)
                                                 gAAAATCACTGGCTTTTCCCATGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTAACTAGCTCATCTTTTAATGTGGTGCGCGTc                                                   >  2:250576/1‑100 (MQ=255)
                                                     aTCACTGGCTTTTCCCGTGAAGTCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCACGTCGCTc                                               <  1:509471/100‑1 (MQ=255)
                                                      tCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTcc                                              >  2:566341/1‑100 (MQ=255)
                                                       cACTGGCTTTTCCCGTGAAATCAGCCCGGCATGTCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTAGCTACa                                             >  2:727545/1‑100 (MQ=255)
                                                       cACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGAgcgacgcgcacca                                             <  1:358704/100‑11 (MQ=255)
                                                        aCTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCggttca                                             >  2:1793926/1‑94 (MQ=255)
                                                        aCTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCAc                                            <  1:848737/100‑1 (MQ=255)
                                                           ggCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCg                                         >  1:2152239/1‑100 (MQ=255)
                                                            gCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCgg                                        <  1:603135/100‑1 (MQ=255)
                                                                 tCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAcc                                   <  2:627101/100‑1 (MQ=255)
                                                                 ncccGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAcc                                   >  2:1415484/2‑100 (MQ=255)
                                                                  cccGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAccc                                  <  2:1517222/100‑1 (MQ=255)
                                                                  cccGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAcc                                   >  1:1146474/1‑99 (MQ=255)
                                                                   ccGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAccc                                  >  1:596222/1‑99 (MQ=255)
                                                                    ccagtaaTCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCa                                <  2:1242982/95‑1 (MQ=255)
                                                                     gTGAAATCAGCCCGGCATCTCGCCAAAAATTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCAt                               <  2:1034080/100‑1 (MQ=255)
                                                                     gTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCAt                               >  1:1461327/1‑100 (MQ=255)
                                                                      tGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCAtt                              >  2:1688947/1‑100 (MQ=255)
                                                                       gAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTc                             <  1:2146565/100‑1 (MQ=255)
                                                                         aaTCATCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCa                            <  1:2271225/99‑1 (MQ=255)
                                                                         aaTCAGCCCGGCCTATCACCAAAAACTGCTTTCTCTTGGCATGTTACCCGGCTCCTCTTTTAATGTGTTACGCATCGCTCCACCCGGCGACCCCATTCa                            >  2:1903381/1‑99 (MQ=255)
                                                                         aaTCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAGTGTGGAGCGCGTCGATCCACTCGGCGACCCCATTCat                           >  2:1619578/1‑100 (MQ=255)
                                                                          aTCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCata                          >  1:1439711/1‑100 (MQ=255)
                                                                            cAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATc                        >  1:1963513/1‑100 (MQ=255)
                                                                              gCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGa                      <  1:1790928/100‑1 (MQ=255)
                                                                                ccGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGaaa                    >  1:973160/1‑100 (MQ=255)
                                                                                 cGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAAc                   <  1:923633/100‑1 (MQ=255)
                                                                                 cGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTGACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAAc                   >  1:1695108/1‑100 (MQ=255)
                                                                                      tatCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCcgtc              >  2:193355/1‑100 (MQ=255)
                                                                                      tatCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCcgtc              >  1:772806/1‑100 (MQ=255)
                                                                                         cGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCgtg           <  1:1280369/100‑1 (MQ=255)
                                                                                             aaaaaCTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAg       >  2:1062566/1‑100 (MQ=255)
                                                                                             aaaaaCAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGCCgtgt          >  1:1283531/1‑97 (MQ=255)
                                                                                                  cTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTgg  <  1:1422313/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TAATGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGG  >  NC_000913/3540140‑3540337

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: