Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,305 C→T R395C (CGT→TGT)  icd → isocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,3050CT97.2% 188.0 / ‑5.9 74R395C (CGT→TGT) icdisocitrate dehydrogenase, e14 prophage attachment site, tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base T (41/28);  minor base G (1/1);  total (45/29)
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

GACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGG  >  NC_000913/1196167‑1196442
                                                                                                                                          |                                                                                                                                         
gACAAAGTAAATCCTGGCTCTATTATTCTCTACGCTGAGATGATGCTGCGCCAAATGGGTTGGACCGAAGCCGCTGACTTAATTGTTAAAGGTATGGCAGGCGCAATCCGCGCGAAAACCCTAACCTAGGACTTCGAGTGTCt                                                                                                                                       >  2:911606/1‑143 (MQ=255)
  cAAAGTAAATCCTGGCTCTATCACCAACACAACACATATTATGCTGCGCCACATGGGTGGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTga                                                                                                                                     <  2:310780/143‑1 (MQ=18)
   aaaGGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgat                                                                                                                                    <  1:736390/143‑1 (MQ=255)
    aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGACAACCGTAACCTCTGACTTCGCGTGTCTgatg                                                                                                                                   >  2:923168/1‑143 (MQ=255)
     aGTAGATCCTGGATCTTTTAGGCTCTCCGGTCATATGATGCGGCGCCTCAGGGGGTGGACCGAAGCGGCTGACTTAATTGTTAACGGTATGGAAGTCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgg                                                                                                                                  <  1:662784/143‑1 (MQ=255)
     aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTgatgg                                                                                                                                  >  2:372677/1‑143 (MQ=255)
      gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCCATCAACGCGACAACCGTAACCTATGACTTCGAGAGtataataga                                                                                                                                 >  2:447591/1‑135 (MQ=25)
         agtCCTGGCTCTAGTATTCTCGCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGatggatgg                                                                                                                              <  1:204088/141‑1 (MQ=255)
         aaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGatggatgg                                                                                                                              <  2:331843/143‑1 (MQ=255)
           tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGCAAACCGTAACCTATGACTTCGAGTGTCTGAGGGATGgcg                                                                                                                            >  1:928052/1‑143 (MQ=255)
                gCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTaaa                                                                                                                       >  2:760063/1‑143 (MQ=255)
                  tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGCTGGCGCTCAAct                                                                                                                     >  1:93647/1‑143 (MQ=255)
                    tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTCATTGTTAAAGGTATGGAAGGCGCAATCAACGCGACAACCGTAACCTATGACTTCGCGTGTCTGATGGATGGCGCTAACctgc                                                                                                                   >  1:768392/1‑143 (MQ=21)
                    tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAActgc                                                                                                                   >  2:371784/1‑143 (MQ=255)
                    tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGGGTCTGATGGATGGCGCTAAActgc                                                                                                                   >  1:175741/1‑143 (MQ=255)
                            tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGGAATgg                                                                                                           >  2:236951/1‑142 (MQ=255)
                              tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGACTGTCTGATGGCTGGCGCTAAACTGCTGAAATGTTc                                                                                                         >  1:176526/1‑143 (MQ=255)
                                 gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCaga                                                                                                      >  1:769369/1‑143 (MQ=255)
                                      gctgatgCTGCGCCACATGAGTTGGACCTCAACAGCTGACTTAACTGTTTGACGTATGAAGGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGTGTGTCTGATTGATCGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                                                                                 <  1:100301/141‑1 (MQ=18)
                                      gatgatGCTGCGCCCCATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                                                                                 <  2:628305/143‑1 (MQ=255)
                                      gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGACATGTTCAGAGTTTg                                                                                                 >  2:90223/1‑143 (MQ=255)
                                      gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                                                                                 <  2:1053902/143‑1 (MQ=255)
                                      gatgagGATGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                                                                                 <  2:1012051/143‑1 (MQ=255)
                                         gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTg                                                                                              >  2:874625/1‑143 (MQ=255)
                                            gctgcGCCACCTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                                                                           <  1:858859/143‑1 (MQ=255)
                                            gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTACACTGCTGCAATGTTCAGAGCTTGGTGAcg                                                                                           >  2:58859/1‑143 (MQ=255)
                                            gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                                                                           >  2:425640/1‑143 (MQ=255)
                                            gctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTAAGAGCGTCTGCCGGATGGCGCTAAACTGCTGCAATGTTCACAGATGGGGGAcg                                                                                           >  1:1130735/1‑143 (MQ=255)
                                             ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAAAGTTCAGAGTTAGGTGAcgc                                                                                          >  2:619323/1‑143 (MQ=255)
                                              tgcTCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGAGGGCGCTAAACTGCTCAAATGTTCAGAGTTTGGTGAcgcg                                                                                         <  2:944784/143‑1 (MQ=255)
                                               gcgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAAGATTCAGAGGTTGGTGGCGCGa                                                                                        >  1:1046560/1‑143 (MQ=18)
                                                cgcCACATCGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGGATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                                                                       <  2:422348/143‑1 (MQ=255)
                                                    acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                                                                   >  2:753732/1‑143 (MQ=255)
                                                     caTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCg                                                                                  >  1:910369/1‑143 (MQ=255)
                                                        gggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaa                                                                               <  2:1009724/143‑1 (MQ=255)
                                                         ggTCGGACCGAAGCGGCTGACCTAATTGTTAAAGGTATGGGAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaaa                                                                              <  1:898652/143‑1 (MQ=255)
                                                               aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGGGTTTGGTGACGCGATCATCGAAACCAGGTa                                                                        >  1:133931/1‑143 (MQ=255)
                                                                ccGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCAAGTGTCTGATGGATGGCGCTAAACTGCTGAAATCTTCAGAGTTTGGGGACGCGagcactgaacaagtgtaa                                                                       >  1:222791/1‑129 (MQ=18)
                                                                  caaGCGCCTAACTTAACTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg                                                                     <  2:399995/142‑1 (MQ=255)
                                                                      cGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGt                                                                 >  2:346966/1‑143 (MQ=255)
                                                                        gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGGTTGGTGACGCGATCACCGAAAACATGGAATGCCGTAg                                                               >  1:737725/1‑143 (MQ=255)
                                                                          tGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGtt                                                             >  2:250056/1‑143 (MQ=255)
                                                                           gatttAATTGTTAATGTGATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGttt                                                            <  2:827113/140‑1 (MQ=18)
                                                                              ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                                         >  2:213920/1‑143 (MQ=255)
                                                                              ttAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                                         >  2:392158/1‑143 (MQ=255)
                                                                                 aTTGTAATAGGTATGGAAGGCGCAATCAACGCGAAACCCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa                                                      <  1:559121/143‑1 (MQ=255)
                                                                                  ttGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAt                                                     >  1:624834/1‑143 (MQ=255)
                                                                                    gTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAttc                                                   >  2:710897/1‑142 (MQ=255)
                                                                                       aagggTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTGCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAACACATGTAATGCCGTAGTTTGTTCAATttatt                                                >  1:1108444/4‑143 (MQ=21)
                                                                                        aaGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATttatta                                               >  1:583995/1‑143 (MQ=255)
                                                                                          ggTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGGGTGTCTGATGGATGGCGCTAAAATCCTGAAATGTTCGGGGTTTGGTGACGCGATCATCGAAAAGAAGTAGAGCGGTAGTTTGttcagcgtattgca                                             >  1:360446/1‑131 (MQ=17)
                                                                                           gTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACg                                            >  2:824142/1‑143 (MQ=255)
                                                                                           gTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACg                                            >  1:916982/1‑143 (MQ=255)
                                                                                                gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGc                                       >  2:1097558/1‑143 (MQ=255)
                                                                                                 aaGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAGGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCg                                      >  1:298306/1‑143 (MQ=255)
                                                                                                  tgtcgcAATGAGCGCGAAAACCGTAACCTCTGACGTCGGCTGTCTTATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAGCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGt                                     <  1:1132280/140‑1 (MQ=21)
                                                                                                    gcgcAATCAACGGGAAACCCGTAACCTATGTCTTCGAGGGTCTGATGGATGCCGGTAAACTGCTGAAATGGTCGGAGTTTGGTGACGCGATAAACGAAAACATGTAATGCCGTAGTGTGTTAAATTTATTAACGGGAGCGTaa                                   <  1:727758/143‑1 (MQ=255)
                                                                                                    gcgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGCAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                   >  1:35185/1‑143 (MQ=255)
                                                                                                        aaTCAACGCGAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCt                               >  1:950580/1‑143 (MQ=255)
                                                                                                         aTCAACGCGAAGACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTGAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGccc                              <  1:965014/143‑3 (MQ=255)
                                                                                                                 gAAAACCGTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAGTTATTAACGGGAGCGTAACGCTGCCgttgtt                      >  2:391329/1‑143 (MQ=255)
                                                                                                                   aaaCCGTAACCTATGACTTCGAGTGCCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTTACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTAGtttt                    <  1:959991/143‑1 (MQ=255)
                                                                                                                    aaCCGTAACCTAGGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGGTAATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGttttt                   <  2:729121/143‑1 (MQ=255)
                                                                                                                     tcggtaacgtcggtcTTCGGGTGTGTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTGATGCCGTCGTTTGTTAAATTTATTAGCGGGAGCGTAACGCTCCCGTTGtttttt                  <  1:29321/129‑1 (MQ=25)
                                                                                                                        gTAACCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGtt               <  2:916982/143‑1 (MQ=255)
                                                                                                                         tatccTCTGACTTCGAGTGTCTGGTGGACGGCGCTAAGCTGCTGAAATGTTCAGCGTTTGGTGACGCGATCATCGAAAACATGTATTGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTa              <  1:330393/140‑1 (MQ=25)
                                                                                                                           aCCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGACAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAGCGCTCCCGTTGTTTTTTGTTAgg            >  1:979315/1‑143 (MQ=255)
                                                                                                                           aCCTATGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg            <  2:562651/143‑1 (MQ=255)
                                                                                                                           aCCTATGACTTCCCGTGTCTGATGGATGGCGCTAAACTGCTTAAATGTTCAGAGTTTGGTGACGCGATCACCGAACACATGTAATGCCGTAGTTTGTTAGATTGATTAACGGGAGCGTAACGATACCTATGttttttactaag            >  1:332962/1‑137 (MQ=21)
                                                                                                                               aTGACTTCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgc        <  2:975688/143‑1 (MQ=255)
                                                                                                                                    ttCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGGGTTTGGTCACACAACCCCCCACTGACGGTGATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACg   <  1:513550/143‑1 (MQ=255)
                                                                                                                                    ttCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAATGGCAGCGGAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACg   >  1:509263/1‑143 (MQ=255)
                                                                                                                                    ttCGAGTGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGTGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACg   >  2:97053/1‑143 (MQ=255)
                                                                                                                                     ccGAGTGTCCGATGGATGCCGCTAAACTCCTAAAATGTTCCCAGCTATGAGCCATTATTATTGAAAACACGTAATGCCTTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg  <  2:842477/142‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
GACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGG  >  NC_000913/1196167‑1196442

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