Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,159,521 C→T A677A (GCG→GCA fhuE ← ferric coprogen/ferric rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,5210CT98.4% 220.6 / ‑4.4 64A677A (GCG→GCAfhuEferric coprogen/ferric rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/31);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATT  >  NC_000913/1159385‑1159653
                                                                                                                                        |                                                                                                                                    
tGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGtgcg                                                                                                                                   >  2:152259/1‑140 (MQ=255)
   tGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAg                                                                                                                                <  1:777710/140‑1 (MQ=255)
    gAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGc                                                                                                                               >  1:1216442/1‑140 (MQ=255)
    gAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGc                                                                                                                               >  1:776108/1‑140 (MQ=255)
                 gtgCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGgc                                                                                                                  <  1:87213/140‑1 (MQ=255)
                  tgCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCt                                                                                                                      <  1:777080/135‑1 (MQ=255)
                  tgCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCt                                                                                                                      >  2:777080/1‑135 (MQ=255)
                     cGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGATCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCgg                                                                                                              >  2:254902/1‑140 (MQ=255)
                      gTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGa                                                                                                             <  1:918631/140‑1 (MQ=255)
                        aGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAg                                                                                                           >  2:967258/1‑140 (MQ=255)
                        aGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAg                                                                                                           <  2:776108/140‑1 (MQ=255)
                          aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGt                                                                                                         >  2:165849/1‑140 (MQ=255)
                          aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGt                                                                                                         >  2:75639/1‑140 (MQ=255)
                          aCGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGt                                                                                                         >  1:1175298/1‑140 (MQ=255)
                            gATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGc                                                                                                       >  2:1231861/1‑140 (MQ=255)
                              tAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCa                                                                                                     >  1:1086664/1‑140 (MQ=255)
                               aGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGa                                                                                                             <  2:715596/131‑1 (MQ=255)
                               aGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGa                                                                                                             >  1:715596/1‑131 (MQ=255)
                                   ccTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACgg                                                                                                >  1:130819/1‑140 (MQ=255)
                                        cACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCa                                                                                           >  1:773988/1‑140 (MQ=255)
                                        cACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCa                                                                                           <  2:867978/140‑1 (MQ=255)
                                          cGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACg                                                                                         <  1:75639/140‑1 (MQ=255)
                                          cGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACg                                                                                         <  2:1175298/140‑1 (MQ=255)
                                          cGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACg                                                                                         <  2:482207/140‑1 (MQ=255)
                                           gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                                                                                        <  2:8078/140‑1 (MQ=255)
                                            ttGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgt                                                                                       >  2:574823/1‑140 (MQ=255)
                                              ggTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                                                                                     <  2:1216442/140‑1 (MQ=255)
                                              ggTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                                                                                     <  1:574823/140‑1 (MQ=255)
                                                tATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCg                                                                                   <  1:261956/140‑1 (MQ=255)
                                                     tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGt                                                                                             <  1:270122/125‑1 (MQ=255)
                                                     tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGt                                                                                             >  2:270122/1‑125 (MQ=255)
                                                       ggTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGTTCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATaca                                                                            >  2:1207146/1‑140 (MQ=255)
                                                       ggTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATaca                                                                            >  2:323042/1‑140 (MQ=255)
                                                       ggTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATaca                                                                            >  1:737482/1‑140 (MQ=255)
                                                                  aCAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGc                                                                 <  1:1231560/140‑1 (MQ=255)
                                                                       ttATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGtt                                                            <  1:186307/140‑1 (MQ=255)
                                                                           tGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAaca                                                        >  2:1098739/1‑140 (MQ=255)
                                                                                  ccccTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCAccgc                                                 >  1:1108050/1‑140 (MQ=255)
                                                                                    ccTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCAccgccg                                               <  2:130819/140‑1 (MQ=255)
                                                                                     cTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGa                                              <  1:967258/140‑1 (MQ=255)
                                                                                         aaCGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGTTCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGt                                          <  1:1207146/140‑1 (MQ=255)
                                                                                         aaCGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGt                                          <  2:1086664/140‑1 (MQ=255)
                                                                                         aaCGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGt                                          <  1:254902/140‑1 (MQ=255)
                                                                                             agaAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAAc                                      >  2:829769/1‑140 (MQ=255)
                                                                                             agaAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAAc                                      >  2:794479/1‑140 (MQ=255)
                                                                                               aaGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAActc                                    >  2:71287/1‑140 (MQ=255)
                                                                                                aGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACACCGCCGACCGTCAACTCTg                                  <  1:323042/140‑1 (MQ=255)
                                                                                                  tttttCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTgg                                 >  2:1131841/1‑140 (MQ=255)
                                                                                                  tttttCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTgg                                 >  1:242500/1‑140 (MQ=255)
                                                                                                  tttttCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCATCAACTCTgg                                 >  1:723647/1‑140 (MQ=255)
                                                                                                           aCCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAgg                        <  2:351657/140‑1 (MQ=255)
                                                                                                           aCCTGGTAGCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAgg                        <  2:587939/140‑1 (MQ=255)
                                                                                                                   gCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGAt                <  1:1098739/140‑1 (MQ=255)
                                                                                                                   gCGGGTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGAt                <  2:242500/140‑1 (MQ=255)
                                                                                                                       gTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCt            >  1:855764/1‑140 (MQ=255)
                                                                                                                       gTGAAGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGATGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCt            >  2:866554/1‑140 (MQ=255)
                                                                                                                          aaGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAg              <  2:54716/135‑1 (MQ=255)
                                                                                                                          aaGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAg              >  1:54716/1‑135 (MQ=255)
                                                                                                                          aaGAGATCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGt         >  2:193482/1‑140 (MQ=255)
                                                                                                                             agaTCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGaa      <  1:71287/140‑1 (MQ=255)
                                                                                                                             agaTCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGaa      <  2:31901/140‑1 (MQ=255)
                                                                                                                                tCCACCAGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACAt   <  2:855764/140‑1 (MQ=255)
                                                                                                                                 ccaccaGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACAtt  >  1:849390/1‑140 (MQ=255)
                                                                                                                                 ccaccaGTGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACAtt  >  1:876412/1‑140 (MQ=255)
                                                                                                                                 ccaccaGTGAGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACAtt  >  1:136121/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATT  >  NC_000913/1159385‑1159653

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