Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,194,577 2 bp→GC coding (97‑98/774 nt) yjaD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,194,5770CG100.0% 729.0 / NA 220R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (73/147);  total (73/147)
*NC_0009134,194,5780GC99.6% 729.7 / NA 220R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (72/147);  total (73/147)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

ATAATTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGC  >  NC_000913/4194490‑4194664
                                                                                       ||                                                                                      
 taATTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAAATGCCATATGGCGAAGGGGCaa                                                                                   >  1:1409032/1‑93 (MQ=255)
  aaTTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                                                  <  2:2433268/93‑1 (MQ=255)
  aaTTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                                                  <  1:3178563/93‑1 (MQ=255)
  aaTTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCaaa                                                                                  <  2:6729127/93‑1 (MQ=255)
   aTTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAt                                                                                 >  1:6989847/1‑93 (MQ=255)
    ttGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                                                >  2:5771378/1‑93 (MQ=255)
      gAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTc                                                                              <  1:5510433/93‑1 (MQ=255)
       aaaaaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCg                                                                             <  1:6169434/93‑1 (MQ=255)
         aaaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGAt                                                                           <  1:6347866/93‑1 (MQ=255)
         aaaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGAt                                                                           <  1:6521392/93‑1 (MQ=255)
         aaaTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGAt                                                                           <  2:1737379/93‑1 (MQ=255)
          aaTTAGACGACGGCTGGTGGGTCGTCAGCCCTCAACAAAAATTATGGTTGCCGAGGGGAGAAATGCCCTATGGGGGAGCGGCAAATTTccaag                                                                          >  1:4259485/1‑89 (MQ=255)
           aTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCt                                                                         >  1:2302328/1‑93 (MQ=255)
           aTTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCt                                                                         <  2:3966720/93‑1 (MQ=255)
             tAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtg                                                                       <  1:6842044/93‑1 (MQ=255)
              aGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgt                                                                      <  2:2199943/93‑1 (MQ=255)
              aGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgt                                                                      <  1:3402886/93‑1 (MQ=255)
               gATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgtg                                                                     <  1:5238370/93‑1 (MQ=255)
               gATCACGGCTGGTGGGTCATCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgtg                                                                     <  1:5545216/93‑1 (MQ=255)
                aTCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                                    <  2:2074257/93‑1 (MQ=255)
                 tCACGGCTTGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   >  2:4861307/1‑93 (MQ=255)
                 tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   >  2:5272292/1‑93 (MQ=255)
                 tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   <  2:5506513/93‑1 (MQ=255)
                 tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   >  2:1378328/1‑93 (MQ=255)
                 tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                                   <  1:1084605/93‑1 (MQ=255)
                   aCGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                                                 <  1:2901695/93‑1 (MQ=255)
                      gCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                                              <  2:1900206/93‑1 (MQ=255)
                      gCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                                              <  2:3512909/93‑1 (MQ=255)
                      gCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                                              >  1:2465318/1‑93 (MQ=255)
                           tggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCa                                                         <  1:2238581/93‑1 (MQ=255)
                           tggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCa                                                         >  1:5184358/1‑93 (MQ=255)
                           tggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCa                                                         <  2:4499076/93‑1 (MQ=255)
                            gggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCAc                                                        <  2:6346308/93‑1 (MQ=255)
                             ggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACt                                                       <  2:3584919/93‑1 (MQ=255)
                              gtcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTa                                                      <  1:6673274/93‑1 (MQ=255)
                               tcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTAc                                                     >  1:2073463/1‑93 (MQ=255)
                                cgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACa                                                    >  1:3333694/1‑93 (MQ=255)
                                 gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaa                                                   >  1:4690106/1‑93 (MQ=255)
                                  tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaaa                                                  <  1:7045912/93‑1 (MQ=255)
                                  tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaaa                                                  >  1:2761902/1‑93 (MQ=255)
                                  tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaaa                                                  >  2:3034303/1‑93 (MQ=255)
                                  tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaaa                                                  <  1:455843/93‑1 (MQ=255)
                                   cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                                                 <  1:59583/93‑1 (MQ=255)
                                   cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                                                 <  1:5795188/93‑1 (MQ=255)
                                    aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATc                                                <  2:7033643/93‑1 (MQ=255)
                                    aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATc                                                >  2:1134178/1‑93 (MQ=255)
                                    aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATc                                                <  2:491859/93‑1 (MQ=255)
                                     gCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCg                                               <  1:254951/93‑1 (MQ=255)
                                     gCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCg                                               >  2:2551703/1‑93 (MQ=255)
                                      ccATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCgg                                              >  2:5682005/1‑93 (MQ=255)
                                      ccATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCgg                                              >  2:7505744/1‑93 (MQ=255)
                                       cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGc                                             <  1:1977292/93‑1 (MQ=255)
                                        aTGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCg                                            >  1:3430200/1‑93 (MQ=255)
                                        aTGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCg                                            >  1:2175565/1‑93 (MQ=255)
                                         tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGa                                           >  2:5614304/1‑93 (MQ=255)
                                           aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                                         >  2:1035924/1‑93 (MQ=255)
                                            aCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATg                                        <  2:4534035/93‑1 (MQ=255)
                                            aCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATg                                        <  2:6073493/93‑1 (MQ=255)
                                              aaaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGc                                      <  1:297445/93‑1 (MQ=255)
                                               aaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCa                                     <  1:884806/93‑1 (MQ=255)
                                                 aaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                                   <  2:7117992/93‑1 (MQ=255)
                                                 aaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                                   >  2:5158742/1‑93 (MQ=255)
                                                  aTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggg                                  >  1:2595615/1‑93 (MQ=255)
                                                  aTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggg                                  <  1:526972/93‑1 (MQ=255)
                                                   ttATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg                                 <  1:2634879/93‑1 (MQ=255)
                                                   ttATGGATGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg                                 <  2:2738702/93‑1 (MQ=255)
                                                    tATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggggg                                <  1:4002029/93‑1 (MQ=255)
                                                      tGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGaa                              <  1:7466187/93‑1 (MQ=255)
                                                      tGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGaa                              <  2:1119419/93‑1 (MQ=255)
                                                      tGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGaa                              >  2:3643822/1‑93 (MQ=255)
                                                       ggTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAc                             <  2:4441854/93‑1 (MQ=255)
                                                        gTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAcc                            <  1:3556716/93‑1 (MQ=255)
                                                        gTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAcc                            >  1:3492336/1‑93 (MQ=255)
                                                         ttGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCt                           >  1:4981991/1‑93 (MQ=255)
                                                          tGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTg                          <  1:1057083/93‑1 (MQ=255)
                                                             cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt                       >  2:7290717/1‑93 (MQ=255)
                                                             cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt                       >  2:2536023/1‑93 (MQ=255)
                                                             cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt                       <  2:4940009/93‑1 (MQ=255)
                                                             cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt                       <  1:2137403/93‑1 (MQ=255)
                                                                aGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGt                    >  2:4709375/1‑93 (MQ=255)
                                                                aGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGt                    >  2:611804/1‑93 (MQ=255)
                                                                 gggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGtt                   >  1:5994533/1‑93 (MQ=255)
                                                                  ggAGAATTGCCATATGGCGGAGCGGGAAATTTCGATCTTGTGGGTCAGGGCGCCCTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTa                  <  1:1657713/93‑1 (MQ=255)
                                                                  ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTa                  <  2:3462132/93‑1 (MQ=255)
                                                                   gagaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAg                 <  1:6141653/93‑1 (MQ=255)
                                                                     gaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTa               <  2:3748900/93‑1 (MQ=255)
                                                                      aaTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTac              <  2:5419662/93‑1 (MQ=255)
                                                                         tGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTacaac           <  1:3065947/93‑1 (MQ=255)
                                                                         tGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTacaac           <  2:1347514/93‑1 (MQ=255)
                                                                         tGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTacaac           <  2:4207400/93‑1 (MQ=255)
                                                                          gCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTacaaca          <  1:6616559/93‑1 (MQ=255)
                                                                           ccATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcag         <  2:5930834/93‑1 (MQ=255)
                                                                           ccATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcag         <  1:4410141/93‑1 (MQ=255)
                                                                            cATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagc        <  2:2384880/93‑1 (MQ=255)
                                                                            cATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagc        <  2:1511545/93‑1 (MQ=255)
                                                                            cATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagc        <  1:1791581/93‑1 (MQ=255)
                                                                             atatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagca       <  1:947376/93‑1 (MQ=255)
                                                                              tatGGCGAAGCGCCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAAcagcag      >  1:511933/1‑93 (MQ=255)
                                                                                 ggCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAgcgg   <  2:1423520/93‑1 (MQ=255)
                                                                                 ggCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAgcgg   <  2:6881505/93‑1 (MQ=255)
                                                                                       ||                                                                                      
ATAATTGAAAAATTAGATCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGGTTAGTACAACAGCAGCGGC  >  NC_000913/4194490‑4194664

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