Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,351,196 (G)7→8 intergenic (‑157/+211) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,1891.G96.3% 192.1 / 1.5 54intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (28/24);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTG  >  NC_000913/1351040‑1351333
                                                                                                                                                      |                                                                                                                                                   
aGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCa                                                                                                                                                      <  1:927627/150‑1 (MQ=255)
 gCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑gg                                                                                                                                                   >  2:1153139/1‑151 (MQ=255)
       aTCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑gggggggg                                                                                                                                             >  2:150056/1‑150 (MQ=255)
          cTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑ggggggggcgc                                                                                                                                          >  2:762021/1‑147 (MQ=255)
          cTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑ggggggggcgc                                                                                                                                          >  1:102347/1‑147 (MQ=255)
            tATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑ggggggggcgc                                                                                                                                          <  2:802702/149‑5 (MQ=255)
             aTTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCttt                                                                                                                                       >  2:366915/1‑151 (MQ=255)
                 ttgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTccc                                                                                                                                   >  1:218793/1‑151 (MQ=255)
                       cTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCAcc                                                                                                                             >  1:787404/1‑151 (MQ=255)
                         tgttgtGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCAccc                                                                                                                            >  1:1130587/1‑150 (MQ=255)
                         tgttgtGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑t                                                                                                                          >  1:227974/1‑151 (MQ=255)
                         tgttgtGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGCCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑t                                                                                                                          >  2:992048/1‑151 (MQ=255)
                           ttgtGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑Ta                                                                                                                         >  1:519729/1‑150 (MQ=255)
                             gtgCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAAc                                                                                                                       <  1:876121/150‑1 (MQ=255)
                                 cTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCct‑c                                                                                                                  >  2:110025/1‑150 (MQ=255)
                                             cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGc                                                                                                     <  1:570602/151‑1 (MQ=255)
                                                  cgtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAggg                                                                                                >  2:455520/1‑151 (MQ=255)
                                                    tcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTTCCCTCACCC‑TAACCCC‑CTCCCCAGAGGGGCgggggg                                                                                               >  1:1138055/1‑150 (MQ=255)
                                                    tcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCCCAACC‑TAACCCC‑CTCCCCCTAGGGGCGAgggg                                                                                               >  2:1068057/1‑150 (MQ=255)
                                                     cgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGa                                                                                              <  1:1171225/150‑1 (MQ=255)
                                                      gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCTCCCTTAACCCC‑TCCCCCAGAgggcgggggg                                                                                                >  1:791223/1‑141 (MQ=255)
                                                           aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                       >  2:405400/1‑151 (MQ=255)
                                                                 cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCaa                                                                                 <  1:661064/151‑1 (MQ=255)
                                                                     aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAatat                                                                             <  1:1096016/151‑1 (MQ=255)
                                                                      cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                            <  1:927333/151‑1 (MQ=255)
                                                                        aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCCCCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGt                                                                          >  2:245590/1‑151 (MQ=255)
                                                                         aCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt                                                                         <  2:755976/151‑1 (MQ=255)
                                                                          cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTa                                                                        <  2:750811/151‑1 (MQ=255)
                                                                              tAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCCGAGGGGCGGGGGGACCGGATTGTGCAAATATTGTTAcccc                                                                    >  1:1171674/1‑151 (MQ=255)
                                                                               aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcccc                                                                    >  1:921919/1‑150 (MQ=255)
                                                                               aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcccc                                                                    >  1:880759/1‑150 (MQ=255)
                                                                                aaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAg                                                                  <  2:328586/151‑1 (MQ=255)
                                                                                 aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAAACCC‑CTCCCCAGAGGGGCGAGGGGACCGGATTGTGCAAAAATTGTTAcccc                                                                    >  2:17090/1‑148 (MQ=38)
                                                                                   ggTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCCCCC‑TAACCCT‑CTCCCCCGAGGGGCGGGGGGACCGGATTGTGCCAAAATTTTTAcccc                                                                    >  2:902708/1‑146 (MQ=255)
                                                                                    gTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaac                                                              <  2:662285/151‑1 (MQ=255)
                                                                                       attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaacaa                                                            <  2:675694/150‑1 (MQ=255)
                                                                                       attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                                           <  2:570370/151‑1 (MQ=255)
                                                                                       attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCCGTTTTTCCCCCCACCC‑TAACCCC‑CCCCCCAAAGGGGGGGGGGGACCGTATTTTGCAAAAATTGTTACCCCACcaacaa                                                            >  2:522547/1‑150 (MQ=255)
                                                                                            aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                                       <  2:102373/150‑1 (MQ=255)
                                                                                              tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCCGAGGGGCGAGGGGACCGCATTGTGCAAATATTGTTACCCCAGCaacaaaaa                                                         >  1:10402/1‑144 (MQ=255)
                                                                                              tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAAACCT‑CTCCCCAGAGGGGGGAGGGGACCGTATTGTGCAAAAATTGTTAACCCCACAAAAAACAGGGTc                                                    >  1:752306/1‑151 (MQ=255)
                                                                                                     ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACa                                               <  2:223975/149‑1 (MQ=255)
                                                                                                     ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACa                                               <  2:177581/149‑1 (MQ=255)
                                                                                                       cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCCGAGGGGCGAGGGGGCCGGATTGTGCAAATATTGTTACCCCCGCAACAAACAGGCGCAAACAGccc                                           >  1:783232/1‑151 (MQ=38)
                                                                                                           gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTa                                        >  1:65640/1‑150 (MQ=255)
                                                                                                                gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TTACCCT‑CTCCCCAGAGGGGCGAGGGGCCCGTATTGTGCAAATATTGTTACCCCCGCAACAAACAGGCTCATACAACCCCTAAcc                                     >  2:178098/1‑148 (MQ=38)
                                                                                                                  tataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TCACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCttt                                 >  1:142286/1‑150 (MQ=255)
                                                                                                                     aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                             <  2:786730/151‑1 (MQ=255)
                                                                                                                     aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCCCATACAGCCCCTAACCCTTTCAt                              <  2:218838/150‑1 (MQ=255)
                                                                                                                      gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATgg                            <  2:160754/151‑1 (MQ=255)
                                                                                                                       ttGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCCGAGGGGCGAGGGGGCCGTATTTTGCAAATATTGTTACCCCAGCAAAAAACAGGCTCATACAGCCCCTAAACCCTTCATgg                            >  1:1051726/1‑150 (MQ=38)
                                                                                                                           cAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATg                       <  2:219551/151‑1 (MQ=255)
                                                                                                                              gTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGG‑CGCTTTTTCCCCCCCCC‑CAACCCA‑CTCCCCCAAGAGGGGGAGGGACCCTCTTGTGCGAAAATTGTTATCCCCACAACAAAAAAGGTCATTCCGACCCCAACCCCTTCAttggcctgg                      >  2:1036489/1‑140 (MQ=255)
                                                                                                                                  gACCGGAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGGATTGTGCAAAAATTGTTACCCCAGCAACAAACAGGCCCATACAGCCCCTAACCCTTTCATGGCGATGGCTgg                  >  1:902375/1‑149 (MQ=255)
                                                                                                                                      ggAGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGt            <  2:221809/151‑1 (MQ=255)
                                                                                                                                        aGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTc          <  2:692792/151‑1 (MQ=255)
                                                                                                                                        aGCAGGCTGCGGCAGGGGGGGG‑CGCTTTTCCCCCCTTCC‑TAAACCTCCTCCCCCG‑GGGGCGGGGGGAACGGATTGTGCAAATATTGTTTACCCCGCAACAAAAAGGCTCAAAAAACCCCCAAACCTTTTATGGCGATGGCTGGGACgg             >  2:722953/1‑148 (MQ=255)
                                                                                                                                                gCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt   <  2:229636/150‑1 (MQ=255)
                                                                                                                                                gCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCtt   <  1:761848/150‑1 (MQ=255)
                                                                                                                                                gCGGCAGGGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTg  <  1:1068941/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AGCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGG‑CGCTTTTCCCCTCACCC‑TAACCCT‑CTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTG  >  NC_000913/1351040‑1351333

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