Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,613,431 A→G 5.8% S155G (AGC→GGC)  acpT → 4'‑phosphopantetheinyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,613,4310AG 5.8% 222.3 / 2.0 73S155G (AGC→GGC) acpT4'‑phosphopantetheinyl transferase
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base G (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCGCCCACCCGACTGCCCATC  >  NC_000913/3613287‑3613578
                                                                                                                                                |                                                                                                                                                   
gCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCAc                                                                                                                                                 <  2:245959/149‑1 (MQ=255)
gCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCAc                                                                                                                                                 <  4:134972/149‑1 (MQ=255)
  cGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTGTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACATCACCt                                                                                                                                               <  8:131400/149‑1 (MQ=255)
       tCAGCCTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCAc                                                                                                                                          >  5:226100/1‑149 (MQ=255)
            cTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCActcctc                                                                                                                                     >  4:105281/1‑149 (MQ=255)
                 ggAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTgt                                                                                                                                <  6:131681/149‑1 (MQ=255)
                  gAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTAt                                                                                                                                             >  2:320543/1‑135 (MQ=255)
                  gAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGGCAGCGTAGACAGCACCTAt                                                                                                                                             <  1:320543/135‑1 (MQ=255)
                     cacGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTgtcagt                                                                                                                            <  8:200049/149‑1 (MQ=255)
                                 gACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCa                                                                                                                <  6:226100/149‑1 (MQ=255)
                                   cGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGc                                                                                                              >  5:183489/1‑149 (MQ=255)
                                      cGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCCGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCtt                                                                                                            >  8:337131/1‑148 (MQ=255)
                                        tGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGa                                                                                                          >  8:333938/1‑148 (MQ=255)
                                          cATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCTTTAAACTGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGaaa                                                                                                        <  4:312355/148‑1 (MQ=255)
                                           aTCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAAt                                                                                                      >  4:308152/1‑149 (MQ=255)
                                           aTCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGGTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAAt                                                                                                      >  3:56252/1‑149 (MQ=255)
                                             ccTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAAtt                                                                                                     >  8:232045/1‑148 (MQ=255)
                                                gATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAg                                                                                                 >  2:137072/1‑149 (MQ=255)
                                                 aTCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGc                                                                                                >  3:264844/1‑149 (MQ=255)
                                                 aTCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCGTGAAAATTTaa                                                                                                  <  4:51231/147‑1 (MQ=255)
                                                     gCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCCCTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTg                                                                                            >  1:347440/1‑149 (MQ=255)
                                                     gCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTTAAAATTTAAGCCTTg                                                                                            >  5:400654/1‑149 (MQ=255)
                                                          tGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATc                                                                                       >  3:126344/1‑149 (MQ=255)
                                                            gAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACGGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTg                                                                                     >  3:226868/1‑149 (MQ=255)
                                                            gAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACGGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTg                                                                                     >  3:231934/1‑149 (MQ=255)
                                                                tGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCAcc                                                                                 <  4:328402/149‑1 (MQ=255)
                                                                 gTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCAccc                                                                                >  2:186060/1‑149 (MQ=255)
                                                                  ttCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTaa                                                                                                  >  2:381814/1‑130 (MQ=255)
                                                                  ttCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTaa                                                                                                  <  1:381814/130‑1 (MQ=255)
                                                                  ttCTGGCGCATCTGGACGAGAAAATAAGCCATCGTTAAACAGCGTGGCGGAAGAGCCTGGCAAATCGTCAGCGTAGACAGCACCTATAACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCAccc                                                                                <  6:400654/148‑1 (MQ=255)
                                                                    cTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCt                                                                              <  3:308152/148‑1 (MQ=255)
                                                                       gcgcATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCa                                                                                   <  6:361971/140‑1 (MQ=255)
                                                                       gcgcATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCa                                                                                   >  5:361971/1‑140 (MQ=255)
                                                                       gcgcATCCGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACt                                                                           <  7:337131/148‑1 (MQ=255)
                                                                        cgcATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTcc                                                                         <  7:232045/149‑1 (MQ=255)
                                                                        cgcATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTcc                                                                         <  3:105281/149‑1 (MQ=255)
                                                                            tCTGGACGCGCAAAGACGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCtt                                                                      >  1:37147/1‑148 (MQ=255)
                                                                            tCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTAt                                                                                                                                             <  4:371542/77‑1 (MQ=255)
                                                                            tCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCtt                                                                      <  2:37147/148‑1 (MQ=255)
                                                                                gACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTAc                                                                   <  6:183489/147‑1 (MQ=255)
                                                                                gACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACg                                                                  <  7:346313/148‑1 (MQ=255)
                                                                                     gcAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCAcc                                                            >  8:137306/1‑149 (MQ=255)
                                                                                               cATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACACCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAgtg                                                  >  5:232484/1‑149 (MQ=255)
                                                                                                aTCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAgtgt                                                 >  7:201671/1‑149 (MQ=255)
                                                                                                 tCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGtgtg                                                <  5:88583/149‑1 (MQ=255)
                                                                                                 tCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGtgtg                                                >  6:88583/1‑149 (MQ=255)
                                                                                                    ttAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCtt                                                                      <  8:132847/124‑1 (MQ=255)
                                                                                                    ttAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCtt                                                                      >  7:132847/1‑124 (MQ=255)
                                                                                                       aaCAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATg                                           <  1:186060/148‑1 (MQ=255)
                                                                                                         cAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGAt                                        <  4:364635/149‑1 (MQ=255)
                                                                                                          aGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCaa                                             <  4:269988/143‑1 (MQ=255)
                                                                                                          aGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATAGa                                         >  3:269988/1‑147 (MQ=255)
                                                                                                            cGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTACCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGa                                     >  4:144109/1‑149 (MQ=255)
                                                                                                               ggcggcAGCGCCTGGCAAATCGTCAGCGTAGCCAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGAtt                                   >  6:176113/1‑148 (MQ=255)
                                                                                                                gcggcAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGAtt                                   >  8:413041/1‑147 (MQ=255)
                                                                                                                gcggcAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCa                                 <  7:333938/149‑1 (MQ=255)
                                                                                                                 cggcAGCGCCTGGCAAATCGTCAGCGTAGACAGCAACTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAg                                <  8:266946/149‑1 (MQ=255)
                                                                                                                   gcAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGtt                              <  8:201671/149‑1 (MQ=255)
                                                                                                                      gcgcCTGGCAAATCGTCAGCGTAGACGGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAAc                           <  4:226868/149‑1 (MQ=255)
                                                                                                                      gcgcCTGGCAAATCGTCAGCGTAGACGGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAAc                           <  4:231934/149‑1 (MQ=255)
                                                                                                                      gcgcCTGGCAAATCGTCAGCGTAGACAGCACCTATAACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCCCCCCTAATCCCTTTACGCTCACCGCCTACAGTGTGCAATGGAt                                        >  1:127786/1‑136 (MQ=255)
                                                                                                                        gggtGGCAATTCGTCAGCGTAGAGTGCACCTATCACTCCTCGCTGTCAGTCAGCGATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACGCCTACTCGCTTTACGCTGACCGACGTCAGTGTGCTTTGGATGGATTCAGTTTACt                          >  7:234478/4‑148 (MQ=255)
                                                                                                                           tGGCAAATCGTCAGAGTAGACCGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATc                      >  4:79726/1‑149 (MQ=255)
                                                                                                                             gCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATc                      <  6:232484/147‑1 (MQ=255)
                                                                                                                             gCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCg                    >  2:237526/1‑149 (MQ=255)
                                                                                                                                    gTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTTTAAGCTTGAAAATTTAAGCATTGAGATCTTCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCGCCCAccc             <  4:264844/149‑1 (MQ=255)
                                                                                                                                    gTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTg                                                                                            >  6:365691/1‑70 (MQ=255)
                                                                                                                                    gTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTg                                                                                            <  5:365691/70‑1 (MQ=255)
                                                                                                                                     tCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGGTAACTGATCCGCCCACCCg            >  6:319500/1‑149 (MQ=255)
                                                                                                                                         cGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCGCCCACCCGACTg        >  1:333954/1‑149 (MQ=255)
                                                                                                                                          gTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCGCCCACCCGACTGc       >  6:338827/1‑149 (MQ=255)
                                                                                                                                              aCAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACg                                                                  <  5:113473/86‑1 (MQ=255)
                                                                                                                                              aCAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACg                                                                  >  6:113473/1‑86 (MQ=255)
                                                                                                                                                aGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCGCCCACCCGACTGCCCATc  >  4:325879/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCGTATTCAGCCTCGGGGAACACGCTGAGATGGACGCCGTGCATCCTGATCAGCAACTGGAAATGTTCTGGCGCATCTGGACGCGCAAAGAAGCCATCGTTAAACAGCGTGGCGGCAGCGCCTGGCAAATCGTCAGCGTAGACAGCACCTATCACTCCTCGCTGTCAGTCAGCCATTGTCAGCTTGAAAATTTAAGCCTTGCGATCTGCACCCCTACTCCCTTTACGCTCACCGCCGACAGTGTGCAATGGATCGATTCAGTTAACTGATCCGCCCACCCGACTGCCCATC  >  NC_000913/3613287‑3613578

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