Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,250,546 C→T intergenic (+170/‑16) yqjK → / → yqjF conserved protein YqjK/DoxX family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,250,5460CT99.6% 1022.0 / NA 252intergenic (+170/‑16)yqjK/yqjFconserved protein YqjK/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (108/143);  minor base G (1/0);  total (109/143)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

TAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTGGGGAAAAA  >  NC_000913/3250448‑3250643
                                                                                                  |                                                                                                 
tAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGgtgc                                                                                                 >  1:4232089/1‑101 (MQ=255)
 aaCAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCa                                                                                                <  1:4231698/101‑1 (MQ=255)
 aaCAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCa                                                                                                >  1:4232093/1‑101 (MQ=255)
  aCAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCa                                                                                                >  2:4232092/1‑100 (MQ=255)
  aCAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAg                                                                                               >  1:4232095/1‑101 (MQ=255)
    aaTCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCa                                                                                                <  2:4231795/98‑1 (MQ=255)
    aaTCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTg                                                                                             <  1:4231598/101‑1 (MQ=255)
    aaTCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTg                                                                                             <  2:4231893/101‑1 (MQ=255)
    aaTCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTg                                                                                             <  1:4231847/101‑1 (MQ=255)
     atcatcATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTgg                                                                                            <  1:4231624/101‑1 (MQ=255)
     atcatcATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTgg                                                                                            <  2:4231812/101‑1 (MQ=255)
     atcatcATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTgg                                                                                            <  2:4231618/101‑1 (MQ=255)
     atcatcATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTgg                                                                                            >  1:4232098/1‑101 (MQ=255)
      tcatcaTTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGga                                                                                           >  1:4232100/1‑101 (MQ=255)
       catcatTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGgag                                                                                          <  2:4231803/101‑1 (MQ=255)
        atcatTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGgaga                                                                                         <  1:4231767/101‑1 (MQ=255)
         tcatTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGagag                                                                                        >  2:4232104/1‑101 (MQ=255)
          catTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGt                                                                                       >  1:4232105/1‑101 (MQ=255)
            ttCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTat                                                                                     <  1:4231818/101‑1 (MQ=255)
            ttCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTat                                                                                     >  2:4232109/1‑101 (MQ=255)
             tCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTatg                                                                                    >  1:4232112/1‑101 (MQ=255)
             tCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTatg                                                                                    <  1:4231655/101‑1 (MQ=255)
               accacGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTatgat                                                                                  <  2:4231706/101‑1 (MQ=255)
                 cacGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGagag                                                                                        <  1:4231820/93‑1 (MQ=255)
                 cacGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTAtgatga                                                                                >  2:4232114/1‑101 (MQ=255)
                 cacGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTAtgatga                                                                                <  2:4231904/101‑1 (MQ=255)
                  acGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAATAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGaa                                                                               <  2:4231781/101‑1 (MQ=255)
                   cGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGaaa                                                                              <  1:4231744/101‑1 (MQ=255)
                     tttATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGaaaaa                                                                            >  1:4232117/1‑101 (MQ=255)
                     tttATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGaaaaa                                                                            <  2:4231670/101‑1 (MQ=255)
                      ttATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGaaaaaa                                                                           >  1:4232120/1‑101 (MQ=255)
                       tATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAAt                                                                          <  1:4232004/101‑1 (MQ=255)
                       tATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAAt                                                                          >  2:4232121/1‑101 (MQ=255)
                       tATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAAt                                                                          <  1:4231719/101‑1 (MQ=255)
                            ttCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTagaa                                                                     >  1:4232125/1‑101 (MQ=255)
                            ttCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTagaa                                                                     <  1:4231775/101‑1 (MQ=255)
                               tctcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGAt                                                                  >  1:4232129/1‑101 (MQ=255)
                               tctcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGAt                                                                  <  1:4231565/101‑1 (MQ=255)
                                ctcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATg                                                                 >  2:4232130/1‑101 (MQ=255)
                                 tcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGt                                                                <  1:4231680/101‑1 (MQ=255)
                                 tcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGt                                                                <  2:4231637/101‑1 (MQ=255)
                                 tcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGt                                                                >  2:4232134/1‑101 (MQ=255)
                                 tcCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGt                                                                >  1:4232133/1‑101 (MQ=255)
                                   cATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTg                                                              <  2:4231749/101‑1 (MQ=255)
                                     tCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgt                                                            <  1:4232046/101‑1 (MQ=255)
                                     tCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgt                                                            <  2:4231825/101‑1 (MQ=255)
                                     tCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgt                                                            <  1:4231854/101‑1 (MQ=255)
                                      cGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgtg                                                           <  2:4231842/101‑1 (MQ=255)
                                       gACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgtgt                                                          >  1:4232140/1‑101 (MQ=255)
                                       gACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgtgt                                                          >  2:4232139/1‑101 (MQ=255)
                                          aGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACt                                                       >  1:4232143/1‑101 (MQ=255)
                                          aGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACt                                                       >  2:4232145/1‑101 (MQ=255)
                                          aGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACt                                                       <  1:4231757/101‑1 (MQ=255)
                                           gCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTg                                                      <  2:4231961/101‑1 (MQ=255)
                                           gCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTg                                                      <  2:4231534/101‑1 (MQ=255)
                                           gCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTg                                                      <  1:4231896/101‑1 (MQ=255)
                                            cAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTgg                                                     >  1:4232147/1‑101 (MQ=255)
                                             aaCGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGt                                                    >  1:4232149/1‑101 (MQ=255)
                                               cgacgCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAg                                                  <  2:4231567/101‑1 (MQ=255)
                                                 acgCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAgcg                                                >  1:4232152/1‑101 (MQ=255)
                                                  cgCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGTCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGAt                                                                  <  1:4231441/82‑1 (MQ=255)
                                                  cgCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGgtg                                                           <  1:4231689/89‑1 (MQ=255)
                                                  cgCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAgcgc                                               <  1:4231877/101‑1 (MQ=255)
                                                    cTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGGGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAgcgcgc                                             >  1:4232154/1‑101 (MQ=255)
                                                         aCCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCAtttt                                        <  2:4231840/101‑1 (MQ=255)
                                                         aCCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCAtttt                                        >  1:4232158/1‑101 (MQ=255)
                                                           cgcgCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTaa                                      >  1:4232161/1‑101 (MQ=255)
                                                           cgcgCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTaa                                      <  2:4231720/101‑1 (MQ=255)
                                                             cgcCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTa                                       >  2:4232164/1‑98 (MQ=255)
                                                              gcCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGc                                   <  1:4231937/101‑1 (MQ=255)
                                                                cATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTa                                                                        <  1:4231858/62‑1 (MQ=255)
                                                                cATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCg                                 >  2:4232167/1‑101 (MQ=255)
                                                                cATTGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCg                                 >  1:4232166/1‑101 (MQ=255)
                                                                   tGCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGAtt                              >  1:4232171/1‑101 (MQ=255)
                                                                    gCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTc                             <  2:4231760/101‑1 (MQ=255)
                                                                    gCACAAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTc                             <  1:4231887/101‑1 (MQ=255)
                                                                      acaAAAAAACAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCAttt                                         <  1:4231827/87‑1 (MQ=255)
                                                                        aaaaaaaCAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGtt                         >  2:4232174/1‑101 (MQ=255)
                                                                        aaaaaaaCAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGtt                         >  1:4232175/1‑101 (MQ=255)
                                                                         aaaaaaCAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGttt                        <  2:4231872/101‑1 (MQ=255)
                                                                         aaaaaaCAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGttt                        <  1:4231860/101‑1 (MQ=255)
                                                                          aaaaaCAATCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTtta                       <  2:4231654/101‑1 (MQ=255)
                                                                                aaTCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCg                 >  2:4232179/1‑101 (MQ=255)
                                                                                 aTCAGCAGCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgc                >  1:4232180/1‑101 (MQ=255)
                                                                                   cagcagCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctg              <  1:4231643/101‑1 (MQ=255)
                                                                                   cagcagCCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctg              <  1:4231901/101‑1 (MQ=255)
                                                                                    agcagcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctgg             >  2:4232185/1‑101 (MQ=255)
                                                                                    agcagcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctgg             >  1:4232183/1‑101 (MQ=255)
                                                                                    agcagcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctgg             <  2:4231975/101‑1 (MQ=255)
                                                                                    agcagcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctgg             <  2:4231909/101‑1 (MQ=255)
                                                                                     gcagcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctggc            >  1:4232188/1‑101 (MQ=255)
                                                                                      cagcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctggct           >  1:4232190/1‑101 (MQ=255)
                                                                                       agcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctggctg          <  2:4231943/101‑1 (MQ=255)
                                                                                       agcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctggctg          <  1:4231755/101‑1 (MQ=255)
                                                                                       agcCTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCgctggctg          >  1:4232194/1‑101 (MQ=255)
                                                                                         ccTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTggg        <  1:4231869/101‑1 (MQ=255)
                                                                                          cTGAGTGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTgggg       >  1:4232197/1‑101 (MQ=255)
                                                                                               tGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTGGGGaaaaa  <  2:4231940/101‑1 (MQ=255)
                                                                                               tGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTGGGGaaaaa  <  2:4231864/101‑1 (MQ=255)
                                                                                               tGGTGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTGGGGaaaaa  >  1:4232200/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TAACAATCATCATTCACCACGTTTATGATTCTCTCCATCGACAGCAACGACGCTAATACCGCGCCATTGCACAAAAAAACAATCAGCAGCCTGAGTGGCGCAGTGGAGAGTATGATGAAAAAATTAGAAGATGTTGGTGTACTGGTAGCGCGCATTTTAATGCCGATTCTGTTTATTACCGCTGGCTGGGGAAAAA  >  NC_000913/3250448‑3250643

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