Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,241 A→T 5.8% intergenic (‑221/+26) ydiV ← / ← nlpC anti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,2410AT 5.8% 210.5 / 2.8 69intergenic (‑221/+26)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATAC  >  NC_000913/1792096‑1792382
                                                                                                                                                 |                                                                                                                                             
gAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTTAGGAGTAAATGGACGGCGATGTatca                                                                                                                                            <  1:48398/149‑1 (MQ=255)
    gAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                                                          >  6:340011/1‑147 (MQ=255)
    gAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                                                          <  5:340011/147‑1 (MQ=255)
       tAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGt                                                                                                                                     >  5:64552/1‑149 (MQ=255)
        aaCTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTc                                                                                                                                    >  2:110057/1‑149 (MQ=255)
           tttCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                                                                 >  5:197632/1‑149 (MQ=255)
               aTTTGTTCCATGTTAACCACTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                                                            <  4:39498/149‑1 (MQ=255)
                tttGTGCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                                                            <  5:7027/149‑1 (MQ=255)
                   gTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGa                                                                                                                          <  3:293159/148‑1 (MQ=255)
                      ccATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                                                                     >  6:97565/1‑134 (MQ=255)
                      ccATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                                                                     <  5:97565/134‑1 (MQ=255)
                       cATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATAc                                                                                                                     >  2:283308/1‑149 (MQ=255)
                              aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                                                             <  7:49244/134‑1 (MQ=255)
                              aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAttttt                                                                                                                             >  8:49244/1‑134 (MQ=255)
                              aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTc                                                                                                              <  7:142058/149‑1 (MQ=255)
                              aCCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTc                                                                                                              <  2:154589/149‑1 (MQ=255)
                               ccACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                                                             >  7:15316/1‑149 (MQ=255)
                                cACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGt                                                                                                            <  4:96880/149‑1 (MQ=255)
                                   tttttCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                         >  7:273676/1‑149 (MQ=255)
                                        cAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                                                                     <  1:110057/148‑1 (MQ=255)
                                         aGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGc                                                                                                   >  2:180363/1‑149 (MQ=255)
                                           gggTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                                                                 <  3:301110/149‑1 (MQ=255)
                                             gTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAg                                                                                                >  1:93863/1‑148 (MQ=255)
                                              tCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATAGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGa                                                                                               >  1:331866/1‑148 (MQ=255)
                                                cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGaaa                                                                                             <  7:296933/148‑1 (MQ=255)
                                                cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                                                            <  1:203690/149‑1 (MQ=255)
                                                cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                                                            >  2:350825/1‑149 (MQ=255)
                                                 tttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtt                                                                                           >  3:243982/1‑149 (MQ=255)
                                                   tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                                                                         >  4:190895/1‑149 (MQ=255)
                                                    ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                                                                         <  4:66311/148‑1 (MQ=255)
                                                         gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                                                                    <  7:315462/148‑1 (MQ=255)
                                                         aaTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCg                                                                                                                                      >  8:278512/2‑98 (MQ=39)
                                                         aaTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCg                                                                                                                                      <  7:278512/97‑1 (MQ=39)
                                                           ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCaa                                                                                 <  2:93863/149‑1 (MQ=255)
                                                             aTCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCaa                                                                                 >  7:294156/1‑147 (MQ=255)
                                                              tCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTTCTGTATTCAGGAGTAAATGGACGGCGATGTAACATCGCCGTCCATTTTTCGACATACTAGATTCGTCGCGCCTGCCAGAAATTTTTCAGCCaa                                                                                 >  6:325974/1‑146 (MQ=255)
                                                                   aTTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAtt                                                                         <  6:64552/149‑1 (MQ=255)
                                                                      ataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                                                                       <  8:347024/148‑1 (MQ=255)
                                                                            gAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                                                                  >  1:78298/1‑147 (MQ=255)
                                                                            gAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                                                                <  2:78298/149‑1 (MQ=255)
                                                                             aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                                                               >  3:211929/1‑149 (MQ=255)
                                                                                 atCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                                                           <  8:338317/149‑1 (MQ=255)
                                                                                  tCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                                                          >  7:16789/1‑149 (MQ=255)
                                                                                  tCCATTTACGGAATTTTTGTCTGCAAAATACTACGGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATAAAGTGAGGAACg                                                          >  1:64123/1‑149 (MQ=255)
                                                                                   ccATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                                                          >  6:10411/1‑148 (MQ=255)
                                                                                     aTTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                         <  7:123075/99‑1 (MQ=255)
                                                                                     aTTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                         >  8:123075/1‑99 (MQ=255)
                                                                                         aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc                                                         >  2:19407/1‑143 (MQ=255)
                                                                                         aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc                                                         <  1:19407/143‑1 (MQ=255)
                                                                                           ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCTTCTCGCCTGACAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAAGCATACCTcc                                                 >  5:41440/1‑149 (MQ=255)
                                                                                           ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTcc                                                 <  5:10411/149‑1 (MQ=255)
                                                                                           ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTcc                                                 <  6:197632/149‑1 (MQ=255)
                                                                                            gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTccc                                                <  5:183366/149‑1 (MQ=255)
                                                                                                ttttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                                                                          >  1:207942/1‑119 (MQ=255)
                                                                                                ttttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                                                                          <  2:207942/119‑1 (MQ=255)
                                                                                                   tGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTg                                          <  8:15316/148‑1 (MQ=255)
                                                                                                         cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAgg                                   >  1:350903/1‑149 (MQ=255)
                                                                                                            aaTACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGt                                <  1:180363/149‑1 (MQ=255)
                                                                                                                ctactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGAt                            <  3:190895/149‑1 (MQ=255)
                                                                                                                        ttCAGGAGTAAATGGAGGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGAAGAGGCGTGGATAAATTgg                     <  4:243982/148‑1 (MQ=255)
                                                                                                                           aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTgg                     <  1:20376/145‑1 (MQ=255)
                                                                                                                           aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTgg                     >  2:20376/1‑145 (MQ=255)
                                                                                                                           aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGgttgt                 >  7:216988/1‑149 (MQ=255)
                                                                                                                           aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCGTGCTGGTAGAGGCGTGGATAAATTGgttgt                 >  3:205603/1‑149 (MQ=255)
                                                                                                                            ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGgttgtt                <  7:68068/149‑1 (MQ=255)
                                                                                                                             gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGttgttg               <  8:294156/149‑1 (MQ=255)
                                                                                                                                  aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATc          <  2:64123/149‑1 (MQ=255)
                                                                                                                                  aaTGGACGGCGATGTATCATCGCCGGCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATc          <  2:350903/149‑1 (MQ=255)
                                                                                                                                      gACGGCGATGAATCATCGCCGTCCATTTTTCGCCATACTAGCTTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAATTGAGGAACGCCGAACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTTGTATCgga       >  2:266874/1‑146 (MQ=255)
                                                                                                                                       aCGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCaa                                                                    >  8:203136/1‑86 (MQ=255)
                                                                                                                                        cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTaaa     >  3:16209/1‑148 (MQ=255)
                                                                                                                                          gagATGAATCATCGCCGTCCCTTTTTCGACATACTAGATTCGTATCGGCTTCCAGAAATTTTTCTGCCAATAGACATTATGAAGTTAGGAACGCATCACTCCCTTGCTGGTCGCGGCGTGGATAAATTGGTTGTTGGGATCGTAAATAc  <  1:266874/147‑1 (MQ=255)
                                                                                                                                                 atcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGt        >  6:200332/1‑136 (MQ=255)
                                                                                                                                                 atcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGt        <  5:200332/136‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATAC  >  NC_000913/1792096‑1792382

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