Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,323,312 G→C 5.6% intergenic (‑82/+24) phnE ← / ← phnD defective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily/phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,3120GC 5.6% 256.7 / 5.6 71intergenic (‑82/+24)phnE/phnDdefective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily/phosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (2/2);  total (36/35)
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.45e-01

GGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCC  >  NC_000913/4323171‑4323455
                                                                                                                                             |                                                                                                                                               
ggCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTc                                                                                                                                          <  1:1121885/149‑1 (MQ=255)
ggCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTc                                                                                                                                          <  1:1301905/149‑1 (MQ=255)
   ccAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCCCGGCGTcgca                                                                                                                                       <  3:996121/149‑4 (MQ=255)
      gCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGt                                                                                                                                           <  2:944754/142‑1 (MQ=255)
      gCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGt                                                                                                                                           >  1:944754/1‑142 (MQ=255)
           gAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCt                                                                                                                               >  4:105931/1‑149 (MQ=255)
              gCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCtt                                                                                                                                        >  3:220561/1‑137 (MQ=255)
              gCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCtt                                                                                                                                        <  4:220561/137‑1 (MQ=255)
               cGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGa                                                                                                                           >  4:750726/1‑149 (MQ=255)
                 agaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAtt                                                                                                                         <  1:297093/149‑1 (MQ=255)
                 agaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAtt                                                                                                                         >  2:297093/1‑149 (MQ=255)
                   aaCCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCgggg                                                                                                                                              <  4:76923/126‑1 (MQ=255)
                   aaCCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCgggg                                                                                                                                              >  3:76923/1‑126 (MQ=255)
                   aaCCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTa                                                                                                                       <  2:352814/149‑1 (MQ=255)
                          tGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCAcc                                                                                                                >  4:342339/1‑149 (MQ=255)
                           gcgcTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCg                                                                                                               >  1:97649/1‑149 (MQ=255)
                           gcgcTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCg                                                                                                               >  3:1042885/1‑149 (MQ=255)
                             gcTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGTTTTACTGCACCGCt                                                                                                             <  3:105931/149‑1 (MQ=255)
                                 gggtgggGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCggg                                                                                                                                               >  3:1133147/1‑111 (MQ=255)
                                 gggtgggGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCggg                                                                                                                                               <  4:1133147/111‑1 (MQ=255)
                                     ggggCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCggg                                                                                                                                               <  3:81804/107‑1 (MQ=255)
                                     ggggCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCggg                                                                                                                                               >  4:81804/1‑107 (MQ=255)
                                            tggtgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTa                                                                                                                                       <  3:164163/108‑1 (MQ=255)
                                            tggtgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTa                                                                                                                                       >  4:164163/1‑108 (MQ=255)
                                               tgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATc                                                                                           >  2:735179/1‑149 (MQ=255)
                                                gatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCAc                                                                                                                 <  3:38195/126‑1 (MQ=255)
                                                gatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCAc                                                                                                                 >  4:38195/1‑126 (MQ=255)
                                                     tttGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCtt                                                                                     >  3:924921/1‑149 (MQ=255)
                                                        gCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAgc                                                                                  <  4:924921/149‑1 (MQ=255)
                                                         cATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAgcg                                                                                 >  3:1403603/1‑149 (MQ=255)
                                                           tGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCAc                                                                                                     <  4:163813/127‑1 (MQ=255)
                                                           tGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCAc                                                                                                     >  3:163813/1‑127 (MQ=255)
                                                                  cTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCa                                                                        >  1:1019162/1‑149 (MQ=255)
                                                                   tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgt                                                                              >  2:282426/1‑142 (MQ=255)
                                                                   tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgt                                                                              <  1:282426/142‑1 (MQ=255)
                                                                    ccGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                      <  4:1042885/149‑1 (MQ=255)
                                                                    ccGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                      <  2:97649/149‑1 (MQ=255)
                                                                            tAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                              <  3:342339/149‑1 (MQ=255)
                                                                            tAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                              <  1:735179/149‑1 (MQ=255)
                                                                             aGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                             <  3:70963/133‑1 (MQ=255)
                                                                             aGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                             >  4:70963/1‑133 (MQ=255)
                                                                                    aTAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTAc                                                                                                                      <  3:615242/85‑1 (MQ=255)
                                                                                    aTAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTAc                                                                                                                      >  4:615242/1‑85 (MQ=255)
                                                                                           gCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCa                                                                        >  1:402910/1‑124 (MQ=255)
                                                                                           gCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCa                                                                        <  2:402910/124‑1 (MQ=255)
                                                                                             ttGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCt                                             <  3:750726/149‑1 (MQ=255)
                                                                                                     cATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCa                                                               >  3:722844/1‑123 (MQ=255)
                                                                                                     cATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCa                                                               <  4:722844/123‑1 (MQ=255)
                                                                                                              tCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGttttt                            <  4:1403603/149‑1 (MQ=255)
                                                                                                                     cTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAAc                                                                                                >  3:541182/1‑74 (MQ=255)
                                                                                                                     cTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAAc                                                                                                <  4:541182/74‑1 (MQ=255)
                                                                                                                     cTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCt                                             >  3:237057/1‑125 (MQ=255)
                                                                                                                     cTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCt                                             <  4:237057/125‑1 (MQ=255)
                                                                                                                      tGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATc                                                                                           >  4:962314/1‑78 (MQ=255)
                                                                                                                      tGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATc                                                                                           <  3:962314/78‑1 (MQ=255)
                                                                                                                            ttGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTTCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTg              >  3:973256/1‑149 (MQ=255)
                                                                                                                             tGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGc             <  2:1019162/149‑1 (MQ=255)
                                                                                                                               tCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCggt                                   <  2:66124/125‑1 (MQ=255)
                                                                                                                               tCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCggt                                   >  1:66124/1‑125 (MQ=255)
                                                                                                                                cGGATAAGACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGaa                                          >  3:378052/1‑117 (MQ=38)
                                                                                                                                  gataagACGCCCCGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTgc                                                 <  3:630981/103‑1 (MQ=255)
                                                                                                                                  gataagACGCCCCGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTgc                                                 >  4:630981/6‑108 (MQ=255)
                                                                                                                                  gATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCgg                                                                   >  3:624028/1‑90 (MQ=255)
                                                                                                                                  gATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCgg                                                                   <  4:624028/90‑1 (MQ=255)
                                                                                                                                  gATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCa                        <  1:1372386/133‑1 (MQ=255)
                                                                                                                                  gATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCa                        >  2:1372386/1‑133 (MQ=255)
                                                                                                                                      cGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTTACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGt    >  2:1338212/1‑149 (MQ=255)
                                                                                                                                      cGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAAc                                                                                                <  1:360986/57‑1 (MQ=255)
                                                                                                                                      cGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAAc                                                                                                >  2:360986/1‑57 (MQ=255)
                                                                                                                                       gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                              <  1:1292366/90‑1 (MQ=255)
                                                                                                                                       gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                              >  2:1292366/1‑90 (MQ=255)
                                                                                                                                       gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAg                       >  4:418834/1‑129 (MQ=255)
                                                                                                                                       gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAg                       <  3:418834/129‑1 (MQ=255)
                                                                                                                                        aCGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTcc  >  2:80652/1‑149 (MQ=255)
                                                                                                                                        aCGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTcc  >  3:668616/1‑149 (MQ=255)
                                                                                                                                          gCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTa                                                                                                          >  2:1293676/1‑43 (MQ=255)
                                                                                                                                          gCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTa                                                                                                          <  1:1293676/43‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCC  >  NC_000913/4323171‑4323455

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