Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,732,361 T→C T232A (ACC→GCC)  yfiH ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,732,3610TC99.4% 550.1 / NA 171T232A (ACC→GCC) yfiHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (79/91);  total (80/91)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

AGATTAATGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGTCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCA  >  NC_000913/2732271‑2732450
                                                                                          |                                                                                         
aGATTAATGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCt                                                                                         <  2:3584710/93‑1 (MQ=255)
  aTTAATGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTg                                                                                       >  1:1638726/1‑93 (MQ=255)
   ttAATGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGt                                                                                      >  1:3408781/1‑93 (MQ=255)
     aaTGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTcg                                                                                    <  2:4659490/93‑1 (MQ=255)
      aTGTGATAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTcgc                                                                                   >  1:3819378/1‑93 (MQ=255)
      aTGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTcgc                                                                                   >  2:3552367/1‑93 (MQ=255)
      aTGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTcgc                                                                                   <  1:2460461/93‑1 (MQ=255)
        gtgAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcga                                                                                 <  2:904405/93‑1 (MQ=255)
        gtgAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcga                                                                                 <  2:3794147/93‑1 (MQ=255)
         tgAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcgac                                                                                <  2:616222/93‑1 (MQ=255)
         tgAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcgac                                                                                <  2:4655216/93‑1 (MQ=255)
         tgAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcgac                                                                                <  1:2814102/93‑1 (MQ=255)
          gaaaagaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcgacg                                                                               >  2:838332/1‑93 (MQ=255)
          gaaaagaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcgacg                                                                               >  2:77197/1‑93 (MQ=255)
           aaaagaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGcgacga                                                                              >  1:507918/1‑93 (MQ=255)
             aagaaaaTCACGTGTACCGGATCGTCTTTATTCTTTAGGTTATACCGGCCAAATGAACCTTGCCATCCGTCCGGTGGCCTTGTCGCGCCGATa                                                                            <  2:1883208/93‑1 (MQ=255)
             aagaaaaTCACGCGTACCGGATCGTGTTGATTCTTTAGGTTATATGAGGCAAATGAAGATTGCCATACGACCGGGGGCCTTGGCGcgacggtg                                                                            >  1:2946080/1‑90 (MQ=255)
             aagaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGGCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATa                                                                            >  1:2946079/1‑93 (MQ=255)
             aagaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATa                                                                            <  1:2467779/93‑1 (MQ=255)
              agaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATaa                                                                           <  1:3753512/93‑1 (MQ=255)
               gaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAag                                                                          >  2:1264121/1‑93 (MQ=255)
               gaaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAag                                                                          <  2:4525457/93‑1 (MQ=255)
                aaaaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAaga                                                                         <  1:4675356/93‑1 (MQ=255)
                  aaTCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAgaga                                                                       <  2:1092800/93‑1 (MQ=255)
                    tCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGagaag                                                                     <  1:328918/93‑1 (MQ=255)
                     cACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGagaaga                                                                    >  1:2754627/1‑93 (MQ=255)
                        gcgTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAg                                                                 >  1:1031509/1‑93 (MQ=255)
                           tACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGtct                                                              >  1:4430162/1‑93 (MQ=255)
                             ccGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCa                                                            <  2:3678086/93‑1 (MQ=255)
                             ccGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCa                                                            <  1:3154869/93‑1 (MQ=255)
                             ccGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCa                                                            >  2:433593/1‑93 (MQ=255)
                              cGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAt                                                           <  2:743336/93‑1 (MQ=255)
                              cGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAt                                                           >  1:2383442/1‑93 (MQ=255)
                              cGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAt                                                           <  2:4473876/93‑1 (MQ=255)
                               ggATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAtt                                                          <  2:3775538/93‑1 (MQ=255)
                                gATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAttt                                                         >  1:4238954/1‑93 (MQ=255)
                                gATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAttt                                                         >  2:382875/1‑93 (MQ=255)
                                 aTCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAtttt                                                        <  1:3348914/93‑1 (MQ=255)
                                 aTCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAtttt                                                        >  1:513413/1‑93 (MQ=255)
                                 aTCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCAtttt                                                        >  1:3178090/1‑93 (MQ=255)
                                 aTCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCATTGTCGCGACGATAAGAGAAGACAGTCTCAtttt                                                        >  1:1556388/1‑93 (MQ=255)
                                  tCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTc                                                       >  1:2984991/1‑93 (MQ=255)
                                   cGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTcc                                                      <  1:834900/93‑1 (MQ=255)
                                     tCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGt                                                    >  2:3123884/1‑93 (MQ=255)
                                     tCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGt                                                    <  2:2016854/93‑1 (MQ=255)
                                     tCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGt                                                    >  2:4413406/1‑93 (MQ=255)
                                     tCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGt                                                    >  1:2356677/1‑93 (MQ=255)
                                     tCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGt                                                    >  1:576102/1‑93 (MQ=255)
                                      cTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGta                                                   >  2:1272630/1‑93 (MQ=255)
                                        tGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGtata                                                 <  1:2914935/93‑1 (MQ=255)
                                          aTTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTTGCGATGATAAGAGAAGAAAGTCTCATTTTCCGTATATg                                               >  2:4202372/1‑93 (MQ=255)
                                           ttCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGt                                              >  1:4088191/1‑93 (MQ=255)
                                           ttCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGt                                              >  1:102010/1‑93 (MQ=255)
                                           ttCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGt                                              >  2:2200540/1‑93 (MQ=255)
                                             cTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTAc                                            <  1:4444273/93‑1 (MQ=255)
                                             cTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTAc                                            <  2:2363658/93‑1 (MQ=255)
                                             cTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGCCCTTTCGCGACCATAAAAGAAGGAAGGCTCATTTTTCGTATATGTACa                                           >  2:1216101/1‑93 (MQ=255)
                                               ttAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACaa                                          <  2:3191070/93‑1 (MQ=255)
                                                  ggTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACgg                                       >  1:3460359/1‑93 (MQ=255)
                                                  ggTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACgg                                       <  2:3270545/93‑1 (MQ=255)
                                                   gTTATATCAGCCAAATGAAACTTGCCATACGACGGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGt                                      >  2:2560717/1‑93 (MQ=255)
                                                   gTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGGCGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGt                                      >  2:1295949/1‑93 (MQ=255)
                                                   gTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGt                                      <  2:2711129/93‑1 (MQ=255)
                                                   gTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGt                                      <  1:5345/93‑1 (MQ=255)
                                                    ttATATCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTc                                     <  2:3065052/93‑1 (MQ=255)
                                                      atatCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTcgc                                   >  2:3493253/1‑93 (MQ=255)
                                                       tatCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTcgcc                                  >  2:491524/1‑93 (MQ=255)
                                                       tatCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTcgcc                                  <  1:2210039/93‑1 (MQ=255)
                                                       tatCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTcgcc                                  >  1:2261472/1‑93 (MQ=255)
                                                       aatCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTcgcc                                  <  2:1538388/92‑1 (MQ=255)
                                                        atCAGCCAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTcgccg                                 >  1:1096823/1‑93 (MQ=255)
                                                              cAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGaaa                           <  2:932254/93‑1 (MQ=255)
                                                              cAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGaaa                           >  1:4647638/1‑93 (MQ=255)
                                                              cAAATGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGaaa                           >  2:1423685/1‑93 (MQ=255)
                                                               aaaTGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGaaaa                          <  2:547703/93‑1 (MQ=255)
                                                               aaaTGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGaaaa                          <  2:2646760/93‑1 (MQ=255)
                                                               aaaTGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGaaaa                          >  2:1388294/1‑93 (MQ=255)
                                                                aaTGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAAt                         <  1:4495889/93‑1 (MQ=255)
                                                                 aTGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAAtt                        <  1:1399634/93‑1 (MQ=255)
                                                                  tGAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAAttt                       <  2:1777588/93‑1 (MQ=255)
                                                                   gAAACTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTg                      >  1:656120/1‑93 (MQ=255)
                                                                    aaaCTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGc                     <  1:1834712/93‑1 (MQ=255)
                                                                    aaaCTTGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGc                     <  1:3636749/93‑1 (MQ=255)
                                                                        ttGCCATACGACCGGTGGTCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCAGAAAATTTGCTCaa                 >  2:4040605/1‑93 (MQ=255)
                                                                         tGCCATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAac                <  2:283225/93‑1 (MQ=255)
                                                                           ccATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAacac              <  1:2479370/93‑1 (MQ=255)
                                                                            cATACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACAcc             <  2:2142551/93‑1 (MQ=255)
                                                                             aTACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACAccc            <  2:3176482/93‑1 (MQ=255)
                                                                             aTACGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACAccc            <  2:4288806/93‑1 (MQ=255)
                                                                               aCGACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCAc          >  1:1207136/1‑93 (MQ=255)
                                                                                 gACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGt        <  2:2856553/93‑1 (MQ=255)
                                                                                 gACCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGt        <  2:3918135/93‑1 (MQ=255)
                                                                                  aCCGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGtt       <  2:1529824/93‑1 (MQ=255)
                                                                                   ccGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTc      <  1:2028353/93‑1 (MQ=255)
                                                                                    cGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCg     <  2:3331877/93‑1 (MQ=255)
                                                                                    cGGTGGCCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCg     >  2:1055202/1‑93 (MQ=255)
                                                                                     ggtggcCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGc    <  2:3450523/93‑1 (MQ=255)
                                                                                      gtggcCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGcc   >  1:2126661/1‑93 (MQ=255)
                                                                                      gtggcCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGcc   >  2:2406420/1‑93 (MQ=255)
                                                                                      gtggcCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGcc   <  2:2481918/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGATTAATGTGAAAAGAAAATCACGCGTACCGGATCGTCTTGATTCTTTAGGTTATATCAGCCAAATGAAACTTGCCATACGACCGGTGGTCTTGTCGCGACGATAAGAGAAGAAAGTCTCATTTTCCGTATATGTACAACGGTCGCCGCCGAAAATTTGCTCAACACCCACGTTCGCCA  >  NC_000913/2732271‑2732450

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