Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,915 G→A 8.8% intergenic (‑92/+39) yddW ← / ← gadC liprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,9150GA 8.8% 392.6 / 22.1 127intergenic (‑92/+39)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (58/56);  new base A (5/6);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 127 total aligned reads displayed.

ACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACG  >  NC_000913/1568768‑1569062
                                                                                                                                                   |                                                                                                                                                    
aCTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGa                                                                                                                                                     >  1:1124054/1‑149 (MQ=255)
 cTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                    <  2:1124054/149‑1 (MQ=255)
 cTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                    <  4:566315/149‑1 (MQ=255)
  tAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                    >  3:772232/1‑148 (MQ=255)
  tAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                    <  4:772232/148‑1 (MQ=255)
      aTCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                              <  2:932917/149‑1 (MQ=255)
         gCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTC‑ATCGCt                                                                                                                                           <  1:1289796/149‑1 (MQ=255)
              tcttctAATCGTTAATTTCTTGTTACGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGACAATACGACAGCAAGCACGAGAAAGAGAGCGATGAATT‑AACGCTCCCtt                                                                                                                                      >  3:916065/1‑149 (MQ=255)
                 tctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑a                                                                                                                                                <  1:874405/136‑1 (MQ=255)
                   tAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAAt                                                                                                                                                   <  4:678227/132‑1 (MQ=255)
                      tCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTtata                                                                                                                              >  2:1338562/1‑149 (MQ=255)
                      tCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTC‑ATCGCTcc                                                                                                                                         >  3:619020/1‑138 (MQ=255)
                        gTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTta                                                                                                                                >  1:939394/1‑145 (MQ=255)
                        gTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTta                                                                                                                                <  2:939394/145‑1 (MQ=255)
                            aTTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCAtt                                                                                                                        >  4:825986/1‑149 (MQ=255)
                             tttCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTc                                                                                                                       >  3:824078/1‑149 (MQ=255)
                                cTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTc                                                                                                                                   <  3:1074049/134‑1 (MQ=255)
                                 ttGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑At                                                                                                                                               >  1:178993/1‑121 (MQ=255)
                                 ttGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑At                                                                                                                                               <  2:178993/121‑1 (MQ=255)
                                 ttGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGAc                                                                                                                   >  1:1308383/1‑149 (MQ=255)
                                     ttCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATgg                                                                                                               <  3:763318/149‑1 (MQ=255)
                                              aGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                    >  3:378612/1‑104 (MQ=255)
                                              aGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGaa                                                                                                                                                    <  4:378612/104‑1 (MQ=255)
                                                atatCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                    <  2:1308383/149‑1 (MQ=255)
                                                     cATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtca                                                                                               >  3:1089583/1‑149 (MQ=255)
                                                           tgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatc                                                                                         >  4:885213/1‑149 (MQ=255)
                                                            gCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                        <  4:824078/149‑1 (MQ=255)
                                                               ggTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                              >  2:1463239/1‑92 (MQ=255)
                                                               ggTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                              <  1:1463239/92‑1 (MQ=255)
                                                                gTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑At                                                                                                                                               <  2:852633/90‑1 (MQ=255)
                                                                      ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACa                                                                                                                  <  4:289041/113‑1 (MQ=255)
                                                                      ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACa                                                                                                                  >  3:289041/1‑113 (MQ=255)
                                                                       gTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattca                                                                                           <  2:1215903/135‑1 (MQ=255)
                                                                       gTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattca                                                                                           >  1:1215903/1‑135 (MQ=255)
                                                                        tAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTC‑ATCGCTc                                                                                                                                          <  3:562173/87‑1 (MQ=39)
                                                                        tAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTC‑ATCGCTc                                                                                                                                          >  4:562173/1‑87 (MQ=39)
                                                                         aaaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtgg                                                                           <  3:825986/149‑1 (MQ=255)
                                                                          aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                          <  3:885213/149‑1 (MQ=255)
                                                                           aCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtc                                                                                                >  2:89635/1‑126 (MQ=255)
                                                                           aCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtc                                                                                                <  1:89635/126‑1 (MQ=255)
                                                                           aCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggtg                                                                         >  3:142279/1‑149 (MQ=255)
                                                                            cAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGa                                                                        >  1:756167/1‑149 (MQ=255)
                                                                             aGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTc                                                                                                     >  2:989049/1‑119 (MQ=255)
                                                                             aGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGAT‑AATTCATCGCTc                                                                                                                                          <  3:870562/82‑1 (MQ=38)
                                                                                tgtgCGCTCATTCTCGTTCTTATCCCCGTTAAGTCAATACGACAGCAAGTACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGGTTCTTGTCATTCATCACAATATAGTGTGGTGAACGt                                                                    <  4:916065/149‑1 (MQ=255)
                                                                                tgtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGt                                                                    <  4:1089583/149‑1 (MQ=255)
                                                                                  tgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGttt                                                                                                      <  1:318354/113‑1 (MQ=255)
                                                                                  tgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGttt                                                                                                      >  2:318354/1‑113 (MQ=255)
                                                                                    cgcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTga                                                                                                 >  1:1129812/1‑115 (MQ=255)
                                                                                    cgcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAc                                                                <  1:545902/149‑1 (MQ=255)
                                                                                     gcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGt                                                                                                                                    <  4:87217/80‑1 (MQ=255)
                                                                                     gcTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTg                                                                                                  >  1:983810/1‑114 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAgt                                                                                                          >  2:1061652/1‑105 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAgt                                                                                                          <  1:1061652/105‑1 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                          <  3:129651/137‑1 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTATCTTGTCATTCATCACAATATAgtgt                                                                             <  2:270201/134‑1 (MQ=255)
                                                                                      cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTATCTTGTCATTCATCACAATATAgtgt                                                                             >  1:270201/1‑134 (MQ=255)
                                                                                       tCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGt                                                                                                 >  4:1416174/1‑113 (MQ=255)
                                                                                       tCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGt                                                                                                 <  3:1416174/113‑1 (MQ=255)
                                                                                        cATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCgg                                                            <  1:1338562/149‑1 (MQ=255)
                                                                                             tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTta                                                                                                                                >  3:223583/1‑76 (MQ=255)
                                                                                             tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTta                                                                                                                                <  4:223583/76‑1 (MQ=255)
                                                                                               gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTc                                                                                                     <  4:566523/101‑1 (MQ=255)
                                                                                               gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTc                                                                                                     >  3:566523/1‑101 (MQ=255)
                                                                                               gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                        >  1:1400219/1‑114 (MQ=255)
                                                                                               gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                        <  2:1400219/114‑1 (MQ=255)
                                                                                               gTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAgg                                                     >  2:904552/1‑149 (MQ=255)
                                                                                                   ttATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGc                                                                                                                                            <  3:1059114/58‑1 (MQ=255)
                                                                                                   ttATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGc                                                                                                                                            >  4:1059114/1‑58 (MQ=255)
                                                                                                     aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                              <  3:816326/54‑1 (MQ=255)
                                                                                                     aTCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                              >  4:816326/1‑54 (MQ=255)
                                                                                                      tCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGgaagaa                                                >  2:1141278/1‑147 (MQ=255)
                                                                                                      tCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                              <  2:756167/149‑1 (MQ=255)
                                                                                                       ccccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTg                                             >  4:249795/1‑149 (MQ=255)
                                                                                                        cccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGa                                            <  4:142279/149‑1 (MQ=255)
                                                                                                         ccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGAc                                           >  3:633820/1‑149 (MQ=255)
                                                                                                               aGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGaa                                                                       <  4:759591/115‑1 (MQ=255)
                                                                                                               aGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGaa                                                                       >  3:759591/1‑115 (MQ=255)
                                                                                                                 tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTgg                                    >  3:1108845/1‑148 (MQ=255)
                                                                                                                 tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTgg                                    <  4:1108845/148‑1 (MQ=255)
                                                                                                                    aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                    >  2:1102492/1‑81 (MQ=255)
                                                                                                                    aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCt                                                                                                    <  1:1102492/81‑1 (MQ=255)
                                                                                                                          cAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGa                          <  4:633820/149‑1 (MQ=255)
                                                                                                                                    gAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                >  1:831760/1‑149 (MQ=255)
                                                                                                                                       aaaGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGc              >  3:389134/1‑148 (MQ=255)
                                                                                                                                         aGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTg                                             <  4:1106889/115‑1 (MQ=255)
                                                                                                                                          gggAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa          >  1:371829/1‑149 (MQ=255)
                                                                                                                                           ggAGCGATAATTC‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTcc           <  3:1465797/147‑1 (MQ=255)
                                                                                                                                           ggAGCGAT‑AATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag         >  2:346101/1‑149 (MQ=255)
                                                                                                                                               cGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTg     <  2:831760/149‑1 (MQ=255)
                                                                                                                                                gATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGt                                                                                                              >  4:1421697/1‑43 (MQ=255)
                                                                                                                                                gATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTgg                                    >  3:1424788/1‑117 (MQ=255)
                                                                                                                                                gATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTgg                                    <  4:1424788/117‑1 (MQ=255)
                                                                                                                                                 aTGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGc              <  2:417668/138‑1 (MQ=255)
                                                                                                                                                 aTGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGc              >  1:417668/1‑138 (MQ=255)
                                                                                                                                                  tGAATT‑ATCGCTCCCTTGTCTTATAACCTTTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcgc                                                              >  3:95178/1‑89 (MQ=255)
                                                                                                                                                  tGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGt                                                          >  1:1428637/1‑93 (MQ=255)
                                                                                                                                                  tGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGt                                                          <  2:1428637/93‑1 (MQ=255)
                                                                                                                                                  tGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGa                                            >  4:690441/1‑107 (MQ=255)
                                                                                                                                                  tGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGa                                            <  3:690441/107‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTAGTATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACG  >  NC_000913/1568768‑1569062

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