Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,661,793 C→A L113F (TTG→TTT)  iscR ← isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,661,7930CA91.8% 389.2 / 26.6 110L113F (TTG→TTT) iscRisc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base A (53/48);  total (58/52)
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.14e-02

Only 100 of 112 total aligned reads displayed.

AGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAA  >  NC_000913/2661646‑2661940
                                                                                                                                                   |                                                                                                                                                   
aGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTaaa                                                                                                                                                   >  2:1101281/1‑150 (MQ=255)
   gcgTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCa                                                                                                                                               <  3:1201885/151‑1 (MQ=255)
    cgTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCAc                                                                                                                                              >  3:279097/1‑151 (MQ=255)
     gTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCAc                                                                                                                                              >  3:944933/1‑150 (MQ=255)
     gTAACCTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACg                                                                                                                                             <  3:953667/151‑1 (MQ=255)
           tAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCcac                                                                                                                                         <  2:246135/149‑1 (MQ=255)
            aaCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgc                                                                                                                                      <  3:552209/151‑1 (MQ=255)
            aaCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgc                                                                                                                                      <  3:83048/151‑1 (MQ=255)
            aaCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgc                                                                                                                                      <  2:1421449/151‑1 (MQ=255)
            aaCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgc                                                                                                                                      >  3:62888/1‑151 (MQ=255)
            aaCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgc                                                                                                                                      >  3:530363/1‑151 (MQ=255)
             aCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGTGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgcg                                                                                                                                     >  3:1468903/1‑151 (MQ=255)
             aCGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAgcg                                                                                                                                     >  2:652335/1‑151 (MQ=255)
               gTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGcgcg                                                                                                                                   <  2:944146/151‑1 (MQ=255)
                tCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGt                                                                                                                                  >  2:305145/1‑151 (MQ=255)
                  gATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTgg                                                                                                                                >  3:1188532/1‑151 (MQ=255)
                      gcgTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCa                                                                                                                            <  3:1245543/151‑1 (MQ=255)
                        gTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAgg                                                                                                                          >  3:1629535/1‑151 (MQ=255)
                            tgtgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGc                                                                                                                         <  3:1112005/148‑1 (MQ=255)
                              tgtgCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCAttt                                                                                                                     >  2:817194/1‑150 (MQ=255)
                                 gCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATc                                                                                                                  <  2:86233/150‑1 (MQ=255)
                                      tGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCgccgcc                                                                                                            >  3:1211094/1‑151 (MQ=255)
                                        cGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCt                                                                                                          >  2:36295/1‑151 (MQ=255)
                                         gTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTg                                                                                                         >  2:1598093/1‑151 (MQ=255)
                                           ggCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGc                                                                                                       <  2:748028/151‑1 (MQ=255)
                                                    tGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCtt                                                                                              >  3:719747/1‑151 (MQ=255)
                                                     gAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCtt                                                                                              <  3:1147629/150‑1 (MQ=255)
                                                       gTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAc                                                                                           >  2:1189546/1‑151 (MQ=255)
                                                        tATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                          >  2:1332723/1‑151 (MQ=255)
                                                         aTGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTAcc                                                                                          <  2:1536305/150‑1 (MQ=255)
                                                              gacgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACa                                                                                    <  2:421541/151‑1 (MQ=255)
                                                               acgacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACaa                                                                                   <  3:1302859/151‑1 (MQ=255)
                                                                 gacCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACg                                                                                 >  2:407557/1‑151 (MQ=255)
                                                                          caTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATc                                                                        <  2:1349485/151‑1 (MQ=255)
                                                                          caTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCAt                                                                         <  2:1333110/150‑1 (MQ=255)
                                                                          caTCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATc                                                                        <  2:340021/151‑1 (MQ=255)
                                                                             ccAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTAc                                                                     <  2:884294/151‑1 (MQ=255)
                                                                              cAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTAc                                                                     >  3:25224/1‑150 (MQ=255)
                                                                              cAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACa                                                                    <  3:1458673/151‑1 (MQ=255)
                                                                               aGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTAc                                                                     >  2:196907/1‑149 (MQ=255)
                                                                               aGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAg                                                                   <  2:1412913/151‑1 (MQ=255)
                                                                                  aCTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAtt                                                                >  2:273248/1‑151 (MQ=255)
                                                                                  aCTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAt                                                                 <  3:767371/150‑1 (MQ=255)
                                                                                   cTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGAtt                                                                <  3:1358177/150‑1 (MQ=255)
                                                                                     tCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGt                                                             >  2:1085334/1‑151 (MQ=255)
                                                                                      ccTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTc                                                            >  3:971354/1‑151 (MQ=255)
                                                                                      ccTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTc                                                            >  3:489070/1‑151 (MQ=255)
                                                                                      ccTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTc                                                            <  3:703543/151‑1 (MQ=255)
                                                                                         ggTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCaa                                                          <  2:185822/150‑1 (MQ=255)
                                                                                         ggTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCa                                                           <  2:181748/149‑1 (MQ=255)
                                                                                          gTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACg                                                        <  3:806879/151‑1 (MQ=255)
                                                                                           ttattaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACgg                                                       >  2:1211396/1‑151 (MQ=255)
                                                                                            tattaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACgg                                                       <  2:771705/150‑1 (MQ=255)
                                                                                            tattaACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGgc                                                      >  2:914784/1‑151 (MQ=255)
                                                                                                  ccAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAt                                                >  2:1398872/1‑151 (MQ=255)
                                                                                                    aGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAAtt                                               >  2:946634/1‑150 (MQ=255)
                                                                                                     gTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAc                                             >  2:97209/1‑151 (MQ=255)
                                                                                                     gTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAc                                             >  2:311409/1‑151 (MQ=255)
                                                                                                     gTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTAc                                             <  3:135172/151‑1 (MQ=255)
                                                                                                       tCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACt                                            <  2:786686/150‑1 (MQ=255)
                                                                                                        cGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTc                                          >  3:542390/1‑151 (MQ=255)
                                                                                                        cGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTc                                          >  3:307610/1‑151 (MQ=255)
                                                                                                        cGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGCCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTc                                          <  3:87346/151‑1 (MQ=255)
                                                                                                         gCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCg                                         <  3:292024/151‑1 (MQ=255)
                                                                                                         gCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCg                                         >  3:333304/1‑151 (MQ=255)
                                                                                                          ccTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGc                                        >  2:42856/1‑151 (MQ=255)
                                                                                                           cTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGcc                                       >  2:528864/1‑151 (MQ=255)
                                                                                                            tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGcc                                       >  2:1307351/1‑150 (MQ=255)
                                                                                                            tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGcc                                       >  2:1324996/1‑150 (MQ=255)
                                                                                                            tAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGcc                                       >  3:902653/1‑150 (MQ=255)
                                                                                                             aaaGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCa                                      >  2:501700/1‑150 (MQ=255)
                                                                                                              aaGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAAc                                    >  2:674407/1‑151 (MQ=255)
                                                                                                               aGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACg                                   >  3:771927/1‑151 (MQ=255)
                                                                                                               aGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACg                                   <  2:1255114/151‑1 (MQ=255)
                                                                                                               aGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACg                                   <  2:1190285/151‑1 (MQ=255)
                                                                                                                gTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACgg                                  >  2:497814/1‑151 (MQ=255)
                                                                                                                  aaggttgttGAGAAACCCGGTGAGACGGTCGCTAAAATCGCGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATGTACAGATGCGTCAACGGCGCTAATTACTTCGCCAACGGc                                 <  3:602347/147‑1 (MQ=255)
                                                                                                                  aaTGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGc                                 <  3:970005/150‑1 (MQ=255)
                                                                                                                    tgttgtTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAt                              >  3:407017/1‑151 (MQ=255)
                                                                                                                    tgttgtTGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAt                              >  3:1187531/1‑151 (MQ=255)
                                                                                                                     gttgttGAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAt                              <  2:197386/150‑1 (MQ=255)
                                                                                                                        gttgAGAAAACCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGAtgctg                          <  2:368807/151‑1 (MQ=255)
                                                                                                                           gagaAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGctgctg                       <  2:1639180/151‑1 (MQ=255)
                                                                                                                               aaaCCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCAt                   >  3:492043/1‑151 (MQ=255)
                                                                                                                                aaCCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATc                  >  3:784255/1‑151 (MQ=255)
                                                                                                                                 aCCGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCt                 >  2:817723/1‑151 (MQ=255)
                                                                                                                                   cGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt               <  3:21139/151‑1 (MQ=255)
                                                                                                                                   cGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt               <  3:1128623/151‑1 (MQ=255)
                                                                                                                                   cGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCttt               <  3:808048/151‑1 (MQ=255)
                                                                                                                                   cGGTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCtt                >  2:61808/1‑150 (MQ=255)
                                                                                                                                     gTGAGACGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGc             >  3:910056/1‑151 (MQ=255)
                                                                                                                                       gagaCGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCt           <  3:334923/151‑1 (MQ=255)
                                                                                                                                       gagaCGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCt           <  2:1390051/151‑1 (MQ=255)
                                                                                                                                       gagaCGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCt           >  3:597509/1‑151 (MQ=255)
                                                                                                                                        agaCGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTa          <  2:1044687/151‑1 (MQ=255)
                                                                                                                                         gaCGGTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTa          >  2:525692/1‑150 (MQ=255)
                                                                                                                                            ggTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAAc        >  3:894067/1‑149 (MQ=255)
                                                                                                                                            ggTCGCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAAc        >  2:346167/1‑149 (MQ=255)
                                                                                                                                                gCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATaa  <  2:622046/151‑1 (MQ=255)
                                                                                                                                                gCTAAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATaa  <  2:778592/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCGCGTAACTTAACGTCGATCGCGTCTTGTGTGCGGGTGCGTGGCGCGTCGTGAGTATGCTGACGACCAGACACATCCAGCACTTCCTGGTTATTAACCAGTTCGCCTAAAGTAATGTTGTTGAGAAAACCGGTGAGACGGTCGCTCAAATCACGCCACAGCGCGTGGGTCAGGCATTTATCGCCGCCCTGGCAGCCGCCTTTACCCTGACAACGGGTGGCATCTACAGATTCGTCAACGGCGCTAATTACTTCGCCAACGGCGATGCTGCTGGCATCTTTGCCTAACAGATAA  >  NC_000913/2661646‑2661940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑