Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,200,901 A→G *466W (TGA→TGG) 
T2A (ACG→GCG) 		 hydH →
hydG →		 HydH
HydG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,200,9010AG84.7% 656.9 / 113.1 287*466W (TGA→TGG) 
T2A (ACG→GCG) 		hydH
hydG		HydH
HydG
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base G (122/121);  total (142/145)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 288 total aligned reads displayed.

ACAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAAC  >  NC_000913/4200812‑4200990
                                                                                         |                                                                                         
aCAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGgcgc                                                                                        >  2:1934689/1‑93 (MQ=255)
 cAATTCGGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCa                                                                                       >  2:3139427/1‑93 (MQ=255)
 cAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCa                                                                                       >  1:3757622/1‑93 (MQ=255)
  aaTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCAc                                                                                      <  1:6334017/93‑1 (MQ=255)
   aTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACg                                                                                     >  2:252110/1‑93 (MQ=255)
   aTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACg                                                                                     >  2:4333900/1‑93 (MQ=255)
    ttCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGa                                                                                    >  1:659826/1‑93 (MQ=255)
     tCAGGTCGCAAGCCGGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGat                                                                                   <  2:6983649/93‑1 (MQ=255)
     tCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGat                                                                                   <  2:7308651/93‑1 (MQ=255)
     tCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGat                                                                                   >  1:2787720/1‑93 (MQ=255)
     tCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGat                                                                                   <  1:3227353/93‑1 (MQ=255)
      cAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGata                                                                                  <  2:2757723/93‑1 (MQ=255)
      cAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGata                                                                                  <  2:2281020/93‑1 (MQ=255)
       aGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGataa                                                                                 >  2:654579/1‑93 (MQ=255)
        ggTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataat                                                                                >  2:6006819/1‑93 (MQ=255)
        ggTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataat                                                                                >  1:7223622/1‑93 (MQ=255)
        ggTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGataat                                                                                >  2:4081594/1‑93 (MQ=255)
         gTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataata                                                                               <  1:4242647/93‑1 (MQ=255)
          tCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCAGCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAatat                                                                              >  1:5385041/1‑93 (MQ=255)
           cGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATc                                                                             <  2:2695035/93‑1 (MQ=255)
             cAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGa                                                                           <  1:6526729/93‑1 (MQ=255)
             cAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGa                                                                           <  1:1418196/93‑1 (MQ=255)
             cAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGa                                                                           >  2:1071893/1‑93 (MQ=255)
               aGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGata                                                                         >  2:468650/1‑93 (MQ=255)
                gCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGatat                                                                        <  1:6005247/93‑1 (MQ=255)
                gCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGatat                                                                        >  2:16712/1‑93 (MQ=255)
                gCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGatat                                                                        <  2:6408612/93‑1 (MQ=255)
                gCCAGGAGGGAAAAAGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGatat                                                                        >  2:5182636/1‑93 (MQ=255)
                  cAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTc                                                                      >  1:840658/1‑93 (MQ=255)
                  cAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTc                                                                      >  1:252989/1‑93 (MQ=255)
                   aggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCt                                                                     >  1:3934625/1‑93 (MQ=255)
                   aggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCt                                                                     <  1:5467180/93‑1 (MQ=255)
                    ggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtg                                                                    <  2:6712765/93‑1 (MQ=255)
                    ggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtg                                                                    >  2:750932/1‑93 (MQ=255)
                    ggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtg                                                                    >  1:4190581/1‑93 (MQ=255)
                    ggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCtg                                                                    <  2:5431009/93‑1 (MQ=255)
                        ggAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtgg                                                                >  1:1259145/1‑93 (MQ=255)
                        ggAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtgg                                                                <  2:6824583/93‑1 (MQ=255)
                        ggAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtgg                                                                >  1:1455994/1‑93 (MQ=255)
                         gAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggt                                                               >  1:4832248/1‑93 (MQ=255)
                         gAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggt                                                               <  2:5665880/93‑1 (MQ=255)
                           aaaGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggtgg                                                             >  1:4573888/1‑93 (MQ=255)
                           aaaGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggtgg                                                             >  2:5733562/1‑93 (MQ=255)
                              ggCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgatg                                                          <  2:1828433/93‑1 (MQ=255)
                                cTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgatgat                                                        <  1:1151801/93‑1 (MQ=255)
                                  cAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGAtgatga                                                      <  2:4835676/93‑1 (MQ=255)
                                  cAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGAtgatga                                                      <  2:3509164/93‑1 (MQ=255)
                                   aaCGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGAc                                                     >  2:7408772/1‑93 (MQ=255)
                                    aCGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACa                                                    <  1:7494291/93‑1 (MQ=255)
                                     cGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACAt                                                   <  2:3887714/93‑1 (MQ=255)
                                       ttCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTa                                                 <  2:4483880/93‑1 (MQ=255)
                                       ttCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTa                                                 <  1:2005843/93‑1 (MQ=255)
                                           cccTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCa                                             <  2:3096457/93‑1 (MQ=255)
                                           cccTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCa                                             <  1:7364885/93‑1 (MQ=255)
                                             ctctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACt                                           >  1:3928470/1‑93 (MQ=255)
                                             ctctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACt                                           <  2:3006275/93‑1 (MQ=255)
                                             ctctGGCTTCCGGCCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGCGGTGGATGATGACATTAGCCACt                                           >  1:5690368/1‑93 (MQ=255)
                                              tctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTg                                          >  2:1735411/1‑93 (MQ=255)
                                              tctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTg                                          <  2:7211747/93‑1 (MQ=255)
                                              tctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTg                                          >  1:1058153/1‑93 (MQ=255)
                                               ctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGc                                         <  2:3539008/93‑1 (MQ=255)
                                               ctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGc                                         >  2:5243796/1‑93 (MQ=255)
                                                tGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCa                                        >  2:2182271/1‑93 (MQ=255)
                                                   cTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTa                                     >  1:274600/1‑93 (MQ=255)
                                                    ttCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAt                                    >  1:4308618/1‑93 (MQ=255)
                                                      ccGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAttt                                  <  1:5301107/93‑1 (MQ=255)
                                                      ccGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAttt                                  <  2:1231526/93‑1 (MQ=255)
                                                      ccGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAttt                                  >  1:3480795/1‑93 (MQ=255)
                                                         gTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGc                               >  1:7441512/1‑93 (MQ=255)
                                                         gTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGc                               >  2:2916458/1‑93 (MQ=255)
                                                         gTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGc                               >  1:2321813/1‑93 (MQ=255)
                                                          tCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCa                              <  1:2576459/93‑1 (MQ=255)
                                                          tCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCa                              <  2:6271866/93‑1 (MQ=255)
                                                          tCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCa                              <  2:4885607/93‑1 (MQ=255)
                                                           cAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAg                             <  2:4678038/93‑1 (MQ=255)
                                                            aaTATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAgg                            >  2:3414687/1‑93 (MQ=255)
                                                            aaTATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAgg                            >  2:1380250/1‑93 (MQ=255)
                                                              tatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCt                          <  1:5952010/93‑1 (MQ=255)
                                                              tatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCt                          >  1:5834319/1‑93 (MQ=255)
                                                                ttACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCttt                        >  1:1924938/1‑93 (MQ=255)
                                                                  aCGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTAc                      <  1:6425538/93‑1 (MQ=255)
                                                                  aCGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTAc                      <  1:886291/93‑1 (MQ=255)
                                                                   cgcgTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACt                     >  2:3755687/1‑93 (MQ=255)
                                                                   cgcgTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACt                     <  2:7099151/93‑1 (MQ=255)
                                                                      gTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcg                  <  1:5985680/93‑1 (MQ=255)
                                                                      gTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcg                  >  2:5016314/1‑93 (MQ=255)
                                                                       tAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcgc                 <  2:3913361/93‑1 (MQ=255)
                                                                       tAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcgc                 >  2:1119646/1‑93 (MQ=255)
                                                                       tAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcgc                 >  2:5752976/1‑93 (MQ=255)
                                                                       tAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcgc                 >  2:5897489/1‑93 (MQ=255)
                                                                       tAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcgc                 <  1:3333334/93‑1 (MQ=255)
                                                                         aGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCgg               >  2:4059968/1‑93 (MQ=255)
                                                                             cccACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgg           <  1:4944888/93‑1 (MQ=255)
                                                                              ccACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTggg          <  2:5622453/93‑1 (MQ=255)
                                                                               cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         >  1:5214819/1‑93 (MQ=255)
                                                                               cacaAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         <  1:4609233/93‑1 (MQ=255)
                                                                                 caAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCt       >  1:5144485/1‑93 (MQ=255)
                                                                                    ggATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCtata    >  2:3236597/1‑93 (MQ=255)
                                                                                     gATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATaa   <  2:1595610/93‑1 (MQ=255)
                                                                                     gATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATaa   <  2:260921/93‑1 (MQ=255)
                                                                                         |                                                                                         
ACAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAAC  >  NC_000913/4200812‑4200990

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