Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 91,043 A→G R4R (AGA→AGG ftsL → cell division protein FtsL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091391,0430AG98.5% 240.3 / ‑3.5 65coding (12/366 nt)ftsLcell division protein FtsL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (39/25);  minor base . (1/0);  total (40/25)
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.96e-01

CAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAC  >  NC_000913/90907‑91177
                                                                                                                                        |                                                                                                                                      
cagcTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTg                                                                                                                                     >  2:982166/1‑140 (MQ=255)
   cTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACa                                                                                                                                  >  2:199218/1‑140 (MQ=255)
    tCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAg                                                                                                                                 <  1:1091403/140‑1 (MQ=255)
     cAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGa                                                                                                                                <  2:1072258/140‑1 (MQ=255)
        aaaaCTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGc                                                                                                                             <  1:199218/140‑1 (MQ=255)
             tGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTaa                                                                                                                        >  1:678756/1‑140 (MQ=255)
              gggTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAg                                                                                                                       <  1:114682/140‑1 (MQ=255)
                gTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGACGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGc                                                                                                                      >  1:1155724/1‑139 (MQ=255)
                gTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGc                                                                                                                      <  2:1155724/139‑1 (MQ=255)
                gTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCa                                                                                                                     <  1:113214/140‑1 (MQ=255)
                 tGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCaa                                                                                                                    >  2:738495/1‑140 (MQ=255)
                  ggCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCaaa                                                                                                                   >  1:1107402/1‑140 (MQ=255)
                      gTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTa                                                                                                               <  2:685930/140‑1 (MQ=255)
                       tCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTaa                                                                                                              >  1:960884/1‑140 (MQ=255)
                            tGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGAt                                                                                                         <  1:1208994/140‑1 (MQ=255)
                              cGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCg                                                                                                       >  1:311036/1‑140 (MQ=255)
                                     tAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGaa                                                                                                <  1:281450/140‑1 (MQ=255)
                                      aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATgg                                                                                                   <  1:613954/136‑1 (MQ=255)
                                      aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATgg                                                                                                   >  2:613954/1‑136 (MQ=255)
                                      aGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAg                                                                                               >  2:1240942/1‑140 (MQ=255)
                                          aaGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCAc                                                                                           <  1:738495/140‑1 (MQ=255)
                                                 tGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGcc                                                                                    >  2:537411/1‑140 (MQ=255)
                                                 tGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGcc                                                                                    >  2:966071/1‑140 (MQ=255)
                                                  gCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCa                                                                                   >  2:1138447/1‑140 (MQ=255)
                                                      ggCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               >  1:450703/1‑140 (MQ=255)
                                                      ggCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               <  1:1240942/140‑1 (MQ=255)
                                                      ggCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               <  1:208026/140‑1 (MQ=255)
                                                      ggCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               >  1:394799/1‑140 (MQ=255)
                                                       gCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGATGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               <  1:307353/139‑1 (MQ=255)
                                                       gCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGATGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               >  2:307353/1‑139 (MQ=255)
                                                       gCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCAt                                                                              >  2:165624/1‑140 (MQ=255)
                                                           agaagaGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCAtt                                                                             <  1:781933/137‑1 (MQ=255)
                                                           agaagaGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCAtt                                                                             >  2:781933/1‑137 (MQ=255)
                                                           agaagaGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGcc                                                                          <  1:982166/140‑1 (MQ=255)
                                                            gaagaGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCt                                                                         <  2:450703/140‑1 (MQ=255)
                                                            gaagaGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCt                                                                         <  1:978588/140‑1 (MQ=255)
                                                                     gCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               >  2:481781/1‑125 (MQ=255)
                                                                     gCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCa                                                                               <  1:481781/125‑1 (MQ=255)
                                                                     gCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATc                                                                >  1:573855/1‑140 (MQ=255)
                                                                     gCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATc                                                                >  2:150644/1‑140 (MQ=255)
                                                                         agaACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTg                                                            >  1:498018/1‑140 (MQ=255)
                                                                          gaACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGa                                                           >  1:1256776/1‑140 (MQ=255)
                                                                            aCCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACg                                                         >  1:1035043/1‑140 (MQ=255)
                                                                                   tGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCtttt                                                  >  1:685966/1‑140 (MQ=255)
                                                                                               agttcTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGaa                                      >  1:795901/1‑140 (MQ=255)
                                                                                                  tcTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCt                                   <  1:367678/140‑1 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGa                                       <  1:1125743/132‑1 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCt                                   >  1:440173/1‑136 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCt                                   <  2:440173/136‑1 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCa                               >  2:930221/1‑140 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCa                               >  2:628659/1‑140 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCa                               >  2:268037/1‑140 (MQ=255)
                                                                                                      cGTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGAACGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGa                                       >  2:1125743/1‑132 (MQ=255)
                                                                                                       gTATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGTCAc                              <  2:685966/140‑1 (MQ=255)
                                                                                                        tATTGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCAct                             <  2:960884/140‑1 (MQ=255)
                                                                                                           tGCAGAGAGGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTg                          >  1:536335/1‑140 (MQ=255)
                                                                                                                agagGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCAct                             >  1:890777/1‑132 (MQ=255)
                                                                                                                agagGACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCAct                             <  2:890777/132‑1 (MQ=255)
                                                                                                                    gACGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCAt                 >  2:334116/1‑140 (MQ=255)
                                                                                                                     aCGAATGCATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCAtt                >  2:1202441/1‑140 (MQ=255)
                                                                                                                          tGCATGATCAGCAGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCatt          >  1:399457/1‑140 (MQ=255)
                                                                                                                            cATGATCAGCAGGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCatta         >  1:1216203/1‑140 (MQ=255)
                                                                                                                            cATGATCAGCAGGGGGACAGAAGCTCTAAGCATAGATAAAGGATCGATGGGAAGCCACGAGCGCCATGCATCGCCTGGGGTTATCGGTGACGATCTTTTGCTTCTTGGAAAGCTGCCACTCTGCCTGTGCAATTGCatta         >  2:1181089/1‑140 (MQ=255)
                                                                                                                                  cagcagGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGa   >  1:937560/1‑140 (MQ=255)
                                                                                                                                   agcagGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAc  >  1:178697/1‑140 (MQ=255)
                                                                                                                                   agcagGGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAc  >  1:410317/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CAGCTCAAAAAACTGGGTGGCCGTCAGCTGCGAGCACTAGGCAAGTTAATGCCGGGCGAAGAAGAGGTGGCTGAGAACCCTCGTGCCCGTAGTTCAGTTCTGCGTATTGCAGAGAGGACGAATGCATGATCAGCAGAGTGACAGAAGCTCTAAGCAAAGTTAAAGGATCGATGGGAAGCCACGAGCGCCATGCATTGCCTGGTGTTATCGGTGACGATCTTTTGCGATTTGGGAAGCTGCCACTCTGCCTGTTCATTTGCATTATTTTGAC  >  NC_000913/90907‑91177

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: