Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,697 A→C 5.8% 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 5.8% 237.5 / 3.6 82intergenic (+24/+30)ileX/nfeFtRNA‑Ile/NADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (44/33);  new base C (3/2);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTAA  >  NC_000913/3215549‑3215844
                                                                                                                                                    |                                                                                                                                                   
cgGTAATGGCTGGATTGCGACACGGAGTTACTTTATAATACGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAAt                                                                                                                                                    >  1:1309309/1‑150 (MQ=255)
      tGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAg                                                                                                                                              <  2:1352419/150‑1 (MQ=255)
       ggCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGGCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGc                                                                                                                                             <  2:163822/150‑1 (MQ=25)
       ggCGGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGTCTACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGc                                                                                                                                             <  1:795232/150‑1 (MQ=255)
         cTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGCct                                                                                                                                           <  1:1136514/150‑1 (MQ=255)
              ttGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttc                                                                                                                                      <  1:251422/150‑1 (MQ=255)
               tGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttct                                                                                                                                     <  2:991967/150‑1 (MQ=255)
                gCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttctt                                                                                                                                    <  1:1573430/150‑1 (MQ=255)
                 cGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttt                                                                                                                                   <  1:843461/150‑1 (MQ=255)
                  gACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttctttt                                                                                                                                  <  2:1244354/150‑1 (MQ=255)
                  gACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttctttt                                                                                                                                  <  2:597737/150‑1 (MQ=255)
                            tACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAg                                                                                                                        <  1:843912/150‑1 (MQ=255)
                            tACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAg                                                                                                                        <  1:810623/150‑1 (MQ=255)
                                ttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTAc                                                                                                                    >  1:1053552/1‑150 (MQ=255)
                                   taATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACttt                                                                                                                 >  2:855373/1‑150 (MQ=255)
                                      tCCGCTACCATGGCCCCTTAGCTCAGTGTTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCATATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCg                                                                                                              <  2:284688/150‑1 (MQ=255)
                                       ccGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGt                                                                                                             >  1:713592/1‑150 (MQ=255)
                                        cGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGTTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTg                                                                                                            <  1:505915/150‑1 (MQ=255)
                                              cATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTa                                                                                                      >  1:141463/1‑150 (MQ=255)
                                                   ccccTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAACTCCAGCAGGGGCCACCAGATATAGCAAAGGCTTACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCg                                                                                                 >  2:119648/1‑150 (MQ=255)
                                                     ccTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCt                                                                                               >  2:771451/1‑150 (MQ=255)
                                                         aGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCAc                                                                                           >  1:964174/1‑150 (MQ=255)
                                                          gCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACg                                                                                          <  2:1168741/150‑1 (MQ=255)
                                                           cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                          >  2:564182/1‑133 (MQ=255)
                                                           cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                          <  1:564182/133‑1 (MQ=255)
                                                           cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGt                                                                                         >  2:1367591/1‑150 (MQ=255)
                                                           cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGt                                                                                         >  2:1057128/1‑150 (MQ=255)
                                                           cTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGACAGTAAGCCGCTGCACGt                                                                                         >  1:1607274/1‑150 (MQ=255)
                                                             cAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAc                                                                                       >  2:27681/1‑150 (MQ=255)
                                                              aGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAcc                                                                                      <  2:367602/150‑1 (MQ=255)
                                                              aGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAcc                                                                                      >  2:519534/1‑150 (MQ=255)
                                                                 ggTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGc                                                                                   >  1:474593/1‑150 (MQ=255)
                                                                      gagCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTgg                                                                              >  2:370796/1‑150 (MQ=255)
                                                                      gagCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTgg                                                                              <  1:150936/150‑1 (MQ=255)
                                                                         cAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTc                                                                           <  2:1536794/150‑1 (MQ=255)
                                                                           ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAt                                                                         >  2:431705/1‑150 (MQ=255)
                                                                             cGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATAc                                                                       <  2:796391/150‑1 (MQ=255)
                                                                             cGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATAc                                                                       <  2:1309309/150‑1 (MQ=255)
                                                                                  cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  >  1:1474195/1‑150 (MQ=255)
                                                                                  cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  >  2:1327822/1‑150 (MQ=255)
                                                                                  cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  >  2:542953/1‑150 (MQ=255)
                                                                                  cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  <  2:581645/150‑1 (MQ=255)
                                                                                   aTAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTcc                                                                 <  1:1286234/150‑1 (MQ=255)
                                                                                     aaTCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGc                                                               >  2:645616/1‑150 (MQ=255)
                                                                                                cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGCCGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcggc                                                    >  1:1136881/1‑150 (MQ=255)
                                                                                                  cTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTTGAAgcggaac                                                  <  2:713592/150‑4 (MQ=255)
                                                                                                      ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTa                                              <  2:379601/150‑1 (MQ=255)
                                                                                                       tCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTaa                                             >  2:1118568/1‑150 (MQ=255)
                                                                                                        cAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTaaa                                            >  1:928628/1‑150 (MQ=255)
                                                                                                            tCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCt                                        >  1:1134358/1‑150 (MQ=255)
                                                                                                            tCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCt                                        >  1:316760/1‑150 (MQ=255)
                                                                                                             ccAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCtt                                       <  1:1057128/150‑1 (MQ=255)
                                                                                                              cagcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTAAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAAACTGTTACGCTTCAAAGCGGCGGCTTAAATTCTTa                                      >  2:462752/1‑150 (MQ=255)
                                                                                                                gcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAAc                                    >  2:918113/1‑150 (MQ=255)
                                                                                                                   ggggCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTGTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGAc                                 >  1:1592888/1‑150 (MQ=255)
                                                                                                                    gggCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACt                                >  1:1544346/1‑150 (MQ=255)
                                                                                                                       ccaccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTa                             <  1:771451/150‑1 (MQ=255)
                                                                                                                       ccaccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTa                             >  2:132850/1‑150 (MQ=255)
                                                                                                                        caccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTAc                            >  1:143487/1‑150 (MQ=255)
                                                                                                                        caccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCAGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTAc                            >  2:89971/1‑150 (MQ=255)
                                                                                                                         accaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTAcc                           <  2:141463/150‑1 (MQ=255)
                                                                                                                                ataGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCg                    <  2:1053552/150‑1 (MQ=255)
                                                                                                                                  aGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  >  1:890074/1‑150 (MQ=255)
                                                                                                                                   gCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGtt                 <  2:964174/150‑1 (MQ=255)
                                                                                                                                   gCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGtt                 >  1:656370/1‑150 (MQ=255)
                                                                                                                                   gCAAAGGCTGACGAGAAATCGTCAGCCATTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGtt                 >  1:569480/1‑150 (MQ=255)
                                                                                                                                    cAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTa                <  2:928628/150‑1 (MQ=255)
                                                                                                                                    aaaaGGCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCCGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTa                >  1:1162043/2‑150 (MQ=255)
                                                                                                                                      aaGGCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATc              >  2:1376514/1‑150 (MQ=255)
                                                                                                                                       aGGCTGACGATGTCTCGTCAGCCTGTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGGCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATcc             >  2:191767/1‑150 (MQ=255)
                                                                                                                                       aGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTAGTTTGCGTGCCAGTAAGCCGGTGCACGTTCCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTATATTCTTAACGACTTTACCTTCGCCGGTTATcc             >  2:1473461/1‑150 (MQ=255)
                                                                                                                                       aGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATcc             <  1:511985/150‑1 (MQ=255)
                                                                                                                                        ggCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCCGCGGCTTAAATTCGTAACGACTTTACCTTCGCCGGTTATCCa            >  1:1217561/1‑150 (MQ=255)
                                                                                                                                         gCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAg           >  1:1374977/1‑150 (MQ=255)
                                                                                                                                         gCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGGGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGACGGTTATCCAg           >  2:871819/1‑150 (MQ=255)
                                                                                                                                            gACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTa                                              <  1:1567852/112‑1 (MQ=255)
                                                                                                                                            gACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTa                                              >  2:1567852/1‑112 (MQ=255)
                                                                                                                                            gACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATg        >  1:1069599/1‑150 (MQ=255)
                                                                                                                                            gACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATg        <  2:1474195/150‑1 (MQ=255)
                                                                                                                                              cGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAGCGACTTTACCTTCGCCGGTTGTCCAGATGaa      >  2:1490867/1‑150 (MQ=255)
                                                                                                                                               gagaAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAg     <  2:143487/150‑1 (MQ=255)
                                                                                                                                                agaAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGGTATCCAGATGAAGt    >  2:823898/1‑150 (MQ=255)
                                                                                                                                                  aaaTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTaa  <  1:119648/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTAA  >  NC_000913/3215549‑3215844

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