Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,191,379 A→G T250A (ACG→GCG)  yqiI → putative fimbrial protein YqiI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,191,3790AG98.5% 220.9 / NA 65T250A (ACG→GCG) yqiIputative fimbrial protein YqiI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/41);  minor base T (1/0);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTAATGGTCTGC  >  NC_000913/3191254‑3191505
                                                                                                                             |                                                                                                                              
ttCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTc                                                                                                                  >  2:89373/1‑140 (MQ=255)
                 aTTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCTTTTGGCGTAGTCAAGACCAAACATCCggg                                                                                                 >  1:241749/1‑140 (MQ=255)
                  ttCATGTCCTTGACTTCTTTGTTAATGTTAATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGt                                                                                                <  1:291739/140‑1 (MQ=255)
                  ttCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGt                                                                                                <  2:576178/140‑1 (MQ=255)
                  ttCATGTCCTTGACTGCTTTGTAAATGTTGATCTGGAACCAGCTGATGACGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGt                                                                                                <  2:307429/140‑1 (MQ=255)
                  ttCATGTCCTTGACTGCCTTGTTAATGTTGACCTGCAACCAGCTGATGCCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGt                                                                                                <  2:296266/140‑1 (MQ=255)
                   tCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTg                                                                                               >  2:140233/1‑140 (MQ=255)
                       gTCCTTGACTGCTTTTTTAATGTTGCTCTGGAACCAGTTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           <  1:259778/140‑1 (MQ=255)
                       gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           >  1:616772/1‑140 (MQ=255)
                       gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           <  2:117121/140‑1 (MQ=255)
                       gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           <  1:607357/140‑1 (MQ=255)
                       gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           <  2:474719/140‑1 (MQ=255)
                       gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           <  1:731853/140‑1 (MQ=255)
                       gTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           >  1:78265/1‑140 (MQ=255)
                       gTCCTTGACTGCTTTGTGAATGGTGATCTGGGACCTGCTTATGGCGTTTTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGc                                                                                           <  2:281281/140‑1 (MQ=255)
                          cTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTg                                                                                        >  1:21169/1‑140 (MQ=255)
                          cTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGATTAGCGTGGTCATGCCCAAAGCTCCGGGTGCGGCCTg                                                                                        >  1:220961/1‑140 (MQ=255)
                            tGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTa                                                                                      >  1:40914/1‑140 (MQ=255)
                             gACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTAc                                                                                     >  2:717259/1‑140 (MQ=255)
                                 gCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAg                                                                                 <  1:103356/140‑1 (MQ=255)
                                 gCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAg                                                                                 <  1:526937/140‑1 (MQ=255)
                                 gCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTCTACTGAg                                                                                 <  2:422888/140‑1 (MQ=255)
                                 gCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTCTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAg                                                                                 <  1:15637/140‑1 (MQ=255)
                                  cTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGTCGTTTAGCTTAGGCATTACCAAAGATCCGGGTGCGGCCTGTACTGAGc                                                                                >  2:247953/1‑140 (MQ=255)
                                    ttGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAg                                                                              >  2:27172/1‑140 (MQ=255)
                                        tAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTa                                                                          <  1:649342/140‑1 (MQ=255)
                                           tGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAata                                                                       <  1:27172/140‑1 (MQ=255)
                                            gTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAatat                                                                      >  2:452694/1‑140 (MQ=255)
                                                aTCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGGCTCTACTGAACAGTTTAATATTTTa                                                                  >  1:312202/1‑140 (MQ=255)
                                                 tCTGGAACCAGCTGATGTCGTTCTCGACTTTGGTACAATCAATTCCCGTACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAg                                                                 <  2:241749/140‑1 (MQ=255)
                                                 tCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAg                                                                 >  1:553831/1‑140 (MQ=255)
                                                   tGGAACCAGCTGATGGCGTTGTGGACTTTGGCACATTAACTTCCCGAATAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAggg                                                               <  2:312202/140‑1 (MQ=255)
                                                   tGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAggg                                                               <  2:40914/140‑1 (MQ=255)
                                                   tGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAggg                                                               <  2:92070/140‑1 (MQ=255)
                                                          aGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGttttt                                                        >  2:161476/1‑140 (MQ=255)
                                                           gCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGtttttt                                                       >  2:505425/1‑140 (MQ=255)
                                                              gATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTGAc                                                    <  2:41569/140‑1 (MQ=255)
                                                              gATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCAc                                                    <  2:738231/140‑1 (MQ=255)
                                                              gATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCAc                                                    <  2:42692/140‑1 (MQ=255)
                                                                     ttGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGat                                             >  1:148735/1‑140 (MQ=255)
                                                                     ttGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGat                                             <  1:717259/140‑1 (MQ=255)
                                                                      tGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGata                                            >  2:519622/1‑140 (MQ=255)
                                                                      tGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGata                                            >  2:498895/1‑140 (MQ=255)
                                                                            cTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCCGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGa                                      <  1:140233/140‑1 (MQ=255)
                                                                             tttGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAg                                     <  2:506176/140‑1 (MQ=255)
                                                                             tttGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTAATGAGCAGTCTAAAATTGTAGGGAGTGTTTTAACTACGGATATTTTGAg                                     >  2:231864/1‑140 (MQ=255)
                                                                              ttGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGt                                    <  1:498895/140‑1 (MQ=255)
                                                                                ggTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGa                                  >  2:40656/1‑140 (MQ=255)
                                                                                       taaatTCCCGAACAATTAAAAATACCAGCGTAAGTTAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGc                           <  2:262344/140‑1 (MQ=255)
                                                                                          atTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCAt                        <  1:40656/140‑1 (MQ=255)
                                                                                          atTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCAt                        <  1:394892/140‑1 (MQ=255)
                                                                                          atTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCAt                        <  2:553831/140‑1 (MQ=255)
                                                                                          atTCCCGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCAt                        <  1:89373/140‑1 (MQ=255)
                                                                                           ttCCCGAACAATTAATAATACCAGCGTGATTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATc                       <  1:303198/140‑1 (MQ=255)
                                                                                            tCCCGAACAATTAAATATACCATCGGGAGTGAGGCGTGTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCt                      <  1:247953/140‑1 (MQ=255)
                                                                                              ccGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATc                       >  1:481766/1‑137 (MQ=255)
                                                                                              ccGAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATc                       <  2:481766/137‑1 (MQ=255)
                                                                                                gAACAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGat                  <  2:475737/140‑1 (MQ=255)
                                                                                                  aCAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGatat                <  1:519622/140‑1 (MQ=255)
                                                                                                  aCAATTAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGatat                <  2:78265/140‑1 (MQ=255)
                                                                                                       tAAAAATACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTa           <  2:21169/140‑1 (MQ=255)
                                                                                                             tACCAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGAGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTAATGGTc     <  2:233968/140‑1 (MQ=255)
                                                                                                               ccAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTAATGGTCTg   <  1:161476/140‑1 (MQ=255)
                                                                                                                cAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTAATGGTCTGc  >  1:566129/1‑140 (MQ=255)
                                                                                                                cAGCGTGAGTGAGGCGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTAATGGTCTGc  >  1:457175/1‑140 (MQ=255)
                                                                                                                             |                                                                                                                              
TTCGCGGATTAAATAACATTCATGTCCTTGACTGCTTTGTTAATGTTGATCTGGAACCAGCTGATGGCGTTGTCGACTTTGGTAAAATAAATTCCCGAACAATTAAAAATACCAGCGTGAGTGAGACGTTTAGCGTAGTCATGACCAAAGATCCGGGTGCGGCCTGTACTGAGCAGTTTAATATTTTAGGGAGTTTTTTCACTACGGATATTTTGAGTGATTATAGCCATCTGGATATAGGTAATGGTCTGC  >  NC_000913/3191254‑3191505

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