Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,586,614 G→A 5.3% L50F (CTT→TTT)  hsdR ← endonuclease R Type I restriction enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,586,6140GA 5.3% 356.1 / 5.3 112L50F (CTT→TTT) hsdRendonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base G (52/53);  new base A (3/3);  total (55/57)
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

Only 100 of 116 total aligned reads displayed.

GTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTG  >  NC_000913/4586465‑4586763
                                                                                                                                                     |                                                                                                                                                     
gTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAg                                                                                                                                                       >  2:658356/1‑150 (MQ=255)
gTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAg                                                                                                                                                       >  2:940953/1‑150 (MQ=255)
 ttATGATATTCTTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATTCAGTGATGATGTTGTCATCAACACAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGa                                                                                                                                                      <  1:1171668/150‑1 (MQ=255)
       tatTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTc                                                                                                                                                >  2:959349/1‑150 (MQ=255)
       tatTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTc                                                                                                                                                >  2:786657/1‑150 (MQ=255)
        atTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCg                                                                                                                                               >  2:957256/1‑150 (MQ=255)
         ttCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                                                              <  1:137093/150‑1 (MQ=255)
             tGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCtgtg                                                                                                                                          <  2:860193/150‑1 (MQ=255)
              gCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTgg                                                                                                                                         <  1:391600/150‑1 (MQ=255)
               cACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGc                                                                                                                                        >  1:1214411/1‑150 (MQ=255)
                          tACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACa                                                                                                                             <  1:363031/150‑1 (MQ=255)
                           aCCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACAt                                                                                                                            >  1:963439/1‑150 (MQ=255)
                            ccACTGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGTATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATa                                                                                                                           <  1:181175/150‑1 (MQ=255)
                             cAATGGTGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATAc                                                                                                                          <  2:1219227/150‑1 (MQ=255)
                                 gcggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAt                                                                                                                      <  2:462986/150‑1 (MQ=255)
                                  cggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAtt                                                                                                                     >  1:1086333/1‑150 (MQ=255)
                                   ggcgTAATTTGTGAAATACAGAGAGGATGTGGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAttt                                                                                                                    <  2:67967/150‑1 (MQ=255)
                                     cgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTa                                                                                                                  <  2:502472/150‑1 (MQ=255)
                                      gTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTaa                                                                                                                 >  1:397626/1‑150 (MQ=255)
                                      gTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTaa                                                                                                                 >  2:1017954/1‑150 (MQ=255)
                                      gTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCAGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGACAAACATACGCATTTTaa                                                                                                                 >  1:202515/1‑150 (MQ=255)
                                        aaTTTGTGAAATACAGAGAGGATGTTGTCATTAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATc                                                                                                               <  2:895775/150‑1 (MQ=255)
                                        aaTTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATc                                                                                                               <  1:443315/150‑1 (MQ=255)
                                          tttGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAg                                                                                                             <  1:593744/150‑1 (MQ=255)
                                          tttGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAg                                                                                                             <  2:1117163/150‑1 (MQ=255)
                                              tgAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtc                                                                                                         <  1:617289/150‑1 (MQ=255)
                                               gAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGCCCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtcg                                                                                                        >  1:712095/1‑150 (MQ=255)
                                               gAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtcg                                                                                                        <  2:1152581/150‑1 (MQ=255)
                                                 aaTACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGCTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTTAGTAAAAGACCAAAATGTTTCGCTGTGGCTTCGCCAAACATAAGCATTTTAATCAGCGTCttg                                                                                                      >  2:361754/1‑150 (MQ=255)
                                                 aaTACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGCTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAAATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCgtg                                                                                                      >  2:149511/1‑150 (MQ=255)
                                                 aaTACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGCTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAAATGTTTCGCTGTGGCTTCGCCAAACATACGCAGTTTAATCAGCGTCgtg                                                                                                      >  2:366664/1‑150 (MQ=255)
                                                    aCAGAGAGGATGTTGTCTTCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTg                                                                                                   <  1:1246714/150‑1 (MQ=255)
                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGa                                                                                                <  1:1055714/150‑1 (MQ=255)
                                                        agagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGat                                                                                               <  2:736792/150‑1 (MQ=255)
                                                         gagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatc                                                                                              >  2:971934/1‑150 (MQ=255)
                                                           ggATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatcat                                                                                            <  1:439839/150‑1 (MQ=255)
                                                              tgttgtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCg                                                                                         <  1:298995/150‑1 (MQ=255)
                                                                    catcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGt                                                                                   >  1:579359/1‑150 (MQ=255)
                                                                     atcaACAAATGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAATGGGGCATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGGAGtt                                                                                  >  1:1256104/1‑150 (MQ=255)
                                                                     atcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGtt                                                                                  >  1:176526/1‑150 (MQ=255)
                                                                       caacaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAt                                                                                <  1:1022559/150‑1 (MQ=255)
                                                                       caacaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAt                                                                                >  1:1274648/1‑150 (MQ=255)
                                                                       aaaaaaaGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAt                                                                                <  1:1181788/146‑1 (MQ=255)
                                                                        aacaaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAtt                                                                               >  1:609026/1‑150 (MQ=255)
                                                                          caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAtttt                                                                             >  2:1010132/1‑150 (MQ=255)
                                                                            aacgCGATTTTGCCGAGTTCCCGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTCACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTcc                                                                           <  2:968318/147‑1 (MQ=255)
                                                                            aaGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTcc                                                                           >  2:1214575/1‑150 (MQ=255)
                                                                            aaGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTcc                                                                           <  1:1253499/150‑1 (MQ=255)
                                                                               gCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTccgcc                                                                        >  1:1169313/1‑150 (MQ=255)
                                                                                 gATTTTGCCGAGTTCAGGCAGGAGTTCGTGTTGATTCTCACTAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCcgccgc                                                                      <  2:1108958/150‑1 (MQ=255)
                                                                                    tttGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGcaca                                                                   >  1:1337083/1‑150 (MQ=255)
                                                                                       gCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGc                                                                <  1:1017954/150‑1 (MQ=255)
                                                                                       gCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGc                                                                <  2:1104916/150‑1 (MQ=255)
                                                                                       gCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGc                                                                >  1:1331013/1‑150 (MQ=255)
                                                                                       gCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGAACAGGc                                                                >  1:220421/1‑150 (MQ=255)
                                                                                           tgTTCACGCAGGAGATCGTGTTGATTCTCAAAAGGGGGGATGTGGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGAGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATg                                                            <  1:807131/149‑1 (MQ=255)
                                                                                           aGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATg                                                            <  1:1040257/150‑1 (MQ=255)
                                                                                              tCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCa                                                         >  1:620595/1‑150 (MQ=255)
                                                                                               cACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCAt                                                        >  1:675651/1‑150 (MQ=255)
                                                                                               cACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGAGTCGTGATGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCAt                                                        >  1:191922/1‑150 (MQ=255)
                                                                                                aCGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATa                                                       <  1:959349/150‑1 (MQ=255)
                                                                                                aCGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATa                                                       <  2:960573/150‑1 (MQ=255)
                                                                                                       agaTCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAg                                                >  2:480640/1‑150 (MQ=255)
                                                                                                        gaTCGTGTTGTTTCTCAAAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGATTCGCCAAACATACGCATTTTATTCAGCGTCGTGTGGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGt                                               <  2:597221/150‑1 (MQ=255)
                                                                                                        gaTCGTGTTGATTCTCACAAGGGTGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGt                                               <  1:896670/150‑1 (MQ=255)
                                                                                                        gaTCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGt                                               <  1:138451/150‑1 (MQ=255)
                                                                                                        gaTCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGt                                               >  1:604071/1‑150 (MQ=255)
                                                                                                        gaTCGTGTGGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGt                                               <  1:951342/150‑1 (MQ=255)
                                                                                                          tCGTGTTGACTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCg                                             <  2:314247/150‑1 (MQ=255)
                                                                                                                tGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATTGCATAAGTGAAGTCGTTGACg                                       >  2:1002854/1‑150 (MQ=255)
                                                                                                                 gaatctcACCAGGGGGGATGTTGACTAACAGACCAAGATGTTTCGCCGTGGCGTCGCCAAACATACGCATTTTAATCAGCCTCGTGTTGGTATCATCCGGGGAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGCTGACGc                                      <  2:6666/147‑1 (MQ=255)
                                                                                                                 gATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGc                                      >  1:277975/1‑150 (MQ=255)
                                                                                                                         caAGTGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAgg                              <  1:1214575/150‑1 (MQ=255)
                                                                                                                         caAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAgg                              <  1:221256/150‑1 (MQ=255)
                                                                                                                         caAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAgg                              >  2:572665/1‑150 (MQ=255)
                                                                                                                                ggATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTTGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATTGCATAAGTGAAGTCGTTGAAGCCATTAAGGAATACaa                       >  1:77340/1‑150 (MQ=255)
                                                                                                                                 gATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGGCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCaaa                      <  2:1214411/150‑1 (MQ=255)
                                                                                                                                  aTGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCaaaa                     >  2:536587/1‑150 (MQ=255)
                                                                                                                                  aTGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCCCACAGGCGATGGCAAAAGTGAAGTATTTGACGCCCTTCAGGAAttacaaa                     >  2:197098/1‑145 (MQ=255)
                                                                                                                                   tGTTGAGAAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAAt                    <  1:751324/150‑1 (MQ=255)
                                                                                                                                     ttGAGTAACAGACCAAGATGGTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTg                  <  1:94189/150‑1 (MQ=255)
                                                                                                                                      tGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATAATCCGTGTAGATATTTTCCGCAGCACCGGCGATGGCATAAGTGAAGACATTGACGCCCTTCAGGAATTCAAAATTgg                 >  2:731470/1‑150 (MQ=255)
                                                                                                                                      tGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATAAGCGTCGTGTTGGGAACATCCGGGTAGTTAATTTCCGCCGCACAGGCGATGGCATAAGTGAAGTAGTTGACGCCCTTCAGGAATTCAAAATTgg                 >  1:1144244/1‑150 (MQ=255)
                                                                                                                                        aGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAt               <  2:29491/150‑1 (MQ=255)
                                                                                                                                        aGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCAACAGGTTAGATATTTTCATCACCACAGTCGATGGAAAAATAGAATTAATTGACGCACTTCAGGAATTCAAAATTGGAt               >  2:785589/1‑150 (MQ=255)
                                                                                                                                         gTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCATGAATTCAAAATTGGAtt              >  1:1079741/1‑150 (MQ=255)
                                                                                                                                         gTAACAGACCAAGATGGTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGGGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAtt              <  2:792974/150‑1 (MQ=255)
                                                                                                                                          tAACAGACCAAGGTGATGCGCTGTGGCTTAGCCAAACATACGCATTTTAATCAGCGTCCTGTTGGGATCATCCGGGTAGTTATTTTCCTCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAttt             <  2:882845/150‑1 (MQ=255)
                                                                                                                                          tAACAGACCAAGATGTTTCGCTGTGGCTTCGACAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGAttt             >  1:71254/1‑150 (MQ=255)
                                                                                                                                             cAGACCAATATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTAtt          <  2:1226408/150‑1 (MQ=255)
                                                                                                                                              aGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTATAATAACCGTCGTGTTGGGATCATCAGTGTAGTTATTTTCCCCCCCCAAGGCGATGGCTAAAGTGAACTCATTGACGCCCTTCAGGAATTCAAAATTGGAATTAtaa         >  2:682924/1‑148 (MQ=255)
                                                                                                                                              aGACCAAAATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtc         <  1:149511/150‑1 (MQ=255)
                                                                                                                                              aGACCAAAATGTATCGCTGTGGCTTCGCCAAACATACGCATTTAAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtc         <  1:361754/150‑1 (MQ=255)
                                                                                                                                               gACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtca        >  2:599622/1‑150 (MQ=255)
                                                                                                                                               gACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtca        >  1:1086335/1‑150 (MQ=255)
                                                                                                                                                  cAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTTGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATtcatca     <  1:940953/150‑1 (MQ=255)
                                                                                                                                                   aaGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTcatcat    >  1:433856/1‑150 (MQ=255)
                                                                                                                                                   aaGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGAATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTcatcat    >  1:998072/1‑150 (MQ=255)
                                                                                                                                                     gATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTg  >  1:536703/1‑150 (MQ=255)
                                                                                                                                                     gATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCAGTGTAGTTATTTTCCGCAGCACAGGCGATGTCATAAGTGAATTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATGAttt  >  2:1229374/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTG  >  NC_000913/4586465‑4586763

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