Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,665 G→A 7.6% G254G (GGG→GGA yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6650GA 7.6% 230.3 / 4.9 94G254G (GGG→GGAyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (27/55);  new base A (3/4);  total (33/61)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGA  >  NC_000913/2357516‑2357810
                                                                                                                                                     |                                                                                                                                                  
aaTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGAGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                                    >  2:959865/1‑150 (MQ=255)
 aTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGc                                                                                                                                                   <  2:516924/150‑1 (MQ=255)
  tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCa                                                                                                                                                  <  2:1334333/150‑1 (MQ=255)
  tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCa                                                                                                                                                  <  1:217014/150‑1 (MQ=255)
  tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCa                                                                                                                                                  >  1:598221/1‑150 (MQ=255)
    aCAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATc                                                                                                                                                <  2:546646/150‑1 (MQ=255)
    aCAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATc                                                                                                                                                <  1:664083/150‑1 (MQ=255)
     cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCa                                                                                                                                               <  2:1216143/150‑1 (MQ=255)
     cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCa                                                                                                                                               <  2:706153/150‑1 (MQ=255)
     cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCa                                                                                                                                               <  2:975577/150‑1 (MQ=255)
       gccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCaaa                                                                                                                                             >  2:586712/1‑150 (MQ=255)
       gccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTACACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCTTCaaa                                                                                                                                             >  1:279922/1‑150 (MQ=255)
        ccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAt                                                                                                                                            >  2:1066443/1‑150 (MQ=255)
         cagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAtt                                                                                                                                           >  1:593927/1‑150 (MQ=255)
         aagcTACAAACAATGTTTAATTATCTGCTGCAATTCGTAGATACCTCCCTTTTCACCCGTTTTATTGAGGAGATTTACGAACTTTCACAACTACAAAAGGAGAGATTAATGAATATTGCTGAACGTCTTCTGCAGGAAGGGAATAAAAtt                                                                                                                                           >  1:60715/2‑150 (MQ=255)
           gcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTTTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtgg                                                                                                                                         <  2:1140666/150‑1 (MQ=255)
           ccTACAAACCCTGATTAATTATCTGCTGCATTGCGGCGATACGTCCCGTTTGACCCGTTTAATTTAGGATATTGCCTAACGTTCACCACTAAAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtgg                                                                                                                                         <  2:1101191/149‑1 (MQ=21)
            cTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggc                                                                                                                                        <  2:534981/150‑1 (MQ=255)
             tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                                       <  2:1689451/150‑1 (MQ=255)
              aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                                      >  2:1637557/1‑150 (MQ=255)
                aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctggc                                                                                                                                    <  1:1053533/150‑1 (MQ=255)
                aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTACACCCGTTTTATTGAGGAGATCTCCGAACGTTCACCACTACAAAAGTAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATtggctggc                                                                                                                                    >  1:880132/1‑150 (MQ=38)
                  acacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAg                                                                                                                                  <  2:1193122/150‑1 (MQ=255)
                       gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                             >  2:1264964/1‑150 (MQ=255)
                       gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                             >  1:947703/1‑150 (MQ=255)
                       gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                             >  1:1619025/1‑150 (MQ=255)
                          tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTCTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaa                                                                                                                          <  2:780015/150‑2 (MQ=255)
                               aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGGGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATg                                                                                                                    >  1:441019/1‑150 (MQ=255)
                                  tgctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGCACGTTCACCACTACAAACGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                                 <  2:1571113/150‑1 (MQ=255)
                                   gctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACa                                                                                                                >  1:1106156/1‑150 (MQ=255)
                                    ctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTATGCTTGAACaa                                                                                                               >  1:410490/1‑150 (MQ=255)
                                     tgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAg                                                                                                              <  1:641672/150‑1 (MQ=255)
                                      gcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATTAATGaagcttt                                                                                                             >  1:1313942/1‑145 (MQ=255)
                                         aTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCGAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAt                                                                                                          <  1:1764115/150‑1 (MQ=255)
                                          tGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAtt                                                                                                         >  2:578708/1‑150 (MQ=255)
                                          tGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAtt                                                                                                         >  1:1478341/1‑150 (MQ=255)
                                            cggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGATAGAACAATCCATTaa                                                                                                       >  1:186813/1‑150 (MQ=255)
                                              gcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGACATTaaaa                                                                                                     >  1:693825/1‑150 (MQ=255)
                                                        ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                                           <  2:1702976/150‑1 (MQ=255)
                                                            ttCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTTCGAATg                                                                                       >  1:1004399/1‑150 (MQ=255)
                                                              aaCCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCt                                                                                     <  2:1278586/149‑1 (MQ=255)
                                                                 ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGa                                                                                  <  2:1494566/150‑1 (MQ=255)
                                                                   gTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGc                                                                                <  1:716668/150‑1 (MQ=255)
                                                                       tATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAAt                                                                                                                                            <  2:126869/87‑1 (MQ=255)
                                                                       tATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAAt                                                                                                                                            >  1:126869/1‑87 (MQ=255)
                                                                        aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGc                                                                           <  2:333723/150‑1 (MQ=255)
                                                                         ttGAGTAGATTTCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCt                                                                          <  1:1460619/150‑1 (MQ=255)
                                                                           gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAATACATTAAAAATGCTATTAGAATTCTGTAGCAggggatt                                                                        >  2:545364/1‑146 (MQ=21)
                                                                           gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAATCCATTAAAATTGCTTTTATCATTCAGGAGCAGAGCatt                                                                        >  2:1650355/1‑150 (MQ=25)
                                                                           gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTAATGCTGAACGGCTACGGCAGGAAGGGCATCAAAATGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTTCGAAAGCTGCAGCATCGCttt                                                                        >  2:102413/1‑150 (MQ=25)
                                                                            aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTg                                                                       <  1:1705185/150‑1 (MQ=255)
                                                                             ggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGa                                                                      <  2:863093/150‑1 (MQ=255)
                                                                               agaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATTAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                                    <  1:1816342/150‑1 (MQ=255)
                                                                               agaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                                    <  1:959865/150‑1 (MQ=255)
                                                                                gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                                   <  1:1066443/150‑1 (MQ=255)
                                                                                gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                                   <  2:317755/150‑1 (MQ=255)
                                                                                  ttGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTg                                                                 <  1:1758506/150‑1 (MQ=255)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGAGATCAGAATGGCTCGCGTCTTGTATAGCGTCAAGAATCCATTAAGATTGCATTGATCATGGTGGAGAATGGCTTTGATCGTgg                                                                >  1:1593677/1‑149 (MQ=12)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGTGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                <  2:1236876/150‑1 (MQ=255)
                                                                                   tGCCGAACGTTCACCACTACAAAACGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGGATCTAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                <  2:1645016/150‑1 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAAtt                                                                                                                                           <  2:79848/75‑1 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAAtt                                                                                                                                           >  1:79848/1‑75 (MQ=255)
                                                                                       gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGCCATTTAAATAGCTATGAGCATGCTTGAGCAGATCATTGAACTTgacaat                                                            >  1:628331/1‑146 (MQ=21)
                                                                                             tCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAACCTAGAACAAGAGATTAAAATTGCTATGCGAATGCTGGAGCAGGGCATTGAACCTGCCATtgatctg                                                      >  1:1170928/1‑145 (MQ=18)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACATGACATTAAAAATGCAATGAGCATTCTGGAGCAGAGCATTGAACGTGACATtgatcgtg                                                     >  2:418408/1‑144 (MQ=18)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATACATGAACAATACAATAAAATAGCATTGATCATTCTGTAGCACTGCATTGAACTTGACATTGAGCTgg                                                     >  2:185867/1‑150 (MQ=18)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAAAGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGTAGGTATGCATGAACACGCCATTAAAATTGCTTTGCGCATGCTGGAGCATGGCTTTGAACGTGAGATTGTGCTgg                                                     <  1:970373/150‑1 (MQ=255)
                                                                                                   atACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACa                                                <  1:94044/149‑1 (MQ=255)
                                                                                                          aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAg                                                                                                                              <  2:364384/66‑1 (MQ=255)
                                                                                                          aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAg                                                                                                                              >  1:364384/1‑66 (MQ=25)
                                                                                                          aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAAc                                         <  1:1637557/150‑1 (MQ=255)
                                                                                                               agaTTAATGACTAATGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACt                                    <  1:989603/150‑1 (MQ=255)
                                                                                                                  ttAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAg                                                                                                                                  >  1:397468/1‑54 (MQ=255)
                                                                                                                    aaTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGc                               <  1:579449/150‑1 (MQ=255)
                                                                                                                    aaTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGc                               <  1:770188/150‑1 (MQ=255)
                                                                                                                      tGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTg                             <  1:1693506/150‑1 (MQ=255)
                                                                                                                        aCTATTGCTGAACTGCTACGTCAGTAAGGTCATCAAATTGGCTGGCATGAAGGTATGCATGAAAAAGCCATTAAAATTGCTTTGCGCATGCTGGATAAGGGCTTTGAACGTTAGATTGTGCTGGCGACAACCCAACTCACTGATGCCGat                           <  2:60715/150‑1 (MQ=37)
                                                                                                                        aCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGat                           <  1:586712/150‑1 (MQ=255)
                                                                                                                        aCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGat                           <  1:1285212/150‑1 (MQ=255)
                                                                                                                        aCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGat                           <  1:1325544/150‑1 (MQ=255)
                                                                                                                        aCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGat                           <  1:1526845/150‑1 (MQ=255)
                                                                                                                        aCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGat                           <  1:758338/150‑1 (MQ=255)
                                                                                                                           aTTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATAtt                        >  2:1199800/1‑150 (MQ=255)
                                                                                                                               cTGAACGGCTTCGGGTGGAAGGGCATCTAATTGGCTGGCAGGAAGGTTTGCATGAACTCGCCTTTAAAATTGCTTTGCGCATGCTGGAGCAGGTCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGa                    <  2:743083/150‑1 (MQ=37)
                                                                                                                               cTGAACGGCTACGGCAGGAAGGGCATCAACTTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGGAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCgt                    <  1:392751/150‑2 (MQ=255)
                                                                                                                                tGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATCGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa                   <  2:922919/150‑1 (MQ=255)
                                                                                                                                 gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAAc                  <  2:279922/150‑1 (MQ=255)
                                                                                                                                      gCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCa             <  1:1214323/150‑1 (MQ=255)
                                                                                                                                      gCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCa             <  1:923804/150‑1 (MQ=255)
                                                                                                                                      gCTAAGGAATGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAATCCATTAAAATTGCATTGCGCATGCTGTAGCAGGGCTTTGTACGTGTGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCa             <  1:875482/150‑1 (MQ=38)
                                                                                                                                       cTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAt            <  1:1264964/150‑1 (MQ=255)
                                                                                                                                        tACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAtt           <  1:1547659/150‑1 (MQ=255)
                                                                                                                                        tACGGCAGGAAGGGCATCAAATTGGCGGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAtt           <  2:1739120/150‑1 (MQ=255)
                                                                                                                                              aGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAgg     <  2:1106156/150‑1 (MQ=255)
                                                                                                                                                 aaGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgaga  <  1:1199800/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGA  >  NC_000913/2357516‑2357810

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: