Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,694 G→T 6.5% intergenic (+21/+33) 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,6940GT 6.5% 361.8 / 8.1 107intergenic (+21/+33)ileX/nfeFtRNA‑Ile/NADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base T (3/4);  total (53/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 = 9.79e-01

Only 100 of 108 total aligned reads displayed.

GTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTA  >  NC_000913/3215551‑3215843
                                                                                                                                               |                                                                                                                                                     
gTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCCGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCg                                                                                                                                                 >  2:1216543/1‑150 (MQ=255)
   aTGGCTGGATTGCGTCACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCa                                                                                                                                              <  1:656270/150‑1 (MQ=255)
      gCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGcc                                                                                                                                           <  1:728682/150‑1 (MQ=255)
      gCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGcc                                                                                                                                           <  2:1553459/150‑1 (MQ=255)
       cTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGCCTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGcat                                                                                                                                          >  2:1085708/1‑148 (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)
                 acacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACGCATAATCGCTTGGTCGCTGGGTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttttt                                                                                                                                <  2:111204/150‑1 (MQ=255)
                  cacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                                                                                                                               >  1:983936/1‑150 (MQ=255)
                  cacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                                                                                                                               <  1:1574570/150‑1 (MQ=255)
                    cGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTtata                                                                                                                             <  1:1034314/150‑1 (MQ=255)
                          tACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAg                                                                                                                       <  1:950293/150‑1 (MQ=255)
                            cTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGtt                                                                                                                     <  2:1395167/150‑1 (MQ=255)
                             tttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTa                                                                                                                    >  2:80862/1‑150 (MQ=255)
                             tttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTa                                                                                                                    <  2:1363440/150‑1 (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:1546731/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                                                                                        >  2:1429800/1‑150 (MQ=255)
                                                         cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGt                                                                                        >  1:802550/1‑150 (MQ=255)
                                                          tCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTa                                                                                       >  2:617637/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                                                                                         >  1:1368606/1‑138 (MQ=255)
                                                                    gagCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACg                                                                                         <  2:1368606/138‑1 (MQ=255)
                                                                      ggaGGCGACTCATTTTCGCTTGGTGGCTGGTTCAAGTGCAGCAGGGGCCACCAGATATAGCTTTGGCTGACGTGAAATCGTCAGGCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCTGTAAGCGGCTGCACGTACGCGGTGTGGGt                                                                           <  2:139606/148‑1 (MQ=255)
                                                                      gCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGt                                                                           >  1:62445/1‑150 (MQ=255)
                                                                      gCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGt                                                                           >  2:561870/1‑150 (MQ=255)
                                                                      gCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGt                                                                           >  2:1766/1‑150 (MQ=255)
                                                                        aGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCa                                                                         <  2:565334/150‑1 (MQ=255)
                                                                         ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAt                                                                        >  2:1243096/1‑150 (MQ=255)
                                                                         ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAt                                                                        >  1:1138592/1‑150 (MQ=255)
                                                                            gACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACt                                                                     <  2:802550/150‑1 (MQ=255)
                                                                              cTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGt                                                                   >  2:1248699/1‑150 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCtgtg                                                                              >  1:652295/1‑137 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCtgtg                                                                              <  2:652295/137‑1 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                 <  2:1020332/150‑1 (MQ=255)
                                                                                cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                 >  1:460052/1‑150 (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)
                                                                                             tCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcgg                                                    >  1:1357129/1‑150 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                            <  1:98153/125‑1 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                            >  2:98153/1‑125 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcggc                                                   >  2:1236027/1‑150 (MQ=255)
                                                                                              cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcggc                                                   >  1:1166836/1‑150 (MQ=255)
                                                                                                 tGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCggcggc                                                <  2:1446586/150‑1 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                               >  2:12737/1‑150 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                               <  1:1429800/150‑1 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                               <  2:140162/150‑1 (MQ=255)
                                                                                                  ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                               >  2:971323/1‑150 (MQ=255)
                                                                                                    ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                         <  2:1570389/90‑1 (MQ=255)
                                                                                                    ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                         >  1:1570389/1‑90 (MQ=255)
                                                                                                      cAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTaaa                                           <  1:241267/150‑1 (MQ=255)
                                                                                                         gTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTc                                        >  1:348032/1‑150 (MQ=255)
                                                                                                         gTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTc                                        >  2:206317/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:1472246/150‑1 (MQ=255)
                                                                                                                       accaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTAcc                          <  2:806592/150‑1 (MQ=255)
                                                                                                                        ccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCt                         >  2:1196707/1‑150 (MQ=255)
                                                                                                                          aGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTc                       >  2:205650/1‑150 (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)
                                                                                                                                      ggCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCa           >  2:487223/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)
                                                                                                                                               |                                                                                                                                                     
GTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTA  >  NC_000913/3215551‑3215843

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: