Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,633,943 C→T P108S (CCA→TCA)  pinQ → Qin prophage; putative site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,633,9430CT98.9% 205.2 / NA 91P108S (CCA→TCA) pinQQin prophage; putative site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (50/40);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCT  >  NC_000913/1633809‑1634071
                                                                                                                                      |                                                                                                                                
gtgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACa                                                                                                                             >  2:299757/1‑140 (MQ=17)
 tgATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAg                                                                                                                            <  1:712236/140‑1 (MQ=18)
   aTCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGa                                                                                                                          >  2:203661/1‑140 (MQ=17)
    tCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGaa                                                                                                                         >  1:37451/1‑140 (MQ=17)
    tCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGaa                                                                                                                         <  2:610016/140‑1 (MQ=18)
     cAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGaaa                                                                                                                        >  1:235354/1‑140 (MQ=17)
      aaTTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGaaaa                                                                                                                       >  1:200200/1‑140 (MQ=17)
        ttgattgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAat                                                                                                                     >  2:598536/1‑140 (MQ=17)
        ttgattgTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAat                                                                                                                     >  1:358733/1‑140 (MQ=17)
             tgtgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAatgatga                                                                                                                >  1:697086/1‑140 (MQ=17)
               tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAatgatgatg                                                                                                              >  2:461098/1‑140 (MQ=17)
               tgACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAatgatgatg                                                                                                              <  1:87722/140‑1 (MQ=18)
                     aaCTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTa                                                                                                        >  2:368402/1‑140 (MQ=17)
                            tCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTcag                                                                                                 <  2:235354/140‑1 (MQ=18)
                               ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTcagcag                                                                                              >  1:199063/1‑140 (MQ=17)
                               ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTcagcag                                                                                              >  1:527621/1‑140 (MQ=17)
                               ccTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTcagcag                                                                                              >  1:162886/1‑140 (MQ=17)
                                 ttGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTc                                                                                            >  1:695602/1‑140 (MQ=17)
                                  tGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCg                                                                                           <  1:598536/140‑1 (MQ=18)
                                    gTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCt                                                                                         <  1:208502/140‑1 (MQ=18)
                                          aTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAAttt                                                                                   <  2:384793/140‑1 (MQ=17)
                                          aTGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTg                                                                                  >  1:126351/1‑140 (MQ=12)
                                               aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                              <  1:368402/140‑1 (MQ=17)
                                               aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                              <  2:200200/140‑1 (MQ=17)
                                               aTGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCAATGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACg                                                                              >  1:277286/1‑140 (MQ=17)
                                                tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                             >  1:256778/1‑140 (MQ=17)
                                                tGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACga                                                                             >  2:614877/1‑140 (MQ=17)
                                                      tCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCtt                                                                       >  2:247814/1‑140 (MQ=17)
                                                            aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACtt                                                                 >  2:258152/1‑140 (MQ=17)
                                                            aaaCAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACtt                                                                 >  2:107690/1‑140 (MQ=17)
                                                               cAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGaa                                                              >  2:308100/1‑140 (MQ=17)
                                                               cAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGaa                                                              <  2:37451/140‑1 (MQ=17)
                                                                 gTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACg                                                            >  1:377733/1‑140 (MQ=18)
                                                                  tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                           >  2:469104/1‑140 (MQ=18)
                                                                  tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                           >  1:124453/1‑140 (MQ=18)
                                                                  tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                           >  1:97941/1‑140 (MQ=18)
                                                                  tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGc                                                           >  2:720792/1‑140 (MQ=17)
                                                                  tGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCa                                                          <  2:126351/140‑1 (MQ=12)
                                                                   ggAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                             <  2:495562/137‑1 (MQ=17)
                                                                   ggAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAAc                                                             >  1:495562/1‑137 (MQ=18)
                                                                    gAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCAc                                                         >  2:168285/1‑140 (MQ=18)
                                                                     aacaacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                        <  2:256778/140‑1 (MQ=17)
                                                                        aacTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAt                                                     <  2:124453/140‑1 (MQ=17)
                                                                         acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                    <  2:358733/140‑1 (MQ=17)
                                                                         acTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                    >  2:221173/1‑140 (MQ=18)
                                                                          cTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTc                                                   <  1:168285/140‑1 (MQ=17)
                                                                            gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACTCACTCATTCCg                                                 >  2:328845/1‑140 (MQ=17)
                                                                            gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                        >  1:606149/1‑133 (MQ=18)
                                                                            gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACt                                                        <  2:606149/133‑1 (MQ=17)
                                                                            gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                    >  1:215029/1‑137 (MQ=18)
                                                                            gACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCAtt                                                    <  2:215029/137‑1 (MQ=17)
                                                                               cGATACAGGTATCAGAGTGCATTGGTTAGCATTGGGGGGCAGTGAACTGACCAGTCCAACAGGAAAAATGATGATGGAAGTAATTTAAGCAGTCGCTGAATTGGAACGAGACCGTTTACTTGAACGCACTCATTCCGGGa                                              <  2:566158/140‑1 (MQ=11)
                                                                               cGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                              <  1:469104/140‑1 (MQ=17)
                                                                               cGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                              <  2:441104/140‑1 (MQ=17)
                                                                               cGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGa                                              >  1:507990/1‑140 (MQ=18)
                                                                                  aaCAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAg                                           >  1:245083/1‑140 (MQ=18)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  2:507990/140‑1 (MQ=17)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  1:461098/140‑1 (MQ=17)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 >  1:522351/1‑140 (MQ=18)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  2:461467/140‑1 (MQ=17)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  2:527621/140‑1 (MQ=17)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  2:555855/140‑1 (MQ=17)
                                                                                            agagTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  2:621793/140‑1 (MQ=17)
                                                                                               gTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcgg                              >  1:559083/1‑140 (MQ=18)
                                                                                                tGCATTGCTTAGCATTGGTGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggc                             >  2:231358/1‑140 (MQ=17)
                                                                                                tGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 <  2:302731/136‑1 (MQ=17)
                                                                                                tGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCcg                                 >  1:302731/1‑136 (MQ=18)
                                                                                                 gCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCgcggcg                            >  2:465552/1‑140 (MQ=18)
                                                                                                  cATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGgcgc                           >  2:365978/1‑140 (MQ=18)
                                                                                                   aTTGCTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCa                          <  2:162886/140‑1 (MQ=17)
                                                                                                       cTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGa                      <  2:522351/140‑1 (MQ=17)
                                                                                                       cTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGa                      <  2:199063/140‑1 (MQ=17)
                                                                                                       cTTAGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGa                      <  1:308100/140‑1 (MQ=17)
                                                                                                          aGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc                   >  1:20091/1‑140 (MQ=18)
                                                                                                          aGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc                   >  1:29971/1‑140 (MQ=18)
                                                                                                          aGCATTGGGGGGCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAAc                   >  1:559429/1‑140 (MQ=18)
                                                                                                                 gggggCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg            >  2:44633/1‑140 (MQ=17)
                                                                                                                 gggggCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg            >  1:347437/1‑140 (MQ=17)
                                                                                                                 gggggCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg            >  2:589727/1‑140 (MQ=17)
                                                                                                                 gggggCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg            >  1:578146/1‑140 (MQ=17)
                                                                                                                 gggggCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg            >  2:247429/1‑140 (MQ=17)
                                                                                                                 gggggCATTGACCTGACCAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTgg            >  2:516128/1‑140 (MQ=17)
                                                                                                                         tgacctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac    <  1:516128/140‑1 (MQ=14)
                                                                                                                         tgacctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac    <  1:458517/140‑1 (MQ=14)
                                                                                                                         tgacctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac    <  2:263207/140‑1 (MQ=14)
                                                                                                                         tgacctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac    <  2:277286/140‑1 (MQ=14)
                                                                                                                         tgacctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac    <  2:29971/140‑1 (MQ=14)
                                                                                                                         tgacctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGaccac    <  1:247429/140‑1 (MQ=14)
                                                                                                                           acctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCt  <  2:559083/140‑1 (MQ=18)
                                                                                                                           acctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCt  <  2:97941/140‑1 (MQ=18)
                                                                                                                           acctgaccAGTTCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCt  <  2:347437/140‑1 (MQ=18)
                                                                                                                                      |                                                                                                                                
GTGATCAATTGATTGTGACAAAACTGGATCGCCTTGGTTGTAATGCAATGGATATCAGGAAAACAGTGGAACAACTGACCGAAACAGGTATCAGAGTGCATTGCTTAGCATTGGGGGGCATTGACCTGACCAGTCCAACAGGAAAAATGATGATGCAAGTAATTTCAGCAGTCGCTGAATTTGAACGAGACCTTTTACTTGAACGCACTCATTCCGGGATAGTAAGAGCCCGCGGCGCAGGGAAACGTTTTGGTCGACCACCT  >  NC_000913/1633809‑1634071

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