Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,665,586 T→C A411A (GCA→GCG yhjA ← putative cytochrome C peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,665,5860TC98.6% 431.7 / ‑3.2 139A411A (GCA→GCGyhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (82/55);  total (82/57)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 144 total aligned reads displayed.

TCCTGCGGCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGTGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAATATC  >  NC_000913/3665497‑3665672
                                                                                         |                                                                                      
tcctGCGGCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGcc                                                                                     >  1:4466825/1‑93 (MQ=255)
   tGCGGCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGtt                                                                                  >  1:559140/1‑93 (MQ=255)
   tGCGGCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGtt                                                                                  >  1:2056054/1‑93 (MQ=255)
      ggCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTaaa                                                                               >  1:4130560/1‑93 (MQ=255)
      ggCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTaaa                                                                               <  1:2811022/93‑1 (MQ=255)
      ggCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTaaa                                                                               >  2:3900719/1‑93 (MQ=255)
      ggCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTaaa                                                                               <  2:4481831/93‑1 (MQ=255)
       gCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAg                                                                              >  2:4540893/1‑93 (MQ=255)
         aGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCa                                                                            >  1:1137170/1‑93 (MQ=255)
          gCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCaa                                                                           >  1:124182/1‑93 (MQ=255)
          gCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCaa                                                                           <  2:1076962/93‑1 (MQ=255)
           ctctTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAAc                                                                          >  1:3247977/1‑93 (MQ=255)
            tctTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACg                                                                         >  2:113910/1‑93 (MQ=255)
            tctTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACg                                                                         >  2:113911/1‑93 (MQ=255)
            tctTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACg                                                                         >  2:2897502/1‑93 (MQ=255)
             ctTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGt                                                                        >  2:1221391/1‑93 (MQ=255)
              tttGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGtt                                                                       >  2:4502300/1‑93 (MQ=255)
               ttGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTa                                                                      >  1:1715938/1‑93 (MQ=255)
               ttGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTa                                                                      <  2:1222065/93‑1 (MQ=255)
               ttGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTa                                                                      >  2:2087750/1‑93 (MQ=255)
                tGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTAc                                                                     >  2:4519186/1‑93 (MQ=255)
                tGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTAc                                                                     <  1:3896371/93‑1 (MQ=255)
                 gCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACg                                                                    >  2:637824/1‑93 (MQ=255)
                 gCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACg                                                                    >  2:3627262/1‑93 (MQ=255)
                    tACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGGAGTACGGCGCCGTTAAAGCCACGTTACGTaa                                                                 >  1:426176/1‑93 (MQ=255)
                     tccTGATAGCGCAGCATCAGTTTCCCCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTaaa                                                                <  1:1034788/92‑1 (MQ=255)
                      ccTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAAc                                                               <  1:2143796/93‑1 (MQ=255)
                       cTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAAcc                                                              >  1:4659848/1‑93 (MQ=255)
                       cTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAAcc                                                              >  2:2474964/1‑93 (MQ=255)
                        tGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCg                                                             <  1:3606239/93‑1 (MQ=255)
                          aTAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGTCACGTCACCGCGATGGAAGTACGGCGCAGTTTAAGCCACTCTACGTAAACCggg                                                           >  2:570741/1‑93 (MQ=255)
                             gcgcAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTAc                                                        <  2:3805613/93‑1 (MQ=255)
                             gcgcAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTAc                                                        >  2:521170/1‑93 (MQ=255)
                              cgcAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACt                                                       <  1:2267243/93‑1 (MQ=255)
                                  gcaTCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCt                                                   >  2:1587618/1‑93 (MQ=255)
                                  gcaTCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCt                                                   >  1:499764/1‑93 (MQ=255)
                                   caTCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTg                                                  <  1:250765/93‑1 (MQ=255)
                                    aTCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACCGCGCCGTTAAAGCAACGTTACGTAAACCGGGGACTTTCTgg                                                 >  2:3694184/1‑92 (MQ=255)
                                     tCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGAc                                                <  1:541070/93‑1 (MQ=255)
                                     tCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGAc                                                <  1:4286944/93‑1 (MQ=255)
                                        gTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCa                                             >  1:3101081/1‑93 (MQ=255)
                                        gTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCa                                             >  1:2892821/1‑93 (MQ=255)
                                            cACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAAttt                                         >  1:4704238/1‑93 (MQ=255)
                                            cACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAAttt                                         <  1:3920704/93‑1 (MQ=255)
                                            cACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAAttt                                         >  1:960286/1‑93 (MQ=255)
                                             aCCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTa                                        >  2:4634856/1‑93 (MQ=255)
                                             aCCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTa                                        <  1:214928/93‑1 (MQ=255)
                                               cGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGTGTACTTTCTGACGCAATTTATc                                      <  2:2990891/93‑1 (MQ=255)
                                               cGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATc                                      >  2:3519227/1‑93 (MQ=255)
                                                 ccccGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCAc                                    <  2:1395792/93‑1 (MQ=255)
                                                 ccccGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCAc                                    <  1:1485834/93‑1 (MQ=255)
                                                 ccccGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCAc                                    <  2:4118734/93‑1 (MQ=255)
                                                   ccGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGc                                  <  1:1056924/93‑1 (MQ=255)
                                                   ccGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGc                                  >  2:2056067/1‑93 (MQ=255)
                                                   ccGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGc                                  <  2:532273/93‑1 (MQ=255)
                                                     gTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGctc                                <  1:2560570/93‑1 (MQ=255)
                                                      tCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGctct                               >  1:793946/1‑93 (MQ=255)
                                                      tCCAGAGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCTATTTATCACGctct                               >  2:1700718/1‑93 (MQ=255)
                                                       ccAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCtctt                              >  2:3428312/1‑93 (MQ=255)
                                                       ccAGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCtctt                              <  1:3682019/93‑1 (MQ=255)
                                                       ccAGCGTCGGAACGTAGCCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTGACCCGGGGCCTTTCTGACGCAATTTATCACGAtctt                              >  2:3428311/1‑93 (MQ=255)
                                                         aGCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCtcttct                            >  1:1338737/1‑93 (MQ=255)
                                                          gCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTcttctt                           >  1:2724163/1‑93 (MQ=255)
                                                          gCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTcttctt                           >  1:320156/1‑93 (MQ=255)
                                                          gCGTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTcttctt                           >  2:1897293/1‑93 (MQ=255)
                                                            gTCGGCACGTCACCGCGATGGAAGTACGGTGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTa                         <  1:740580/93‑1 (MQ=255)
                                                            gTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTa                         <  1:3280479/93‑1 (MQ=255)
                                                            gTCGGCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTa                         >  1:2653904/1‑93 (MQ=255)
                                                             tCGGCACGTCACCGCAATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAg                        <  1:322231/93‑1 (MQ=255)
                                                               ggCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTc                      <  1:1697871/93‑1 (MQ=255)
                                                               ggCACGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTc                      >  2:425188/1‑93 (MQ=255)
                                                                  aCGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTcaca                   <  2:1667282/93‑1 (MQ=255)
                                                                  aCGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTcaca                   <  2:2709723/93‑1 (MQ=255)
                                                                  aCGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTcaca                   >  1:716781/1‑93 (MQ=255)
                                                                  aCGTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGCTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTcaca                   >  2:2477842/1‑93 (MQ=255)
                                                                    gTCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAACCCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACAtt                 <  1:4133053/93‑1 (MQ=255)
                                                                     tCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTc                >  1:4479155/1‑93 (MQ=255)
                                                                     tCACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTc                <  1:1700480/93‑1 (MQ=255)
                                                                      cACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCa               <  1:4430497/93‑1 (MQ=255)
                                                                      cACCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCa               <  2:1076935/93‑1 (MQ=255)
                                                                       aCCGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCAt              >  1:4021041/1‑93 (MQ=255)
                                                                        ccGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCAtc             >  1:3680233/1‑92 (MQ=255)
                                                                        ccGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATa             <  1:4041124/93‑1 (MQ=255)
                                                                        ccGCGATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATa             <  2:4105516/93‑1 (MQ=255)
                                                                         cgcgATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATAc            >  1:963749/1‑93 (MQ=255)
                                                                          gcgATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACg           >  1:3868223/1‑93 (MQ=255)
                                                                            gATGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGAc         <  1:4571781/93‑1 (MQ=255)
                                                                             aTGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGAcc        >  1:4556269/1‑93 (MQ=255)
                                                                              tGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCa       >  1:868066/1‑93 (MQ=255)
                                                                              tGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCa       >  1:613728/1‑93 (MQ=255)
                                                                              tGGAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTAGTCAACTTCATAAGACCaa      >  1:98476/1‑93 (MQ=255)
                                                                               ggAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCaa      <  2:4205458/93‑1 (MQ=255)
                                                                               ggAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCaa      >  2:2555273/1‑93 (MQ=255)
                                                                                gAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAat     <  1:1924920/93‑1 (MQ=255)
                                                                                gAAGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAat     <  1:4657634/93‑1 (MQ=255)
                                                                                 aaGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAata    >  1:3272943/1‑93 (MQ=255)
                                                                                  aGTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAatat   <  1:364420/93‑1 (MQ=255)
                                                                                   gTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAATATc  >  1:1556933/1‑93 (MQ=255)
                                                                                   gTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAATATc  >  2:2832870/1‑93 (MQ=255)
                                                                                   gTACGGCGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAATATc  >  1:2745563/1‑93 (MQ=255)
                                                                                         |                                                                                      
TCCTGCGGCAGCTCTTTGCCTACCTGATAGCGCAGCATCAGTTTCACCGCCCCGTCCAGCGTCGGCACGTCACCGCGATGGAAGTACGGTGCCGTTAAAGCAACGTTACGTAAACCGGGTACTTTCTGACGCAATTTATCACGCTCTTCTTTAGTCACATTCATACGACCAATATC  >  NC_000913/3665497‑3665672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: