Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,269,012 A→G S200G (AGT→GGT)  yhaC → pentapetide repeats‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,0120AG98.8% 297.7 / NA 81S200G (AGT→GGT) yhaCpentapetide repeats‑related protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (42/38);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGA  >  NC_000913/3268879‑3269142
                                                                                                                                     |                                                                                                                                  
tctGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCa                                                                                                                               >  1:1257306/1‑139 (MQ=255)
  tGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAac                                                                                                                             >  2:420659/1‑139 (MQ=255)
     aCTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAAcaca                                                                                                                          >  2:564000/1‑139 (MQ=255)
          cGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACg                                                                                                                     >  2:541668/1‑139 (MQ=255)
            tAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGcc                                                                                                                   >  1:481841/1‑139 (MQ=255)
               tGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  1:420659/139‑1 (MQ=255)
                 tCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAAt                                                                                                              <  2:1257306/139‑1 (MQ=255)
                   gATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcac                                                                                                                 >  2:54595/1‑134 (MQ=255)
                   gATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcac                                                                                                                 <  1:54595/134‑1 (MQ=255)
                   gATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAAtct                                                                                                            >  1:209421/1‑139 (MQ=255)
                     tttAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTcc                                                                                                          <  2:168238/139‑1 (MQ=255)
                         gATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                >  2:1102784/1‑129 (MQ=255)
                         gATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  1:1102784/129‑1 (MQ=255)
                              ttatGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                 <  1:541668/139‑1 (MQ=255)
                                   aTTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACt                                                                                            >  2:371449/1‑139 (MQ=255)
                                                  atatGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGt                                                                             <  1:65538/139‑1 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                >  1:608883/1‑100 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  2:608883/100‑1 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCca                                                                                                                  >  2:499356/1‑98 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCca                                                                                                                  <  1:499356/98‑1 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                         >  2:235846/1‑139 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                         >  1:290170/1‑139 (MQ=255)
                                                      ggCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                         >  2:288868/1‑139 (MQ=255)
                                                       gCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                 >  2:1315666/1‑114 (MQ=255)
                                                       gCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattat                                                                                                 <  1:1315666/114‑1 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTACTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattatt                                                                                                <  2:928832/111‑1 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGattatt                                                                                                >  1:928832/1‑111 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAAc                                                                                   <  2:1486/124‑1 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAAc                                                                                   >  1:1486/1‑124 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaa                                                                        <  1:534029/135‑1 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaa                                                                        >  2:534029/1‑135 (MQ=255)
                                                           aTACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTCAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAg                                                                    >  1:1097044/1‑139 (MQ=255)
                                                            tACTATTTTAAATCCTTTTAAGTTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGACACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGc                                                                   >  2:790954/1‑139 (MQ=255)
                                                            tACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                >  2:159925/1‑94 (MQ=255)
                                                            tACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  1:159925/94‑1 (MQ=255)
                                                            tACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTaa                                                                                          >  1:215904/1‑116 (MQ=255)
                                                            tACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTaa                                                                                          <  2:215904/116‑1 (MQ=255)
                                                            tACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGc                                                                   <  2:481841/139‑1 (MQ=255)
                                                                aTTTTTAATCCTTGTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCATGTTAATTTTGTTAGTTCAACACAAATCGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTGAATAAGCATTa                                                               <  2:255158/139‑1 (MQ=255)
                                                                 ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                >  1:1029969/1‑89 (MQ=255)
                                                                 ttttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  2:1029969/89‑1 (MQ=255)
                                                                  tttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGc                                                             <  2:133313/139‑1 (MQ=255)
                                                                  tttAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGc                                                             <  2:290170/139‑1 (MQ=255)
                                                                     aaaTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTCGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGat                                                          <  1:790954/139‑1 (MQ=255)
                                                                      aaTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGata                                                         >  1:487785/1‑139 (MQ=255)
                                                                       aTCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGatat                                                        >  2:910473/1‑139 (MQ=255)
                                                                          cTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTa                                                     <  1:371449/139‑1 (MQ=255)
                                                                                   tGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTgaaga                                            <  1:910473/139‑1 (MQ=255)
                                                                                    gAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCa                                                                  >  2:1037800/1‑116 (MQ=255)
                                                                                    gAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCa                                                                  <  1:1037800/116‑1 (MQ=255)
                                                                                       ttCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  2:497466/67‑1 (MQ=255)
                                                                                       ttCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                >  1:497466/1‑67 (MQ=255)
                                                                                          aacaaATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATc                                     >  2:1214412/1‑139 (MQ=255)
                                                                                                   tACTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTc                            <  1:564000/139‑1 (MQ=255)
                                                                                                    aCTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                <  2:898477/54‑1 (MQ=255)
                                                                                                    aCTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCcaca                                                                                                                >  1:898477/1‑54 (MQ=255)
                                                                                                    aCTAATGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCg                           >  1:1068147/1‑139 (MQ=255)
                                                                                                        aTGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAg                       <  2:487785/139‑1 (MQ=255)
                                                                                                         tGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                         <  2:683255/88‑1 (MQ=255)
                                                                                                         tGCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTta                                                                         >  1:683255/1‑88 (MQ=255)
                                                                                                          gCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGc                                                                                 >  2:1243396/1‑79 (MQ=255)
                                                                                                          gCCAACATGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGc                                                                                 <  1:1243396/79‑1 (MQ=255)
                                                                                                                 tGTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTaaa              <  1:1214412/139‑1 (MQ=255)
                                                                                                                  gTTTCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAAt             <  2:1068147/139‑1 (MQ=255)
                                                                                                                   tttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATg            >  1:634877/1‑139 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCAt                                                                 >  2:1033432/1‑85 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCAt                                                                 <  1:1033432/85‑1 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGc                                                             >  1:1166035/1‑89 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGc                                                             <  2:1166035/89‑1 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTaaa              <  2:472242/136‑1 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTaaa              >  1:472242/1‑136 (MQ=255)
                                                                                                                    ttCAGGTTAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGa           >  1:1166125/1‑139 (MQ=255)
                                                                                                                     tCAGGTTAATTTTGTTGGTTCAACACAAAATGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGatat                                                        <  1:782154/93‑1 (MQ=255)
                                                                                                                     tCAGGTTAATTTTGTTGGTTCAACACAAAACGTCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGatat                                                        >  2:782154/1‑93 (MQ=255)
                                                                                                                          ttAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTa                                                               >  2:553237/1‑81 (MQ=255)
                                                                                                                          ttAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTa                                                               <  1:553237/81‑1 (MQ=255)
                                                                                                                           tAATTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGga    <  1:288868/139‑1 (MQ=255)
                                                                                                                             aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaat                                                                       <  1:309621/70‑1 (MQ=255)
                                                                                                                             aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTtaat                                                                       >  2:309621/1‑70 (MQ=255)
                                                                                                                             aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGgaga  >  1:151609/1‑139 (MQ=255)
                                                                                                                             aTTTTGTTGGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGgaga  >  2:192690/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                  
TCTGGACTTACGTAATGTCGATTTAGATTATTATGATTTCACAGATAAACATATGGCTAATACTATTTTAAATCCTTTTAAATTGAATTCAACAAATTTTACTAATGCCAACATGTTTCAGGTTAATTTTGTTAGTTCAACACAAAACGCCACAATCTCCTGGGATTATTTACTAAAAATAACGCCTGTTTTAATAAGCATTAGCGATATGTATTCTGAAGAAAAAATCAAGTTTGTCGAAAGTTGTTTAAATGAGCCTGGAGA  >  NC_000913/3268879‑3269142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑