Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,915 T→C 5.2% intergenic (‑190/+20) uof ← / ← fldA RyhB‑regulated fur leader peptide/flavodoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,9150TC 5.2% 243.9 / 2.4 77intergenic (‑190/+20)uof/fldARyhB‑regulated fur leader peptide/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base C (2/2);  total (37/40)
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.94e-01

AATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCATCT  >  NC_000913/710769‑711064
                                                                                                                                                  |                                                                                                                                                     
aaTATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                                                                                                                    <  2:118952/150‑1 (MQ=255)
   atTGTTAATAAAAACGTGGAACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCCGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCt                                                                                                                                                 >  2:1289978/1‑150 (MQ=255)
      gTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcg                                                                                                                                              <  2:1277192/150‑1 (MQ=255)
       ttAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGCCGTAGACATGAGTCTAcgc                                                                                                                                             >  1:830626/1‑150 (MQ=255)
        taataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                                                                                                                            >  1:627883/1‑150 (MQ=255)
        taataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                                                                                                                            >  1:632377/1‑150 (MQ=255)
        taataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                                                                                                                            >  1:752785/1‑150 (MQ=255)
          ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgccgc                                                                                                                                          >  2:401512/1‑150 (MQ=255)
           taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgccgta                                                                                                                                         <  1:1282743/150‑3 (MQ=255)
           taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                         >  2:1108957/1‑150 (MQ=255)
              aaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatca                                                                                                                                      <  1:1235726/150‑1 (MQ=255)
                aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGCCATTATCTATTAATACGGCGTAGACTCATGTCTACGCCGtattaat                                                                                                                                    <  1:530155/150‑5 (MQ=255)
                     gCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                         <  2:606516/140‑1 (MQ=255)
                     gCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                         >  1:606516/1‑140 (MQ=255)
                      cacaGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGc                                                                                                                              <  2:473613/150‑1 (MQ=255)
                          gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                                                                                                                          >  1:374823/1‑150 (MQ=255)
                            ttACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTgag                                                                                                                        <  2:890159/150‑1 (MQ=255)
                               cATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAt                                                                                                                     >  2:21281/1‑150 (MQ=255)
                                   tACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                                                                                                                 <  1:568849/150‑1 (MQ=255)
                                     cAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtc                                                                                                               >  2:407067/1‑150 (MQ=255)
                                      aaCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtcg                                                                                                              <  2:567844/150‑1 (MQ=255)
                                        cGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag                                                                                                            <  2:107324/150‑1 (MQ=255)
                                         ggCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgaga                                                                                                           >  1:754923/1‑150 (MQ=255)
                                            aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATCATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGc                                                                                                        <  2:97732/150‑1 (MQ=255)
                                             aGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCa                                                                                                       <  1:1253589/150‑1 (MQ=255)
                                              ggCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCaa                                                                                                      >  1:153651/1‑150 (MQ=255)
                                               gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                     <  1:129300/150‑1 (MQ=255)
                                               gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                     <  2:1046808/150‑1 (MQ=255)
                                               gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                     <  2:738179/150‑1 (MQ=255)
                                                 cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAActc                                                                                                   >  1:797297/1‑150 (MQ=255)
                                                 cACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACACGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAActc                                                                                                   >  2:121980/1‑150 (MQ=255)
                                                    ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc                                                                                                <  2:689523/150‑1 (MQ=255)
                                                    ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc                                                                                                <  2:694921/150‑1 (MQ=255)
                                                    ttGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTc                                                                                                <  2:700154/150‑1 (MQ=255)
                                                     tGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCa                                                                                               <  1:37025/150‑1 (MQ=255)
                                                       aaCACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGa                                                                                             <  2:771098/150‑1 (MQ=255)
                                                         cacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGaaa                                                                                           >  2:814624/1‑150 (MQ=255)
                                                             gggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTg                                                                                       >  1:259245/1‑150 (MQ=255)
                                                                 cATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTa                                                                                   <  2:754923/150‑1 (MQ=255)
                                                                   ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAAc                                                                                 <  2:1105508/150‑1 (MQ=255)
                                                                    tGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAAcc                                                                                <  2:423385/150‑1 (MQ=255)
                                                                         aGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCAttt                                                                           <  2:752785/150‑1 (MQ=255)
                                                                              caacaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTa                                                                      <  1:498032/150‑1 (MQ=255)
                                                                                acaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTaca                                                                    >  1:412426/1‑150 (MQ=255)
                                                                                 caGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTacac                                                                   >  2:145227/1‑150 (MQ=255)
                                                                                 caGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTaaaa                                                                   >  1:1299398/1‑147 (MQ=255)
                                                                                        aGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGc                                                            <  1:1289978/150‑1 (MQ=255)
                                                                                         gTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCg                                                           >  2:1280930/1‑150 (MQ=255)
                                                                                            aCAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTc                                                        <  1:1108957/150‑1 (MQ=255)
                                                                                            aCAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTc                                                        <  2:1394254/150‑1 (MQ=255)
                                                                                            aCAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTc                                                        <  2:627883/150‑1 (MQ=255)
                                                                                            aCAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTc                                                        <  2:632377/150‑1 (MQ=255)
                                                                                             cAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCa                                                       >  2:219535/1‑150 (MQ=255)
                                                                                              aaaTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAg                                                      >  1:170468/1‑150 (MQ=255)
                                                                                               aaTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCa                                                                                                       >  2:500383/1‑100 (MQ=255)
                                                                                               aaTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCa                                                                                                       <  1:500383/100‑1 (MQ=255)
                                                                                                 tttGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTc                                                   <  2:566746/150‑1 (MQ=255)
                                                                                                    gTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCgg                                                >  1:976726/1‑150 (MQ=255)
                                                                                                      aGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCt                                              <  2:412426/150‑1 (MQ=255)
                                                                                                       gCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTc                                                               >  2:296754/1‑132 (MQ=255)
                                                                                                       gCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTc                                                               <  1:296754/132‑1 (MQ=255)
                                                                                                       gCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTg                                             >  2:1166679/1‑150 (MQ=255)
                                                                                                                ttGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtc                                    >  1:1088976/1‑150 (MQ=255)
                                                                                                                ttGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtc                                    <  1:21281/150‑1 (MQ=255)
                                                                                                                     tGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGa                               >  2:1391338/1‑150 (MQ=255)
                                                                                                                         aTTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGc                           <  1:145227/150‑1 (MQ=255)
                                                                                                                               tATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGAcc                     >  2:395403/1‑150 (MQ=255)
                                                                                                                                  tAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGAc                  >  1:1015193/1‑150 (MQ=255)
                                                                                                                                   aaTACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACa                 >  2:69488/1‑150 (MQ=255)
                                                                                                                                   aaTACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACa                 >  2:272839/1‑150 (MQ=255)
                                                                                                                                      gcGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACaaat              >  2:1104558/2‑149 (MQ=255)
                                                                                                                                       cGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGt             <  2:1064952/150‑1 (MQ=255)
                                                                                                                                          cGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgt          >  2:1163992/1‑150 (MQ=255)
                                                                                                                                            tAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTACCCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgtcg        <  2:1000993/150‑1 (MQ=255)
                                                                                                                                            tAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgtcg        <  1:1420842/150‑1 (MQ=255)
                                                                                                                                              gACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgtcgtc      >  2:42750/1‑150 (MQ=255)
                                                                                                                                                cATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCAt    >  1:526671/1‑150 (MQ=255)
                                                                                                                                                  tGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCATCt  >  1:1431695/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCATCT  >  NC_000913/710769‑711064

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: