Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,592,199 C→T 6.7% A181A (GCC→GCT tsr → methyl‑accepting chemotaxis protein I, serine sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,1990CT 6.7% 298.9 / 9.2 90A181A (GCC→GCTtsrmethyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base C (42/42);  new base T (3/3);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGC  >  NC_000913/4592068‑4592340
                                                                                                                                   |                                                                                                                                             
gAGCTGATGCAACTCTTAGGTGCAGGCAAAAGCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGCTTCGAGAAGCAATATGTAGCTTACATGGAGGAAACCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATg                                                                                                                               <  3:956177/148‑1 (MQ=255)
gAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATg                                                                                                                               <  4:1106457/148‑1 (MQ=255)
                            aaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                 >  2:710650/1‑118 (MQ=255)
                            aaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                 <  1:710650/118‑1 (MQ=255)
                                   cGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTgg                                                                                                        <  2:612960/136‑1 (MQ=255)
                                   cGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTgg                                                                                                        >  1:612960/1‑136 (MQ=255)
                                     aGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGt                                                                                          >  2:1116699/1‑148 (MQ=255)
                                       ttCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                 >  4:86920/1‑107 (MQ=255)
                                       ttCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                 <  3:86920/107‑1 (MQ=255)
                                           ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                 <  4:827650/103‑1 (MQ=255)
                                           ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACaa                                                                                                                                 >  3:827650/1‑103 (MQ=255)
                                           ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTgatgat                                                                                    >  4:1095539/1‑148 (MQ=255)
                                           ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTgatgat                                                                                    >  1:887222/1‑148 (MQ=255)
                                              aTCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgt                                                                                 <  2:887222/148‑1 (MQ=255)
                                                         cAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgtg                                                                                                         <  2:418637/113‑1 (MQ=255)
                                                         cAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGAtgtg                                                                                                         >  1:418637/1‑113 (MQ=255)
                                                           gggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                             >  1:322870/1‑91 (MQ=255)
                                                           gggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                             <  2:322870/91‑1 (MQ=255)
                                                           gggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTAc                                                                                                                       <  1:882215/97‑1 (MQ=255)
                                                           gggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTAc                                                                                                                       >  2:882215/1‑97 (MQ=255)
                                                           gggATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTc                                                                    >  4:307586/1‑148 (MQ=255)
                                                                atCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCtt                                                               <  1:1116699/148‑1 (MQ=255)
                                                                atCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCtt                                                               <  3:307586/148‑1 (MQ=255)
                                                                    ggACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGt                                                                     <  3:728462/138‑1 (MQ=255)
                                                                    ggACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGt                                                                     >  2:298704/1‑138 (MQ=255)
                                                                    ggACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGt                                                                     <  1:298704/138‑1 (MQ=255)
                                                                    ggACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGcc                                                           >  1:289282/1‑148 (MQ=255)
                                                                    ggACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACACCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGt                                                                     >  4:728462/1‑138 (MQ=255)
                                                                     gACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                               <  3:621721/111‑1 (MQ=255)
                                                                     gACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCtggtg                                                                                               >  4:621721/1‑111 (MQ=255)
                                                                     gACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCa                                                                   >  3:193240/1‑139 (MQ=255)
                                                                     gACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCa                                                                   <  4:193240/139‑1 (MQ=255)
                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGACGGCGATATCatgg                                                                                                                               <  1:1020716/78‑4 (MQ=37)
                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGACGGCGATATCatgg                                                                                                                               >  2:1020716/1‑75 (MQ=38)
                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGACGGCGATATCatgg                                                                                                                               <  1:187632/78‑4 (MQ=38)
                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCTGACGGCGATATCatgg                                                                                                                               >  2:187632/1‑75 (MQ=37)
                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCg                                                                                           <  3:217174/114‑1 (MQ=255)
                                                                      aCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCg                                                                                           >  4:217174/1‑114 (MQ=255)
                                                                           ttCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCa                                                                   <  4:53698/133‑1 (MQ=255)
                                                                           ttCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCa                                                                   >  3:53698/1‑133 (MQ=255)
                                                                                       tATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                             >  4:578503/1‑63 (MQ=255)
                                                                                       tATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGcc                                                                                                                             <  3:578503/63‑1 (MQ=255)
                                                                                       tATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAg                                        <  3:1095539/148‑1 (MQ=255)
                                                                                            ggCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgt                                                                                 <  3:151426/102‑1 (MQ=255)
                                                                                            ggCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGAtcgt                                                                                 >  4:151426/1‑102 (MQ=255)
                                                                                                tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATaa                                                                                                                                    >  4:886894/1‑47 (MQ=255)
                                                                                                tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATaa                                                                                                                                    <  3:886894/47‑1 (MQ=255)
                                                                                                tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTg                                                                          >  3:736686/1‑105 (MQ=255)
                                                                                                tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTg                                                                          <  4:736686/105‑1 (MQ=255)
                                                                                                tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTgg                               >  1:140021/1‑148 (MQ=255)
                                                                                                 aCATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGt                              <  2:289282/148‑1 (MQ=255)
                                                                                                    tGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAtt                                                                                                     <  1:615332/74‑1 (MQ=255)
                                                                                                    tGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAtt                                                                                                     >  2:615332/1‑74 (MQ=255)
                                                                                                     ggAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAgcg                          >  3:238341/1‑148 (MQ=255)
                                                                                                      gAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAgcgc                         >  4:822908/1‑148 (MQ=255)
                                                                                                       aGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCtt                                                               <  2:132971/109‑1 (MQ=255)
                                                                                                       aGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCtt                                                               >  1:132971/1‑109 (MQ=255)
                                                                                                        gCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCa                       >  3:1031847/1‑148 (MQ=255)
                                                                                                         cAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaa                                          <  3:433500/128‑1 (MQ=255)
                                                                                                         cAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTaa                                          >  4:433500/1‑128 (MQ=255)
                                                                                                           aaaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATg                    >  4:834085/1‑148 (MQ=255)
                                                                                                            aaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTc                                                  >  1:809327/1‑117 (MQ=255)
                                                                                                            aaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTc                                                  <  2:809327/117‑1 (MQ=255)
                                                                                                            aaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAg                                        >  3:680436/1‑127 (MQ=255)
                                                                                                            aaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAg                                        <  4:680436/127‑1 (MQ=255)
                                                                                                            aaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGa                   <  3:822908/148‑1 (MQ=255)
                                                                                                             aCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCg                                                 >  1:880850/1‑117 (MQ=255)
                                                                                                             aCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCg                                                 <  2:880850/117‑1 (MQ=255)
                                                                                                              cGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCt                                                                                                   <  2:491335/66‑1 (MQ=255)
                                                                                                              cGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCt                                                                                                   >  1:491335/1‑66 (MQ=255)
                                                                                                                  cGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAgg                                                                                                                <  3:907682/49‑1 (MQ=255)
                                                                                                                  cGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAgg                                                                                                                >  4:907682/1‑49 (MQ=255)
                                                                                                                  cGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGcc             <  3:834085/148‑1 (MQ=255)
                                                                                                                    gCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCg                                                 >  3:843440/1‑110 (MQ=255)
                                                                                                                    gCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTGCGCCGTCTGGTTCg                                                 <  4:843440/110‑1 (MQ=255)
                                                                                                                      tCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGAt         <  2:140021/148‑1 (MQ=255)
                                                                                                                      tCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGAt         <  4:1031847/148‑1 (MQ=255)
                                                                                                                         ttgttatCGCTGACGGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCg                                                                       <  1:332092/79‑1 (MQ=255)
                                                                                                                         ttgttatCGCTGACGGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCg                                                                       >  2:332092/5‑83 (MQ=255)
                                                                                                                         aTGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGt                                               >  4:501311/1‑107 (MQ=255)
                                                                                                                         aTGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGt                                               <  3:501311/107‑1 (MQ=255)
                                                                                                                         aTGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGa      >  1:292285/1‑148 (MQ=255)
                                                                                                                         aTGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGa      <  4:238341/148‑1 (MQ=255)
                                                                                                                          tGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGt                                                    <  4:993503/101‑1 (MQ=255)
                                                                                                                          tGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGt                                                    >  3:993503/1‑101 (MQ=255)
                                                                                                                             tatCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGc  >  4:453893/1‑148 (MQ=255)
                                                                                                                               tCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTc                                                        <  4:664818/92‑1 (MQ=255)
                                                                                                                               tCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTc                                                        >  3:664818/1‑92 (MQ=255)
                                                                                                                               tCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTg       >  2:1063497/1‑141 (MQ=255)
                                                                                                                               tCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTg       <  1:1063497/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                             
GAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGC  >  NC_000913/4592068‑4592340

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