Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,516,463 C→A G72V (GGC→GTC)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,516,4630CA98.8% 291.0 / ‑3.2 86G72V (GGC→GTC) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (37/48);  total (37/49)
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.76e-01

GACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACG  >  NC_000913/4516324‑4516595
                                                                                                                                           |                                                                                                                                    
gACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGc                                                                                                                                      >  1:618480/1‑140 (MQ=255)
     cGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTg                                                                                                                                 >  1:329597/1‑140 (MQ=255)
     cGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTg                                                                                                                                 >  1:975117/1‑140 (MQ=255)
            ggCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCt                                                                                                                          >  1:1556620/1‑140 (MQ=255)
            ggCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCt                                                                                                                          >  1:1399805/1‑140 (MQ=255)
            ggCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCt                                                                                                                          <  2:1399805/140‑1 (MQ=255)
                tcttctTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGcc                                                                                                                      >  1:1483780/1‑140 (MQ=255)
                 cttcttTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcg                                                                                                                     <  1:242283/140‑1 (MQ=255)
                 cttcttTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcg                                                                                                                     <  1:582383/140‑1 (MQ=255)
                  ttcttTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcgc                                                                                                                    >  2:1667106/1‑140 (MQ=255)
                  ttcttTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcgc                                                                                                                    >  2:1015784/1‑140 (MQ=255)
                   tcttTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcgcg                                                                                                                   <  1:290161/140‑1 (MQ=255)
                      tttGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGt                                                                                                                >  2:1055495/1‑140 (MQ=255)
                          gTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAgg                                                                                                            <  1:46823/140‑1 (MQ=255)
                              gggcgcgggcTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTgg                                                                                                        >  2:1404441/1‑140 (MQ=255)
                               ggcgcgggcTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGTTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGc                                                                                                       >  1:868674/1‑140 (MQ=255)
                               ggcgcgggcTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcg                                                                                                                     >  1:1169015/1‑126 (MQ=255)
                               ggcgcgggcTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCcg                                                                                                                     <  2:1169015/126‑1 (MQ=255)
                                     ggcTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAAc                                                                                                 >  2:226044/1‑140 (MQ=255)
                                      gcTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAAcc                                                                                                <  2:941669/140‑1 (MQ=255)
                                       cTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCg                                                                                               <  1:1049013/140‑1 (MQ=255)
                                        tCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCga                                                                                              <  1:1015784/140‑1 (MQ=255)
                                          cAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCgaga                                                                                            >  1:752912/1‑140 (MQ=255)
                                             cGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAggg                                                                                         <  1:467683/140‑1 (MQ=255)
                                                ccAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTaa                                                                                      >  2:28595/1‑140 (MQ=255)
                                                  aGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCga                                                                                              >  2:1304366/1‑130 (MQ=255)
                                                  aGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCga                                                                                              <  1:1304366/130‑1 (MQ=255)
                                                  aGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAAt                                                                                    >  2:238995/1‑140 (MQ=255)
                                                  aGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTATCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAAt                                                                                    >  2:1323804/1‑140 (MQ=255)
                                                     tGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCg                                                                                 <  2:957310/140‑1 (MQ=255)
                                                       ttATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCt                                                                               <  1:28595/140‑1 (MQ=255)
                                                         aTTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgt                                                                             >  2:1302054/1‑140 (MQ=255)
                                                          tttCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtg                                                                            <  2:1361617/140‑1 (MQ=255)
                                                           ttCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgt                                                                           <  2:618480/140‑1 (MQ=255)
                                                           ttCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgt                                                                           <  1:1670309/140‑1 (MQ=255)
                                                            tCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgtg                                                                          <  2:752912/140‑1 (MQ=255)
                                                            tCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgtg                                                                          <  1:226044/140‑1 (MQ=255)
                                                            tCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCtgtgtg                                                                          >  2:485453/1‑140 (MQ=255)
                                                               cAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCgg                                                                       <  1:1302054/140‑1 (MQ=255)
                                                               cAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCgg                                                                       <  1:1404441/140‑1 (MQ=255)
                                                               cAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCgg                                                                       <  1:106768/140‑1 (MQ=255)
                                                               cAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCgg                                                                       <  2:295982/140‑1 (MQ=255)
                                                                aGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGc                                                                      >  2:91904/1‑140 (MQ=255)
                                                                 gCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCa                                                                     >  2:270509/1‑140 (MQ=255)
                                                                 gCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCa                                                                     >  2:1673488/1‑140 (MQ=255)
                                                                 gCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCa                                                                     <  1:1313462/140‑1 (MQ=255)
                                                                  cGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCAt                                                                    >  2:344879/1‑140 (MQ=255)
                                                                    gTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATAc                                                                  <  1:485453/140‑1 (MQ=255)
                                                                     ttttACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACt                                                                 <  1:270509/140‑1 (MQ=255)
                                                                           cTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGa                                                           <  1:344879/140‑1 (MQ=255)
                                                                              cAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAgcg                                                        <  2:975117/140‑1 (MQ=255)
                                                                                 tCCGCTGCCGTCCAGCCGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCtt                                                     >  2:955911/1‑140 (MQ=255)
                                                                                   cGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGAtt                                                             >  2:514739/1‑130 (MQ=255)
                                                                                   cGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGAtt                                                             <  1:514739/130‑1 (MQ=255)
                                                                                   cGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTa                                                   >  2:218248/1‑140 (MQ=255)
                                                                                      tGCCGTCCAGCAGTTGTTGCAGGCCGCTATCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCg                                                <  1:1323804/140‑1 (MQ=255)
                                                                                        ccGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAg                                              <  1:218248/140‑1 (MQ=255)
                                                                                        ccGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAg                                              >  2:188475/1‑140 (MQ=255)
                                                                                              agcagTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATcc                                        <  2:868674/140‑1 (MQ=255)
                                                                                                  gttgttGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGt                                    >  2:690390/1‑140 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTTTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  2:753322/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  2:1579868/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  2:1483780/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  1:690390/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  2:329597/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  1:238995/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  1:1667106/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  1:714327/140‑1 (MQ=255)
                                                                                                   ttgttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  1:188475/140‑1 (MQ=255)
                                                                                                     gttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                 <  1:1055495/140‑1 (MQ=255)
                                                                                                     gttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                 >  2:1415943/1‑140 (MQ=255)
                                                                                                     gttgCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                 <  2:1556620/140‑1 (MQ=255)
                                                                                                             ccGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAAc                         <  1:1332193/140‑1 (MQ=255)
                                                                                                              cGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   >  1:243594/1‑129 (MQ=255)
                                                                                                              cGCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  2:243594/129‑1 (MQ=255)
                                                                                                               gCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   >  2:85437/1‑128 (MQ=255)
                                                                                                               gCTCTCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTg                                   <  1:85437/128‑1 (MQ=255)
                                                                                                                ctctCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                 >  2:1624470/1‑129 (MQ=255)
                                                                                                                ctctCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCg                                 <  1:1624470/129‑1 (MQ=255)
                                                                                                                ctctCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCtg                      <  2:736429/140‑1 (MQ=255)
                                                                                                                ctctCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCtg                      >  1:1205036/1‑140 (MQ=255)
                                                                                                                ctctCGACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCtg                      >  1:1088129/1‑140 (MQ=255)
                                                                                                                     gACGTCGTAATCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTgcag                 >  1:1139302/1‑140 (MQ=255)
                                                                                                                             aaTCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCg         >  1:396613/1‑140 (MQ=255)
                                                                                                                              aTCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCgg        <  1:1415943/140‑1 (MQ=255)
                                                                                                                               tCGCCGTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGa       <  1:91904/140‑1 (MQ=255)
                                                                                                                                    gTGCAGGACGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACg  >  1:688354/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                    
GACCACGGTCAGGGCATCTTCTTTTGGTGCGGGCGCGGGCTCCAGCGTCCAGCTGTTATTTCCCAGCGGTTTTACCTGCAGTCCGCTGCCGTCCAGCAGTTGTTGCAGGCCGCTCTCGACGTCGTAATCGCCGTGCAGGCCGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACG  >  NC_000913/4516324‑4516595

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: