Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,511,424 A→G Y12H (TAC→CAC)  fecE ← ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,4240AG98.6% 256.3 / ‑4.4 73Y12H (TAC→CAC) fecEferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/35);  total (38/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 = 9.07e-01

CCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGC  >  NC_000913/4511288‑4511556
                                                                                                                                        |                                                                                                                                    
ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACt                                                                                                                                   >  1:467670/1‑140 (MQ=255)
ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACt                                                                                                                                   >  2:716382/1‑140 (MQ=255)
ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACt                                                                                                                                   >  2:935145/1‑140 (MQ=255)
ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACt                                                                                                                                   >  2:1592868/1‑140 (MQ=255)
ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACt                                                                                                                                   >  1:330676/1‑140 (MQ=255)
      tGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGt                                                                                                                             <  1:935145/140‑1 (MQ=255)
      tGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGAGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGt                                                                                                                             <  2:51613/140‑1 (MQ=255)
        cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCa                                                                                                                           >  1:1074981/1‑140 (MQ=255)
        cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCa                                                                                                                           >  2:470479/1‑140 (MQ=255)
        cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCa                                                                                                                           >  2:183118/1‑140 (MQ=255)
        cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGATAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCa                                                                                                                           >  2:719610/1‑140 (MQ=255)
        cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGATAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCa                                                                                                                           >  1:179082/1‑140 (MQ=255)
           cATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGAc                                                                                                                                    >  1:1170119/1‑128 (MQ=255)
             ttAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGAttt                                                                                                                      >  2:802920/1‑140 (MQ=255)
             ttAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGAttt                                                                                                                      >  2:1163729/1‑140 (MQ=255)
                       cgAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTa                                                                                                            <  1:1523403/140‑1 (MQ=255)
                       cgAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTa                                                                                                            <  2:137156/140‑1 (MQ=255)
                          aaaaCAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAg                                                                                                         <  2:888122/140‑1 (MQ=255)
                              cAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCAt                                                                                                     <  1:1348322/140‑1 (MQ=255)
                                  ttAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAt                                                                                                 >  2:1021680/1‑140 (MQ=255)
                                  ttAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTAt                                                                                                 <  2:467670/140‑1 (MQ=255)
                                  ttAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGCAAAGTCATTTAt                                                                                                 <  1:432771/140‑1 (MQ=255)
                                   tAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATc                                                                                                <  2:1188634/140‑1 (MQ=255)
                                   tAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATc                                                                                                >  2:1546037/1‑140 (MQ=255)
                                     aCAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGc                                                                                              <  2:1096684/140‑1 (MQ=255)
                                         cGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCa                                                                                             >  1:214144/1‑137 (MQ=255)
                                         cGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCa                                                                                             <  2:214144/137‑1 (MQ=255)
                                              aTTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTcaca                                                                                     <  2:580787/140‑1 (MQ=255)
                                                ttCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAg                                                                                   <  1:1021680/140‑1 (MQ=255)
                                                 tCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGc                                                                                  <  2:481210/140‑1 (MQ=255)
                                                        aCCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGa                                                                           >  1:86369/1‑140 (MQ=255)
                                                        aCCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTCGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGa                                                                           <  2:82174/140‑1 (MQ=255)
                                                         ccGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGAc                                                                          >  1:288174/1‑140 (MQ=255)
                                                                 aCCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAAc                                                                               >  1:1451493/1‑127 (MQ=255)
                                                                 aCCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAAc                                                                               <  2:1451493/127‑1 (MQ=255)
                                                                   cGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGgc                                                                >  1:955867/1‑140 (MQ=255)
                                                                    gATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGgcg                                                               <  1:740386/140‑1 (MQ=255)
                                                                     aTCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGgcgc                                                              >  1:224046/1‑140 (MQ=255)
                                                                     aTCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGgcgc                                                              >  1:1416974/1‑140 (MQ=255)
                                                                      tCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCa                                                             >  2:392593/1‑140 (MQ=255)
                                                                          ggCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGa                                                         >  1:738973/1‑140 (MQ=255)
                                                                               tGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATg                                                    >  2:955692/1‑140 (MQ=255)
                                                                               tGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATg                                                    >  2:15381/1‑140 (MQ=255)
                                                                                  tCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCg                                                 <  1:1027363/140‑1 (MQ=255)
                                                                                      ccccGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCa                                             >  2:296464/1‑140 (MQ=255)
                                                                                        ccGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGc                                           >  2:104052/1‑140 (MQ=255)
                                                                                         cGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCa                                          >  1:145549/1‑140 (MQ=255)
                                                                                         cGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCa                                          <  1:179520/140‑1 (MQ=255)
                                                                                         cGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCa                                          <  2:738973/140‑1 (MQ=255)
                                                                                           ttGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCAcc                                        >  2:1193208/1‑139 (MQ=255)
                                                                                           ttGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACg                                        >  1:1403616/1‑140 (MQ=255)
                                                                                               cAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCaa                                    <  2:179082/140‑1 (MQ=255)
                                                                                                 gTGAGAGTGAAACGTCGATAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAAcc                                  <  1:116086/140‑1 (MQ=255)
                                                                                                    agagTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCggg                               <  2:140575/140‑1 (MQ=255)
                                                                                                    agagTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCggg                               <  1:470479/140‑1 (MQ=255)
                                                                                                    agagTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCggg                               <  1:716382/140‑1 (MQ=255)
                                                                                                      agTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAg                             <  2:330676/140‑1 (MQ=255)
                                                                                                      agTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAg                             >  2:1281608/1‑140 (MQ=255)
                                                                                                      agTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAg                             <  2:827396/140‑1 (MQ=255)
                                                                                                        tGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCt                           >  2:83797/1‑140 (MQ=255)
                                                                                                        tGAAACGTCGTTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAACGGCGGTCAGCACGCCAACCGGGAGCt                           <  2:288174/140‑1 (MQ=255)
                                                                                                              gtcgtTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTg                     >  2:1243968/1‑140 (MQ=255)
                                                                                                               tcgtTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTgg                    <  1:296464/140‑1 (MQ=255)
                                                                                                               tcgtTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTgg                    <  2:1402718/140‑1 (MQ=255)
                                                                                                                cgtTAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggg                   >  2:845219/1‑140 (MQ=255)
                                                                                                                  ttAAGTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                 <  2:224046/140‑1 (MQ=255)
                                                                                                                      gTACCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGa             <  1:83797/140‑1 (MQ=255)
                                                                                                                        aCCTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGaat           >  1:778939/1‑140 (MQ=255)
                                                                                                                          cTTGTCTGTCCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGaataa         >  1:894752/1‑140 (MQ=255)
                                                                                                                              tctgtcCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtct     <  2:1416974/140‑1 (MQ=255)
                                                                                                                                tgtcCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                 <  1:1597799/126‑1 (MQ=255)
                                                                                                                                tgtcCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTggggg                 >  2:1597799/1‑126 (MQ=255)
                                                                                                                                tgtcCCGTGACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCg   <  2:1244649/140‑1 (MQ=255)
                                                                                                                                 gtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTTTCGc  >  1:378221/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGC  >  NC_000913/4511288‑4511556

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: