Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,084,051 T→C N591S (AAC→AGC)  fdoG ← formate dehydrogenase O subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,084,0510TC98.4% 222.8 / ‑4.4 64N591S (AAC→AGC) fdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/32);  total (31/33)
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.52e-01

TGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCTTGCCCT  >  NC_000913/4083916‑4084186
                                                                                                                                       |                                                                                                                                       
tGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTaga                                                                                                                                     >  2:670365/1‑140 (MQ=255)
     cTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAgg                                                                                                                                   <  1:678402/137‑1 (MQ=255)
     cTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAgg                                                                                                                                   >  2:678402/1‑137 (MQ=255)
             aGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTttcagtgttc                                                                                                                        >  1:1065738/1‑140 (MQ=255)
                    tGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGt                                                                                                                 <  1:5678/140‑1 (MQ=255)
                    tGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGt                                                                                                                 <  1:348003/140‑1 (MQ=255)
                     gAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTc                                                                                                                >  1:452971/1‑140 (MQ=255)
                      aGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCg                                                                                                               <  2:837584/140‑1 (MQ=255)
                       gTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGa                                                                                                              <  1:378667/140‑1 (MQ=255)
                            cGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTg                                                                                                         >  2:226077/1‑140 (MQ=255)
                            cGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTg                                                                                                         >  1:69868/1‑140 (MQ=255)
                                     cGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGt                                                                                                                 >  2:325801/1‑123 (MQ=255)
                                     cGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGt                                                                                                                 <  1:325801/123‑1 (MQ=255)
                                             tCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGttt                                                                                        <  1:877005/140‑1 (MQ=255)
                                              ccGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTc                                                                                       >  1:367727/1‑140 (MQ=255)
                                                gcgAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCga                                                                                     <  1:926374/140‑1 (MQ=255)
                                                 cgAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCgac                                                                                    >  2:536615/1‑140 (MQ=255)
                                                 cgAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCgac                                                                                    >  2:23907/1‑140 (MQ=255)
                                                 cgAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCgac                                                                                    >  2:587514/1‑140 (MQ=255)
                                                   aaGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCgacag                                                                                  >  2:1100665/1‑140 (MQ=255)
                                                         gTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAAcc                                                                            >  1:881879/1‑140 (MQ=255)
                                                         gTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAAcc                                                                            <  2:881879/140‑1 (MQ=255)
                                                         gTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAAcc                                                                            >  1:1119711/1‑140 (MQ=255)
                                                           cGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGa                                                                          <  1:1100665/140‑1 (MQ=255)
                                                           cGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGa                                                                          >  2:773746/1‑140 (MQ=255)
                                                            gATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGAt                                                                         <  1:773746/140‑1 (MQ=255)
                                                              tGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCa                                                                       >  2:157861/1‑140 (MQ=255)
                                                               ggCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCAc                                                                      <  1:226077/140‑1 (MQ=255)
                                                                    aCGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgt                                                                 >  2:355866/1‑140 (MQ=255)
                                                                    aCGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgt                                                                 >  1:247071/1‑140 (MQ=255)
                                                                    aCGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgt                                                                 >  1:282892/1‑140 (MQ=255)
                                                                       gAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgtttt                                                              <  2:435267/140‑1 (MQ=255)
                                                                       gAACACTTCGGTCTGGATCTTCAACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgtttt                                                              >  2:878079/1‑140 (MQ=255)
                                                                        aaCACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgttttt                                                             >  2:471509/1‑140 (MQ=255)
                                                                           acTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACtttgtttttgtt                                                          <  1:893880/140‑1 (MQ=255)
                                                                             ttCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCg                                                        >  1:804652/1‑140 (MQ=255)
                                                                               cGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCggg                                                      <  2:1110045/140‑1 (MQ=255)
                                                                                ggTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGa                                                     <  1:157861/140‑1 (MQ=255)
                                                                                 gTCTGGATCTTCAACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGaa                                                    <  1:878079/140‑1 (MQ=255)
                                                                                   cTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATg                                                  <  2:282892/140‑1 (MQ=255)
                                                                                   cTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATg                                                  <  2:247071/140‑1 (MQ=255)
                                                                                      gATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAgg                                               >  1:23939/1‑140 (MQ=255)
                                                                                         cTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCaa                                            <  2:1065738/140‑1 (MQ=255)
                                                                                         cTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCaa                                            <  1:551762/140‑1 (MQ=255)
                                                                                           tCGATGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACa                                          <  2:1117973/140‑1 (MQ=255)
                                                                                           tCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACa                                          <  1:23907/140‑1 (MQ=255)
                                                                                                 aGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTa                                    <  1:670365/140‑1 (MQ=255)
                                                                                                     aaCTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGc                                <  2:139781/140‑1 (MQ=255)
                                                                                                      aCTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGcc                               <  1:782296/140‑1 (MQ=255)
                                                                                                          cGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTgg                           <  1:355866/140‑1 (MQ=255)
                                                                                                          cGTTCAGCTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTgg                           <  1:536615/140‑1 (MQ=255)
                                                                                                                 cTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTgg                           >  2:159578/1‑133 (MQ=255)
                                                                                                                 cTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTgg                           <  1:159578/133‑1 (MQ=255)
                                                                                                                 cTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGatat                    <  1:471509/140‑1 (MQ=255)
                                                                                                                 cTCACCGTGGTTCTGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGatat                    >  2:270473/1‑140 (MQ=255)
                                                                                                                              tGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCtt       >  2:1000218/1‑140 (MQ=255)
                                                                                                                              tGCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCtt       >  1:905193/1‑140 (MQ=255)
                                                                                                                               gCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCTTg      <  1:1000218/140‑1 (MQ=255)
                                                                                                                               gCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCTTg      >  1:603957/1‑140 (MQ=255)
                                                                                                                               gCCAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGGTAAAGCCCTGGCAGATATAGCCATTGACCTTg      >  1:47753/1‑140 (MQ=255)
                                                                                                                                ccAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTg            >  1:1169528/1‑133 (MQ=255)
                                                                                                                                ccAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTg            <  2:1169528/133‑1 (MQ=255)
                                                                                                                                ccAGAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCTTGc     <  1:475187/140‑1 (MQ=255)
                                                                                                                                   gAAGCTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCTTGCCct  >  1:886501/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGCCACTGTAACCAGCGACCTGAGTTAACGATAGAACCGTTCTCTTCCGCGAAGCAGGTCGATGGCAGACGGAACACTTCGGTCTGGATCTTCGACGAGTCAACTTCGTTCAGCTCACCGTGGTTCTGCCAGAAGTTAGAGGTTTCAGTGTTCAGTGGGTCGATGGTGACGAGGAACTTCAGTTTCGACAGACAACCGATCACTTTGTTTTTGTTCGGGAATGAGGCAACAGGGTTAAAGCCCTGGCAGATATAGCCATTGACCTTGCCCT  >  NC_000913/4083916‑4084186

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: