Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,697 A→C 6.5% intergenic (+24/+30) ileX → / ← nfeF tRNA‑Ile/NADPH‑dependent ferric chelate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,215,6970AC 6.5% 320.8 / 13.7 108intergenic (+24/+30)ileX/nfeFtRNA‑Ile/NADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (51/50);  new base C (3/4);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

Only 100 of 110 total aligned reads displayed.

GTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTAA  >  NC_000913/3215551‑3215844
                                                                                                                                                  |                                                                                                                                                   
gTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCCGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCg                                                                                                                                                  >  2:1216543/1‑150 (MQ=255)
   aTGGCTGGATTGCGTCACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCa                                                                                                                                               <  1:656270/150‑1 (MQ=255)
     ggCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGc                                                                                                                                             <  2:628688/150‑1 (MQ=255)
      gCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGcc                                                                                                                                            <  1:728682/150‑1 (MQ=255)
      gCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGcc                                                                                                                                            <  1:878190/150‑1 (MQ=255)
      gCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGcc                                                                                                                                            <  2:1553459/150‑1 (MQ=255)
       cTGGATGGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCct                                                                                                                                           <  2:198016/150‑1 (MQ=255)
          gATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCAAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttt                                                                                                                                        >  2:1050754/1‑150 (MQ=255)
            ttGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCCCCAGATATAGCAAAGGCTGACGAGAAATCGGCAGCctttttg                                                                                                                                      >  2:1284513/1‑149 (MQ=255)
            ttGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGGCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttc                                                                                                                                      <  2:750822/150‑1 (MQ=255)
               ggACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTGCAAGTGCAGCAGGGGCGACCAGATATAGCAAAGGCTGACGAGAAATCGGCAGCctttttcttt                                                                                                                                   <  2:612725/149‑1 (MQ=255)
               cGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttt                                                                                                                                   >  1:1508462/1‑150 (MQ=255)
                gACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttctttt                                                                                                                                  >  1:806592/1‑150 (MQ=255)
                 acacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttttt                                                                                                                                 <  2:105138/150‑1 (MQ=255)
                 acacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACGCATAATCGCTTGGTCGCTGGGTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttttt                                                                                                                                 <  2:111204/150‑1 (MQ=255)
                  cacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                                                                                                                                >  1:983936/1‑150 (MQ=255)
                  cacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                                                                                                                                <  1:1574570/150‑1 (MQ=255)
                    cGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCGTTTTTCTTTTtata                                                                                                                              <  2:1379429/150‑1 (MQ=255)
                    cGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTtata                                                                                                                              <  1:1034314/150‑1 (MQ=255)
                            cTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGtt                                                                                                                      <  2:1395167/150‑1 (MQ=255)
                             tttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTa                                                                                                                     >  2:80862/1‑150 (MQ=255)
                                   aTCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGc                                                                                                               >  1:1003470/1‑150 (MQ=255)
                                    tCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCg                                                                                                              <  2:611762/150‑1 (MQ=255)
                                      cGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTg                                                                                                            <  1:214100/150‑1 (MQ=255)
                                       gCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGc                                                                                                           <  1:1044859/150‑1 (MQ=255)
                                          aCCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAg                                                                                                        <  2:1001136/150‑1 (MQ=255)
                                                gCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGcc                                                                                                  <  2:1143525/150‑1 (MQ=255)
                                                     ttAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGc                                                                                             <  2:983936/150‑1 (MQ=255)
                                                     ttAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGc                                                                                             <  1:509825/150‑1 (MQ=255)
                                                       aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCAc                                                                                           >  2:1209929/1‑150 (MQ=255)
                                                         cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGt                                                                                         >  1:802550/1‑150 (MQ=255)
                                                         cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGt                                                                                         >  2:1429800/1‑150 (MQ=255)
                                                          tCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTa                                                                                        >  2:617637/1‑150 (MQ=255)
                                                           cAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAc                                                                                       >  1:1472246/1‑150 (MQ=255)
                                                            aGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAcc                                                                                      <  1:80862/150‑1 (MQ=255)
                                                              gggTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCg                                                                                    <  2:1492953/149‑1 (MQ=255)
                                                                gTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCt                                                                                  >  1:565334/1‑150 (MQ=255)
                                                                  tAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCtgt                                                                                >  1:767207/1‑150 (MQ=255)
                                                                    gagCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACg                                                                                          <  2:1368606/138‑1 (MQ=255)
                                                                    gagCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACg                                                                                          >  1:1368606/1‑138 (MQ=255)
                                                                      ggaGGCGACTCATTTTCGCTTGGTGGCTGGTTCAAGTGCAGCAGGGGCCACCAGATATAGCTTTGGCTGACGTGAAATCGTCAGGCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCTGTAAGCGGCTGCACGTACGCGGTGTGGGt                                                                            <  2:139606/148‑1 (MQ=255)
                                                                      gCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGt                                                                            >  1:62445/1‑150 (MQ=255)
                                                                      gCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGt                                                                            >  2:1766/1‑150 (MQ=255)
                                                                      gCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGt                                                                            >  2:561870/1‑150 (MQ=255)
                                                                        aGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCa                                                                          <  2:565334/150‑1 (MQ=255)
                                                                         ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAt                                                                         >  1:1138592/1‑150 (MQ=255)
                                                                         ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAt                                                                         >  2:1243096/1‑150 (MQ=255)
                                                                            gACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACt                                                                      <  2:802550/150‑1 (MQ=255)
                                                                              cTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGt                                                                    >  2:1248699/1‑150 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCtgtg                                                                               <  2:652295/137‑1 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  >  1:460052/1‑150 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  <  2:1020332/150‑1 (MQ=255)
                                                                                 aTAATCGCTTGGTCGCGGGTTCAAGTCCAGCAGTGGCCACCAGATATAGCAAAGGCGGACGAGAAATCGACAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGGCATACTGTTcc                                                                 <  1:1328773/150‑1 (MQ=255)
                                                                                   aaTCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGc                                                               >  1:238358/1‑150 (MQ=255)
                                                                                      cGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTc                                                            <  1:561870/150‑1 (MQ=255)
                                                                                       gCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCa                                                           >  2:597939/1‑150 (MQ=255)
                                                                                        cTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCaa                                                          <  2:1019208/150‑1 (MQ=255)
                                                                                           ggTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGGCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgc                                                       >  1:699695/1‑150 (MQ=255)
                                                                                            gTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcg                                                      >  2:462911/1‑150 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                             >  2:98153/1‑125 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                             <  1:98153/125‑1 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcggc                                                    >  1:1166836/1‑150 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcggc                                                    >  2:1236027/1‑150 (MQ=255)
                                                                                                 tGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCggcggc                                                 <  2:1446586/150‑1 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                                <  2:140162/150‑1 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                                >  2:12737/1‑150 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                                <  1:1429800/150‑1 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                                >  2:971323/1‑150 (MQ=255)
                                                                                                    ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                          >  1:1570389/1‑90 (MQ=255)
                                                                                                    ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                          <  2:1570389/90‑1 (MQ=255)
                                                                                                         gTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTc                                         >  1:348032/1‑150 (MQ=255)
                                                                                                         gTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTc                                         >  2:206317/1‑150 (MQ=255)
                                                                                                           ccAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCtt                                       >  2:187630/1‑150 (MQ=255)
                                                                                                           ccAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCtt                                       >  1:1100275/1‑150 (MQ=255)
                                                                                                              gcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAAc                                    <  1:187630/150‑1 (MQ=255)
                                                                                                                 ggggCCACCAGATATAGCAAAGGCTGACGAGAGATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGAc                                 >  2:1559245/1‑150 (MQ=255)
                                                                                                                    gCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACttt                              >  2:66413/1‑150 (MQ=255)
                                                                                                                     ccaccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCATAAATTCTTAACGACTTTa                             <  1:617637/150‑1 (MQ=255)
                                                                                                                       accaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTAcc                           <  2:806592/150‑1 (MQ=255)
                                                                                                                       accaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTAcc                           <  2:1472246/150‑1 (MQ=255)
                                                                                                                        ccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCt                          >  2:1196707/1‑150 (MQ=255)
                                                                                                                          aGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTc                        >  2:205650/1‑150 (MQ=255)
                                                                                                                             tataGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGcc                     <  1:462911/150‑1 (MQ=255)
                                                                                                                                gaaaaaGGCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  2:583336/147‑1 (MQ=255)
                                                                                                                                aGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  2:460052/150‑1 (MQ=255)
                                                                                                                                aGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  1:971323/150‑1 (MQ=255)
                                                                                                                                aGCAAAGGCTGACGAGAAATCGGCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  1:936097/150‑1 (MQ=255)
                                                                                                                                 gCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACttt                              >  1:992571/1‑137 (MQ=255)
                                                                                                                                 gCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGtt                 >  2:917241/1‑150 (MQ=255)
                                                                                                                                 gCAAAGGCTGACGAGAAATCGACAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGGACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGCCTTAAATTCTTAACGACTATACCTTCGCCGGtt                 >  1:1139141/1‑150 (MQ=255)
                                                                                                                                   aaaGGCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTAt               >  1:409272/1‑150 (MQ=255)
                                                                                                                                     aGGCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATcc             >  1:1190924/1‑150 (MQ=255)
                                                                                                                                      ggCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCa            >  2:1032812/1‑150 (MQ=255)
                                                                                                                                       gCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  >  1:909248/1‑143 (MQ=255)
                                                                                                                                       gCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  2:909248/143‑1 (MQ=255)
                                                                                                                                        cTGACGAGAAATCGTCAGCATTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGa          <  1:1085708/150‑1 (MQ=255)
                                                                                                                                          gACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATg        <  1:239611/150‑1 (MQ=255)
                                                                                                                                            cGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGATTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCTCCGTTTATCCAGATGaa      >  1:373490/1‑150 (MQ=255)
                                                                                                                                               gaAATCGTCAGCCTTTTTCTTTTTATATATCGGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTa   >  2:81917/1‑150 (MQ=255)
                                                                                                                                                aaaTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCGGATGAAGTaa  <  2:1508462/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTAA  >  NC_000913/3215551‑3215844

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