Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,218,001 T→C N359S (AAC→AGC)  aer ← aerotaxis sensor receptor, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,218,0010TC98.8% 297.1 / ‑4.9 85N359S (AAC→AGC) aeraerotaxis sensor receptor, flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (37/47);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGACAGTTTAT  >  NC_000913/3217866‑3218127
                                                                                                                                       |                                                                                                                              
gcagcaTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAAt                                                                                                                            <  1:396749/140‑1 (MQ=255)
gcagcaTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAAt                                                                                                                            <  2:100058/140‑1 (MQ=255)
   gcaTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCCAATCAg                                                                                                                         <  1:1206337/140‑1 (MQ=255)
    caTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAg                                                                                                                         >  2:1017023/1‑139 (MQ=255)
    caTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAg                                                                                                                         <  1:1017023/139‑1 (MQ=255)
     aTTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCg                                                                                                                       >  1:589338/1‑140 (MQ=255)
              gCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATgg                                                                                                              <  2:136528/140‑1 (MQ=255)
               cGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGt                                                                                                             >  2:1244950/1‑140 (MQ=255)
               cGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGt                                                                                                             >  1:1147459/1‑140 (MQ=255)
                   tggctAAATGACGCCCTCCCCCTGCCCCCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCa                                                                                                         <  2:577736/140‑1 (MQ=255)
                   tggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCa                                                                                                         <  1:1379457/140‑1 (MQ=255)
                   tggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCa                                                                                                         <  1:649329/140‑1 (MQ=255)
                    ggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAg                                                                                                                         <  2:350218/123‑1 (MQ=255)
                    ggctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAg                                                                                                                         >  1:350218/1‑123 (MQ=255)
                     gctAAATGACTCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCGGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAAt                                                                                                       <  2:1438076/140‑1 (MQ=255)
                     gctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAAt                                                                                                       <  2:177967/140‑1 (MQ=255)
                     gctAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAAt                                                                                                       <  1:186055/140‑1 (MQ=255)
                        aaaTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATgcg                                                                                                    <  1:1637112/140‑1 (MQ=255)
                        aaaTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATgcg                                                                                                    <  1:1194862/140‑1 (MQ=255)
                        aaaTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATgcg                                                                                                    <  2:137491/140‑1 (MQ=255)
                         aaTGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATgcgc                                                                                                   >  1:1116111/1‑140 (MQ=255)
                               gCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTa                                                                                             <  1:1229882/140‑1 (MQ=255)
                                   ttCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCa                                                                                                         >  2:474761/1‑124 (MQ=255)
                                   ttCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCa                                                                                                         <  1:474761/124‑1 (MQ=255)
                                     ccccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAAt                                                                                                       >  1:851758/1‑124 (MQ=255)
                                     ccccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAAt                                                                                                       <  2:851758/124‑1 (MQ=255)
                                     ccccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCg                                                                                       >  2:288592/1‑140 (MQ=255)
                                      cccTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCgg                                                                                      <  1:257971/140‑1 (MQ=255)
                                        cTGCCCCCCCTGAAAACCCTTTCCCCTGTTCGCCGGCACGCGCCGCTCCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCg                                                                                    <  2:1662136/140‑1 (MQ=37)
                                             accacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATc                                                                               >  2:1069600/1‑140 (MQ=255)
                                             accacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATc                                                                               <  1:1048771/140‑1 (MQ=255)
                                             accacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGACTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATc                                                                               >  1:1407879/1‑140 (MQ=255)
                                               cacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGt                                                                             >  2:935265/1‑140 (MQ=255)
                                               cacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGt                                                                             >  2:1134043/1‑140 (MQ=255)
                                               cacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGt                                                                             >  1:1537109/1‑140 (MQ=255)
                                               cacTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGt                                                                             >  1:1662136/1‑140 (MQ=255)
                                                  tGCAAACCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCa                                                                          <  2:1157431/140‑1 (MQ=255)
                                                  tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCa                                                                          <  2:589338/140‑1 (MQ=255)
                                                  tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCa                                                                          <  2:1355176/140‑1 (MQ=255)
                                                  tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCa                                                                          <  1:457665/140‑1 (MQ=255)
                                                  tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCa                                                                          <  2:1537109/140‑1 (MQ=255)
                                                  tGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCa                                                                          <  1:1069600/140‑1 (MQ=255)
                                                    cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCg                                                                                    >  2:1519130/1‑128 (MQ=255)
                                                    cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCg                                                                                    <  1:1519130/128‑1 (MQ=255)
                                                    cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCAtt                                                                        >  1:273801/1‑140 (MQ=255)
                                                    cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCAtt                                                                        >  1:465361/1‑140 (MQ=255)
                                                    cAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCAtt                                                                        <  2:781499/140‑1 (MQ=255)
                                                     aaaaCCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTg                                                                       >  2:1630345/1‑140 (MQ=255)
                                                         ccTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCtt                                                                   >  2:230259/1‑140 (MQ=255)
                                                             tGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATc                                                               <  1:507214/140‑1 (MQ=255)
                                                                   gTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCg                                                           >  1:793233/1‑138 (MQ=255)
                                                                   gTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCg                                                           <  2:793233/138‑1 (MQ=255)
                                                                    ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCtt                                                                   >  2:1420847/1‑129 (MQ=255)
                                                                    ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCtt                                                                   <  1:1420847/129‑1 (MQ=255)
                                                                    ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTgg                                                        >  2:295527/1‑140 (MQ=255)
                                                                    ttCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTgg                                                        >  2:518549/1‑140 (MQ=255)
                                                                        ccGGCCCGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCAt                                                    <  1:622770/140‑1 (MQ=255)
                                                                             aCGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCt                                               <  2:273801/140‑1 (MQ=255)
                                                                              cgcgCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTc                                              >  1:1190841/1‑140 (MQ=255)
                                                                              cgcgCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTc                                              >  2:1515767/1‑140 (MQ=255)
                                                                              cgcgCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTc                                              >  1:1511286/1‑140 (MQ=255)
                                                                               gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTcc                                             <  2:1179814/140‑1 (MQ=255)
                                                                               gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTcc                                             <  1:1515767/140‑1 (MQ=255)
                                                                               gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTcc                                             <  1:295527/140‑1 (MQ=255)
                                                                               gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTcc                                             <  2:1190841/140‑1 (MQ=255)
                                                                               gcgcCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATAGCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTcc                                             <  2:109948/140‑1 (MQ=255)
                                                                                      ttCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCt                                      >  2:470123/1‑140 (MQ=255)
                                                                                       tCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTg                                     <  2:1116111/140‑1 (MQ=255)
                                                                                               cgcATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCg                             >  2:586441/1‑140 (MQ=255)
                                                                                                      aGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACt                      >  1:962718/1‑140 (MQ=255)
                                                                                                      aGGGCCAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACt                      <  2:1147459/140‑1 (MQ=255)
                                                                                                          ccAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCa                  <  2:962718/140‑1 (MQ=255)
                                                                                                           cAGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAg                 <  1:586441/140‑1 (MQ=255)
                                                                                                            aGAATATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGt                >  1:561130/1‑140 (MQ=255)
                                                                                                              aaTATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTg                                >  2:480612/1‑122 (MQ=255)
                                                                                                              aaTATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTg                                <  1:480612/122‑1 (MQ=255)
                                                                                                              aaTATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGa              >  1:694548/1‑140 (MQ=255)
                                                                                                              aaTATTGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGa              <  1:518549/140‑1 (MQ=255)
                                                                                                                  ttGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGa          >  2:383408/1‑140 (MQ=255)
                                                                                                                  ttGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGa          >  2:802819/1‑140 (MQ=255)
                                                                                                                   tGGTCTGAAACGCAATATCGCTAATCCGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTTCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGAc         <  2:1468486/140‑1 (MQ=255)
                                                                                                                   tGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGAc         >  2:486858/1‑140 (MQ=255)
                                                                                                                   tGGTCTGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGAc         <  2:1511286/140‑1 (MQ=255)
                                                                                                                        tGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGACAGttt    >  2:1266678/1‑140 (MQ=255)
                                                                                                                        tGAAACGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCAGCATTACTGGCAGTGATTGACAGttt    >  2:325559/1‑140 (MQ=255)
                                                                                                                          aaaCGCAATATCGCTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCACACCCTGCACTGCCGCAATACTGGCAGTGATTGACAGTTTAt  <  1:689610/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                              
GCAGCATTAGCGCTGCGGCTGGCTAAATGACGCACTTCCCCTGCCACCACTGCAAAACCTTTGCCCTGTTCGCCGGCACGCGCCGCTTCCACCGCCGCATTCAGGGCCAGAATATTGGTCTGAAACGCAATATCGTTAATCAGCGAAGTAATGGTGCCAATGCGCTGGGTACTGTCGGCGATATCGTCCATTGTCTTGATCACCGTGGTCATCGCCTCCCCACCCTGCACTGCCGCATTACTGGCAGTGATTGACAGTTTAT  >  NC_000913/3217866‑3218127

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