Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,137,180 C→G A160P (GCA→CCA)  yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,1800CG96.9% 352.0 / 0.9 98A160P (GCA→CCA) yijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base G (42/53);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCA  >  NC_000913/4137083‑4137276
                                                                                                 |                                                                                                
aaCAGCACATCTTCTTCCTGCGCGCCTCCGCCAATGTTATGTATCACCAGCGGCCTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCa                                                                                                <  2:1013980/100‑1 (MQ=255)
   aGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGcc                                                                                             <  2:1563071/100‑1 (MQ=255)
    gCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGccc                                                                                            >  2:1696348/1‑100 (MQ=255)
     cacaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCg                                                                                           >  1:1531679/1‑100 (MQ=255)
      acaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGt                                                                                          >  2:1990773/1‑100 (MQ=255)
      acaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGt                                                                                          <  1:1167656/100‑1 (MQ=255)
       caTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCTCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGt                                                                                          <  1:956728/99‑1 (MQ=255)
       caTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGtt                                                                                         >  1:868089/1‑100 (MQ=255)
        aTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTa                                                                                        <  1:40862/100‑1 (MQ=255)
         tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTAt                                                                                       >  1:1024301/1‑100 (MQ=255)
         tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGAGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTAt                                                                                       >  1:87050/1‑100 (MQ=255)
          cttcttCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATc                                                                                      >  1:478104/1‑100 (MQ=255)
          cttcttCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATc                                                                                      <  2:1893260/100‑1 (MQ=255)
           ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATcc                                                                                     <  1:1626082/100‑1 (MQ=255)
           ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATcc                                                                                     <  2:638248/100‑1 (MQ=255)
           tctttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATcc                                                                                     <  2:246906/97‑1 (MQ=255)
              ttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAga                                                                                  >  1:246907/1‑99 (MQ=255)
                 cTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGc                                                                               >  2:45820/1‑100 (MQ=255)
                 cTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGc                                                                               >  1:1591572/1‑100 (MQ=255)
                 cTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGc                                                                               >  1:573463/1‑100 (MQ=255)
                  tGCGCGCCTGCGCCAATGTTATGTAGCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGcc                                                                              <  1:902870/100‑1 (MQ=255)
                  tGCGCGCCTGCGCCAATGTTATGTACCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGcc                                                                              <  2:523513/100‑1 (MQ=255)
                   gcgcgcCTCCGCCAATGTTATGTATCCCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCa                                                                             <  1:959018/100‑1 (MQ=255)
                    cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGcccg                                                                            <  1:1387880/100‑3 (MQ=255)
                        ccTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGgaga                                                                        >  2:1179770/1‑100 (MQ=255)
                          tGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACa                                                                      <  1:75010/100‑1 (MQ=255)
                          tGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACa                                                                      >  2:1777684/1‑100 (MQ=255)
                          tGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACa                                                                      <  2:817728/100‑1 (MQ=255)
                           gcgcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACaa                                                                     <  1:1089104/100‑1 (MQ=255)
                            cgcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAat                                                                    >  2:22883/1‑100 (MQ=255)
                            cgcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAat                                                                    >  2:61844/1‑100 (MQ=255)
                             gcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAata                                                                   <  1:1460022/100‑1 (MQ=255)
                              ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGCCAatat                                                                  >  2:1802443/1‑100 (MQ=255)
                              ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAatat                                                                  >  2:318333/1‑100 (MQ=255)
                              ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAatat                                                                  >  2:661546/1‑100 (MQ=255)
                               cAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATc                                                                 <  1:22883/100‑1 (MQ=255)
                                  tGTTATGTATCACCATCGTCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                                                              <  2:26150/100‑1 (MQ=255)
                                   gTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                                                             <  2:986035/100‑1 (MQ=255)
                                    ttATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCCGTCGCCAGGAGACAATATCGCCCg                                                            >  1:240549/1‑100 (MQ=255)
                                     tATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                                                           >  2:1870946/1‑100 (MQ=255)
                                     tATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                                                           >  2:817631/1‑100 (MQ=255)
                                      aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCt                                                          >  2:1110116/1‑100 (MQ=255)
                                       tGTATCACCAGCGGCGTCACGTAGAGGGCGACGCCATATGAACCCACCAAAATATGTCGCAGCCCTTTATCTAGCCGCGAGGAGACAATATCTCCCTCtt                                                         <  2:1704159/100‑1 (MQ=255)
                                       tGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCtt                                                         <  2:12237/100‑1 (MQ=255)
                                        gTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATGTCGCCCGCttn                                                        <  2:1420507/100‑2 (MQ=255)
                                        gTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTg                                                        <  2:1478391/100‑1 (MQ=255)
                                        gTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTg                                                        <  2:868086/100‑1 (MQ=255)
                                         tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                                                       >  1:400636/1‑100 (MQ=255)
                                          aTCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                                                      <  1:28530/100‑1 (MQ=255)
                                          aTCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                                                      <  2:972235/100‑1 (MQ=255)
                                          aTCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                                                      <  2:87049/100‑1 (MQ=255)
                                           tCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa                                                     >  2:529468/1‑100 (MQ=255)
                                           tCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa                                                     >  1:482816/1‑100 (MQ=255)
                                            cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                    <  1:7705/100‑1 (MQ=255)
                                             aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGt                                                   >  1:268121/1‑100 (MQ=255)
                                              ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTc                                                  <  2:478103/100‑1 (MQ=255)
                                                aGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCAc                                                <  1:2020888/100‑1 (MQ=255)
                                                    gcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg                                            <  1:1677656/100‑1 (MQ=255)
                                                     cgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg                                           <  2:172878/100‑1 (MQ=255)
                                                       tACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGtt                                         >  2:1933794/1‑100 (MQ=255)
                                                         ccGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTcn                                       <  2:1796657/100‑2 (MQ=255)
                                                            tCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGc                                    <  1:1777698/100‑1 (MQ=255)
                                                             cgcgGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCt                                   >  2:1897753/1‑100 (MQ=255)
                                                               cgGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtgg                                 >  2:1454557/1‑100 (MQ=255)
                                                                   cGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTggg                             >  2:1393425/1‑100 (MQ=255)
                                                                   cGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTggg                             <  2:393770/100‑1 (MQ=255)
                                                                   cGAAGCCATCTGAAACCACCCAAATATTTGGAAGCCCTTTATCTAGTGGCCAGGAGACAAAAGCGCCCGCTTGTTAGTCACTGGGGTTCTTGTTGGTggg                             >  2:1321654/1‑100 (MQ=255)
                                                                    gAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTgggg                            <  1:601700/100‑1 (MQ=255)
                                                                    gAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTgggg                            <  1:1870963/100‑1 (MQ=255)
                                                                     aaGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCACCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGc                           >  1:1801141/1‑100 (MQ=255)
                                                                         cATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgtgt                       <  1:1933812/100‑1 (MQ=255)
                                                                         cATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgagt                       >  1:1132380/1‑100 (MQ=255)
                                                                          aTCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGtt                      <  2:400635/100‑1 (MQ=255)
                                                                          aTCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGtt                      <  1:61844/100‑1 (MQ=255)
                                                                           tCTGAAACCACCCCCCTATGTGGCCGCCCGTTATAAAGTCGCCAGGAGACAATACCGCCCGCTTGATAGTAACTCGGGTCATTGCTGGTGGGGCGTGttn                     <  2:808646/100‑2 (MQ=255)
                                                                            cTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGtttt                    <  1:1955670/100‑1 (MQ=255)
                                                                             tGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                   <  1:1990791/100‑1 (MQ=255)
                                                                             tGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                   >  1:17369/1‑100 (MQ=255)
                                                                              gAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTAt                  <  1:1464016/100‑1 (MQ=255)
                                                                              gAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTAt                  >  1:1065783/1‑100 (MQ=255)
                                                                              gAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTAt                  <  1:260584/100‑1 (MQ=255)
                                                                               aaaCCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATc                 <  1:1696361/100‑1 (MQ=255)
                                                                                aaCCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCa                >  1:1506207/1‑100 (MQ=255)
                                                                                aaCCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCa                >  1:425144/1‑100 (MQ=255)
                                                                                aaCCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCa                <  1:143454/100‑1 (MQ=255)
                                                                                aaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCa                >  2:286310/1‑100 (MQ=255)
                                                                                 accaccCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCAt               <  1:1802459/100‑1 (MQ=255)
                                                                                   caccCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATgg             <  2:1537157/100‑1 (MQ=255)
                                                                                     ccccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg           >  2:1062014/1‑100 (MQ=255)
                                                                                       ccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcggc         <  2:729132/100‑1 (MQ=255)
                                                                                       ccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcggc         <  1:1393434/100‑1 (MQ=255)
                                                                                       ccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcggc         >  1:1012411/1‑100 (MQ=255)
                                                                                          atatGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTaa      <  2:1117825/100‑1 (MQ=255)
                                                                                             tgtgGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc   >  2:1867511/1‑100 (MQ=255)
                                                                                             tgtgGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc   >  1:1836282/1‑100 (MQ=255)
                                                                                             tgtgGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc   <  1:1897771/100‑1 (MQ=255)
                                                                                              gtgGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTTCTGGTGGGGCGTGTTTTATCATGGCGGCTCAAcc   >  2:1953557/1‑99 (MQ=255)
                                                                                              gtgGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa  >  2:1379621/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCA  >  NC_000913/4137083‑4137276

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