Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 187,251 G→A T467M (ACG→ATG)  glnD ← uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913187,2510GA98.7% 293.1 / ‑4.1 76T467M (ACG→ATG) glnDuridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (40/35);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

ATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTT  >  NC_000913/187104‑187395
                                                                                                                                                   |                                                                                                                                                
atatCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCatgt                                                                                                                                               >  2:557883/1‑151 (MQ=255)
       gAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTg                                                                                                                                        >  2:678380/1‑151 (MQ=255)
          aCAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCTAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGaa                                                                                                                                     <  2:1004583/151‑1 (MQ=255)
               gcggcgATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGAt                                                                                                                                <  2:125991/151‑1 (MQ=255)
                cggcgATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATc                                                                                                                               <  2:390947/151‑1 (MQ=255)
                 ggcgATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGAt                                                                                                                                <  2:938566/149‑1 (MQ=255)
                  gcgATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCa                                                                                                                              >  2:1012360/1‑150 (MQ=255)
                      tGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTg                                                                                                                         >  1:595865/1‑151 (MQ=255)
                          aaTCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATc                                                                                                                     <  2:414351/151‑1 (MQ=255)
                           aTCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGAGCAAACTGCATCt                                                                                                                    <  2:988740/151‑1 (MQ=255)
                            tCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTg                                                                                                                   >  1:955126/1‑151 (MQ=255)
                             cAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTg                                                                                                                   <  1:363653/150‑1 (MQ=255)
                               gCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGccc                                                                                                                <  2:472679/151‑1 (MQ=255)
                                  cAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgac                                                                                                             <  1:389272/151‑1 (MQ=255)
                                  cAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCcg                                                                                                               <  2:534280/149‑1 (MQ=255)
                                   aGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCcgacg                                                                                                            >  1:888414/1‑151 (MQ=255)
                                      tGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGat                                                                                                          <  2:405154/150‑1 (MQ=255)
                                      tGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGat                                                                                                          >  2:507106/1‑150 (MQ=255)
                                       gACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGatat                                                                                                        >  2:1006725/1‑151 (MQ=255)
                                        aCGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATg                                                                                                       <  1:41907/151‑1 (MQ=255)
                                           gCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCg                                                                                                     >  1:192575/1‑150 (MQ=255)
                                             aGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGatagatcgga                                                                                                  <  1:961070/151‑8 (MQ=255)
                                             aGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGAcc                                                                                                  >  1:619420/1‑151 (MQ=255)
                                                cgcgGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCAtt                                                                                               <  1:384071/151‑1 (MQ=255)
                                                 gcgGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCATTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTg                                                                                              <  1:15638/151‑1 (MQ=255)
                                                  cgGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGc                                                                                             <  2:465473/151‑1 (MQ=255)
                                                    gCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCgg                                                                                           <  2:60742/151‑1 (MQ=255)
                                                    gCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCgg                                                                                           >  2:150941/1‑151 (MQ=255)
                                                     ccAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGc                                                                                          <  2:546364/151‑1 (MQ=255)
                                                     ccAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGc                                                                                          >  1:642146/1‑151 (MQ=255)
                                                       aGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCat                                                                                        >  2:1091746/1‑151 (MQ=255)
                                                       aGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCat                                                                                        >  2:309727/1‑151 (MQ=255)
                                                       aGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCat                                                                                        >  1:841758/1‑151 (MQ=255)
                                                         aCGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCatat                                                                                      >  1:164873/1‑151 (MQ=255)
                                                          cGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCAtata                                                                                     <  2:856687/151‑1 (MQ=255)
                                                           gTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCAtata                                                                                     >  2:5977/1‑150 (MQ=255)
                                                           gTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATAc                                                                                    >  2:692994/1‑151 (MQ=255)
                                                              cacacacaACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATAcgcg                                                                                 >  2:789878/1‑151 (MQ=255)
                                                                    caACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGc                                                                           <  2:986299/151‑1 (MQ=255)
                                                                           tggcgctggcgcGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTAt                                                                    >  2:880149/1‑151 (MQ=255)
                                                                            ggcgctggcgcGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATg                                                                   <  2:978702/151‑1 (MQ=255)
                                                                              cgctggcgcGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGc                                                                 <  1:702300/151‑1 (MQ=255)
                                                                                 tggcgcGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGc                                                                 <  2:1178285/148‑1 (MQ=255)
                                                                                  ggcgcGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATg                                                             <  1:514433/151‑1 (MQ=255)
                                                                                     gcgTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCAt                                                          <  1:1081839/151‑1 (MQ=255)
                                                                                      cgTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCAtt                                                         >  2:1177254/1‑151 (MQ=255)
                                                                                      cgTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCAtt                                                         >  1:468726/1‑151 (MQ=255)
                                                                                       gTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTg                                                        <  1:392695/151‑1 (MQ=255)
                                                                                            ttcACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAAt                                                   >  2:826534/1‑151 (MQ=255)
                                                                                              cACTGGCAGAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAAt                                                   >  1:455514/1‑149 (MQ=255)
                                                                                                 tGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTAGGGCGATGCATTGGCAATAGccc                                              >  1:466597/1‑151 (MQ=255)
                                                                                                      aaaCTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGcgcc                                         >  1:768367/1‑151 (MQ=255)
                                                                                                           ctcCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCAcc                                    >  1:223269/1‑151 (MQ=255)
                                                                                                             ccAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCg                                   <  1:483340/150‑1 (MQ=255)
                                                                                                                 tttCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCg                              >  2:964362/1‑151 (MQ=255)
                                                                                                                   tCAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCgg                             <  1:631722/150‑1 (MQ=255)
                                                                                                                    cAGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCggg                            <  1:870377/150‑1 (MQ=255)
                                                                                                                     aGCATCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGt                          >  2:241515/1‑151 (MQ=255)
                                                                                                                        aTCACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGAc                       >  1:224899/1‑151 (MQ=255)
                                                                                                                          cACGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGc                     >  2:945537/1‑151 (MQ=255)
                                                                                                                           aCGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGc                     >  2:1010812/1‑150 (MQ=255)
                                                                                                                           aCGCGGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCa                    <  2:919530/151‑1 (MQ=255)
                                                                                                                            cgcgGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGc                     >  2:189025/1‑149 (MQ=255)
                                                                                                                              cgGATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGaa                 <  1:150905/151‑1 (MQ=255)
                                                                                                                               ggATAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAAt                >  2:985787/1‑151 (MQ=255)
                                                                                                                                 aTAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGc              >  2:632881/1‑151 (MQ=255)
                                                                                                                                  tAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAg                   <  2:181217/145‑1 (MQ=255)
                                                                                                                                  tAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAg                   >  1:181178/1‑145 (MQ=255)
                                                                                                                                  tAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGc              <  1:903468/150‑1 (MQ=255)
                                                                                                                                  tAGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCt             >  2:241/1‑151 (MQ=255)
                                                                                                                                   aGTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTc            >  2:509029/1‑151 (MQ=255)
                                                                                                                                    gTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCa           <  2:1081222/151‑1 (MQ=255)
                                                                                                                                    gTATGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCa           <  2:32296/151‑1 (MQ=255)
                                                                                                                                      aTGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCaaa         <  2:455625/151‑1 (MQ=255)
                                                                                                                                      aTGTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCaa          >  1:730050/1‑150 (MQ=255)
                                                                                                                                        gTTCATCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCaaaa        >  1:396977/1‑150 (MQ=255)
                                                                                                                                             tCCACCATGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGtt  >  2:506877/1‑151 (MQ=255)
                                                                                                                                              ccaccaTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGtt  <  1:689981/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATATCGTGAAACAGCGCGGCGATGAAAATCAGCTCAGTTGACGGCAGGCGCGGCCAGACGTCCACACACAACGGATGGCGCTGGCGCGTTTCTTCACTGGCAAAACTCTCCAGTTTCAGCATCACGCGGATAGTATGTTCATCCACCGTGTAGGCGTGGAACAGATCAAACTGCATCTGCCCGACGATATGCGACCATTGCGGCATATACGCGCCGAGCACGCTATGGCGATGCATTGGCAATAGCCCGCGCCGCACCGCTCCGGGGTGACGCAGAATGCTCAAAAACAGTT  >  NC_000913/187104‑187395

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